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I.  INTRODUCTION 


A.  BACKGROUND 

Over  5,  000  years  ago,  evolution  of  the  post  and  lintel 
structures  of  the  stone  age  gave  rise  to  the  arch.  Highly 
regarded  for  its  graceful  shape  and  design  suitability,  the 
simple  arch  structure  has  been  applied  to  engineering  and 
architectural  designs  ever  since.  The  ancient  Roman  Coliseum 
and  aqueducts,  great  cathedrals  of  the  Middle  Ages,  and 
railway  bridges  of  modern  history  are  just  a  few  of  the  many 
examples  of  structures  comprised  of  arches  standing  today 
(Figure  1.1).  Throughout  its  history,  engineers  and  architects 
have  labored  to  improve  the  design  of  the  arch  in  order  to 
enhance  the  overall  design  structure.  This  desire  for 
perfection  has  led  engineers  to  devise  a  rational,  directed 
design  procedure  and  hence,  the  concept  of  optimization  was 
created. 

The  advent  of  the  computer  era  has  lead  to  20  years  of 
extensive  development  in  the  use  of  numerical  optimization 
techniques.  These  techniques  offer  a  logical  approach  to 
design  decisions  where  intuition  and  experience  previously 
prevailed.  Coupled  with  trends  toward  material  and  cost 
efficiency,  numerical  optimization  has  prompted  considerable 
research  in  the  field  of  automated  design  [Ref.  1].  As 
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Figure  1 . 1  (a)  The  Coliseun 

(b)  Firth  of  Forth  Railway  Bridge 

(c)  Isernia,  Italy,  railway  bridge 

Photos  by  Mrs.  P.  Menzlos  and  COR  D.C.  Hamer 
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a  step  in  design  optimization  of  structures,  the  arch  has  been 
the  subject  of  numerous  optimality  studies  to  enhance 
applicability  in  engineering  and  architectural  designs. 

One  such  study  was  performed  by  Farshad  in  1976 
[Ref.  2]  .  Using  calculus  of  variations,  he  derived 
optimality  conditions  for  nonlinear  partial  differential 
equations  for  hinged-hinged  arches.  The  total  potential  energy 
of  the  system,  augmented  with  several  objective  functions  via 
Lagrange  multipliers,  was  minimized  with  respect  to  design  and 
state  variables  to  achieve  equilibrium  and  optimality.  The 
nonlinear  systems  of  equations  for  optimal  thrust,  minimum 
length  of  the  arch,  and  minimum  volume  were  presented  but  not 
solved. 

In  1980,  Rozvany  et  al.  [Ref.  3]  used  the  Prager- 
Shield  criteria  to  optimize  statically  determinate  arches.  His 
'arch'  consisted  of  two  inclined  funicular  frame  beams 
ridgedly  interconnected  with  a  concentrated  load  applied  at 
the  joints.  In  the  optimal  'arch'  only  bending  or  axial  forces 
develop  depending  on  the  ratio  of  4L/D,  where  L  is  the  span  of 
the  structure  and  D  is  the  depth  of  the  cross  section.  Ratios 
greater  than  eight  to  one  produced  axial  forces  only  and  the 
optimal  shape  has  a  height  of  half  the  span.  Ratios  smaller 
than  eight  to  one  develop  only  bending  and  the  optimum 
structure  is  a  straight  beam.  In  each  case,  the  width  of  the 
beam  segments  for  the  optimal  'arch'  varied  linearly  from  the 
hinged  support  to  the  axis  of  symmetry. 
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That  same  year,  Lipson  et  al.  [Ref.  4]  used  the 
'complex'  method  to  optimize  parabolic  arches  subject  to 
uniform  loading.  His  'arch'  was  comprised  of  equal  length 
straight  beam  sections  of  thin  walled  rectangular  tubes. 
Maintaining  constant  depth  and  width  for  each  segment,  the 
vertical  and  horizontal  wall  thicknesses  determined  the  arch 
shape  which  was  optimized  for  minimum  total  weight.  An  arch 
with  a  rise  of  0.342  times  the  span  length  proved  to  be  the 
optimum. 

In  1988,  Ang  et  al.  [Ref.  5]  solved  the  arch 
optimization  problem  by  parametricing  the  unspecified  arch 
axis  using  spline  functions  and  employing  a  smoothing  function 
to  approximate  the  non-smooth  objective  function.  The  'arch' 
was  considered  to  be  a  'plastic'  design  of  rectangular  cross 
section  subject  to  bending  and  axial  compression.  Three  types 
of  boundary  conditions  were  imposed,  simply  supported-simply 
supported,  clamped-clamped,  and  simply  supported-clamped.  The 
optimum  shape  of  the  arch  is  claimed  to  be  a  parabola  with  a 
rise  of  0.433  times  the  span  length.  Apparently,  there  is  some 
disagreement  between  these  results  and  those  previously  noted. 

In  addition  to  arch  optimization  studies.  Ding  and  Esping 
[Ref.  6]  solved  the  minimum  weight  design  problem  for 
frame  structures  when  stress  and  displacement  constraints  are 
considered.  Using  dual  numerical  methods,  seven  cross- 
sectional  shapes  were  treated  by  approximating  the  stresses 
with  pseudo  and  virtual  load  techniques.  Results  were 
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presented  for  a  beam  clamped  at  both  ends,  a  portal  frame,  a 
2X5  grillage,  and  a  helicopter  tail  boom  structure.  Although 
Ding  and  Esping's  investigation  does  not  specifically  solve 
for  arch  structures,  the  approximations  used  are  completely 
detailed  with  convincing  results. 

In  December  of  1990,  Charles  Scott  McDavid  of  the  Naval 
Postgraduate  School  presented  his  thesis,  "Weight  Optimum  Arch 
Structures, "  which  optimized  circular  arches  subject  to 
various  loadings  and  end  conditions.  Specifically,  he 
optimized  arches  segmented  into  rectangular  boxes  that  varied 
in  width  only.  Through  his  research  he  concluded  that  a 
bar/beam  element  model  is  a  viable  technique  for  the 
approximation  of  arch  structures,  and  that  an  arch  structure 
that  is  more  statically  indeterminate  is  more  efficient  under 
identical  loading.  Additionally,  he  proposed  possibilities  for 
future  research  which  includes  varying  both  the  height  and 
width  dimension,  the  major  thrust  of  this  investigation. 

B .  PROBLEM  DEFINITION 

In  order  to  provide  an  in  depth  study,  each  of  the  cited 
investigations  began  with  a  problem  definition  and  specific 
assumptions  about  the  type  of  arch  to  be  considered.  For  this 
investigation,  the  arch  is  defined  as  a  structure  of  constant 
curvature  (i.e.,  circular  arches)  which  when  supported  at  both 
ends  and  loaded  laterally  develops  perpendicular  reactions. 
This  is  intended  to  eliminate  thick  walled  curved  beams  and 
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straight  beams  which  develop  virtually  no  perpendicular 
reactions  when  loaded  laterally.  Additionally,  the  cross- 
section  dimensions  are  small  relative  to  the  radius  of 
curvature  and  therefore  the  centroidal  and  neutral  axes  are 
assumed  to  coincide.  Without  the  thin  depth  assumption, 
complications  arise  in  the  calculations  of  the  displacements 
and  the  slopes  because  the  arch  no  longer  behaves  as  predicted 
by  the  beam  equilibrium  equation: 

(EIv"  )"  =Py(s)  (1.1) 

and  the  bar  equilibrium  equation: 

{AEu'  )'  =-P^{s)  (1.2) 

where  the  prime  superscript  notation  denotes  differentiation 
with  respect  to  the  independent  variable,  s,  and 
E  =  Young's  Modulus 

I  =  Cross-sectional  Moment  of  Inertia 

V  =  Lateral  Displacement 

Py  =  Lateral  Loading 

A  =  Cross-sectional  Area 

u  =  Axial  Displacement 

Px  =  Axial  Loading 

s  =  the  Independent  Variables 

In  order  to  facilitate  the  development  of  a  finite  element 

code  to  approximate  the  local  displacements,  the  arch  is 

approximated  by  a  series  of  straight  segments.  From  the  local 

displacements,  the  virtual  load  techniques,  as  described  in 

the  Ding  and  Esping  paper,  are  applied  to  determine  the 

internal  psuedostresses .  Once  the  stress  distribution  is 

determined,  the  arch  volume  is  minimized  to  a  structure  that 
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maintains  the  developed  stresses  below  the  predefined  maximum 
allowable  stress. 

The  thrust  of  this  investigation  is  to  minimize  the  total 
weight  of  a  linearly  elastic,  isotropic,  and  homogeneous  arch 
under  a  variety  of  loadings  and  end  conditions.  Optimization 
in  this  investigation  refers  to  the  variance  of  the  cross- 
sectional  dimensions  (that  is,  the  design  variables)  to  obtain 
optimum  least  weight  structures.  Design  Optimization  Tool 
(DOT)  software  [Ref.  7]  is  used  to  perform  the 
optimization  subject  to  prescribed  constraints  on  the  design 
variables  as  well  as  on  the  stress  limitations.  The  objective 
is  to  minimize  the  total  volume  of  the  arch  while  maintaining 
stresses  below  the  yield  strength  of  the  arch  material.  The 
intent  of  this  study  is  to  provide  direction  and  guidance  on 
which  further  research  for  weight  optimization  may  be 
developed. 
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II.  PROBLEM  FORMaiATION 


A.  PROBLEM  STATEMENT  AMD  ASSUMPTIONS 

As  noted  in  the  introduction,  the  purpose  of  this 
investigation  is  to  optimize  arch  structures  to  form  a 
foundation  upon  which  further  research  can  be  based.  These 
arch  structures,  subject  to  specified  loadings  and  end 
conditions,  vary  in  cross  sectional  geometry  to  minimize  the 
weight.  In  order  to  limit  the  scope  of  this  study, 
approximations  and  specific  assumptions  are  made  as  follows: 

•  The  arch  maintains  a  constant  radius  of  curvature. 

•  The  arch  is  approximated  by  a  series  of  straight  segments 
of  a  solid  rectangular  cross  sectional  geometry. 

•  Cross  section  design  is  restricted  to  ensure  the 
applicability  of  beam  and  bar  equilibrium  equations  (1.1 
and  1.2). 

•  To  prevent  failure  the  internal  stresses  developed  due  to 
the  loading  must  not  exceed  the  yield  strength  of  the 
material . 

•  The  arch  structure  is  composed  of  a  linearly  elastic, 
isotropic,  homogeneous  material. 

To  begin  the  design  optimization  process,  the  arch 
structure  is  approximated  by  contiguous  straight  line 
segments.  Each  segment  is  modeled  by  a  bar-beam  structure 
connects  to  the  adjacent  segment  at  a  point  defined  as  the 
nodal  point.  At  each  nodal  point,  the  cross  section  base  and 
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height  dimensions  are  selected  as  the  design  variables.  From 
this  model,  the  optimization  problem  can  be  formulated  into 
objective  and  constraint  functions  which  are  functions  of 
these  design  variables, 

B.  MUTBEMZITICAL  MODEL 

Due  to  the  complex  nature  of  this  problem,  the  constant 
radius  arch  structure  is  modeled  by  a  series  of  straight 
contiguous  elements  where  the  arch  radius  of  curvature,  R,  and 
the  number  of  elements  used  to  approximate  the  arch,  NEL,  is 
specified.  (Figure  2.1)  For  simplicity,  the  length  of  each 
element  is  constant  such  that: 

L  =  6R/  NEL 

where  6  represents  the  subtended  arc  of  the  arch. 


9 


At  each  nodal  point,  there  exists  a  base  and  height 
dimension  such  that  the  cross  sectional  dimensions  from  one 
element  to  the  adjacent  element  maintains  smooth  piecewise 
continuity.  (Figure  2.2)  The  resultant  element  shape  is  that 
of  a  three  dimensional  trapezoid  whereby  the  volume  is 
calculated  by  multiplying  the  average  base  and  height  with  the 
length  of  the  element.  In  mathematical  terms,  the  volume  of 
the  ith  element  is  calculated  as  follows: 


Vol ume  (  i  )  =  (i)  *  i  )  *  1- 

(2.1) 

where 

(B  (i)  +  B  (i+l))/2 

(2.2) 

U 

= 

(H  (i)  +  H  a+l))/2 

(2.3) 

B 

= 

the  Nodal  Base  Dimension 

H 

= 

the  Nodal  Height  Dimension 

L 

the  Element  Length 

i 

=r 

the  ith  Element 

Figur*  2.2  Arch  El«m«nts 
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Defined  in  the  problem  statement,  the  optimal  arch  is 
achieved  by  varying  the  cross  sectional  dimensions,  the  base 
and  height,  in  order  to  minimize  the  weight.  Thus  the  nodal 
base  and  height  dimensions  are  the  design  variables  for  which 
the  objective  function  is  defined. 

C.  OPTIMIZATION  PROBLEM 

The  objective  of  this  study  is  to  minimize  the  weight  of 
an  arch  structure  while  maintaining  a  stress  distribution 
which  does  not  exceed  the  yield  strength  of  the  material. 
Additionally,  other  constraints  on  the  design  variables  are 
imposed.  Since  the  arch  is  composed  of  a  homogenous  material, 
the  weight  of  the  arch  is  directly  proportional  to  the  volume 
of  the  arch.  Thus,  the  objective  of  this  investigation  is 
satisfied  by  minimizing  the  total  arch  volume.  The  total  arch 
volume,  Vtot,  is  the  sum  of  the  elemental  volume,  v(i).  Thus 
in  mathematical  form,  the  objective  function  is  as  follows: 

V  ( i ) 

where  the  elemental  volumes,  v(i),  calculated  by  Equation 
(2.1),  is  summed  for  all  elements  to  compute  the  total  arch 
volume . 

In  keeping  with  the  assumptions  made  in  the  problem 
statement,  the  objective  function  is  constrained  in  order  to 
impose  practical  and  important  physical  restriction  on  the 


Object i  ve  =  MIN  ( 


(2.4) 
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problem.  Properly  defined,  the  constraints  are  used  to  avoid 
undesirable  behavior  such  as  yielding,  to  ensure  validity  of 
the  governing  equilibrium  equations,  and  to  provide  a 
realistic  design.  For  this  study,  the  constraints  fall  into 
three  categories,  strength  criteria,  geometric  limitations, 
and  side  constraints. 

First,  for  specified  loadings  and  end  conditions,  the 
optimized  arch  must  not  'fail  by  yielding.'  Assuming  the  arch 
material  to  be  linearly  elastic,  the  applied  loading  must  not 
cause  the  structure  to  exceed  the  elastic  limit  of  the 
selected  material.  Therefore,  the  internal  stresses  developed 
must  remain  below  the  yield  strength  of  the  material. 
Mathematically,  the  strength  criteria  is  as  follows: 

O  (i)  ^Sy 

or  in  normalized  form; 

(  o  (i)  /  Sy)  -1.0  <  0.0  (2.5) 

where  a(i)  is  the  maximum  stress  developed  at  the  ith  nodal 
point  of  the  arch  and  Sy  is  the  yield  strength  of  the  arch 
material  selected  by  the  designer.  Unfortunately,  the  stress 
distribution,  in  terms  of  the  design  variables  is  not  readily 
available.  However,  using  the  beam  and  bar  equilibrium 
equations  (1.1  and  1.2),  a  finite  element  scheme  based  on  the 
model  can  be  developed  to  determine  the  arch's  displacements 
and  slopes  due  to  a  given  loading.  Knowing  how  the 
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displacements  and  slopes  change  throughout  the  arch,  the 
stresses  at  the  nodal  points  can  be  calculated. 

Secondly,  limits  must  be  imposed  on  the  cross  sectional 
geometry  in  order  to  ensure  applicability  of  the  bar  and  beam 
equilibrium  equations.  Limiting  the  cross  section  base  and 
height  dimensions  relative  to  one  another  prevents  the 
structure  from  becoming  either  a  shell-like  or  deep  curved 
beam  structure.  To  maintain  the  geometry  of  the  arch,  the 
following  conditions  are  imposed: 

B  (i)  -3.0  »  W  (i)  <  0.0  (2.6) 

and 

H(i)  -10.0  *B  (j)  <0.0  (2.7) 

Finally,  the  side  constraints  are  imposed  to  ensure  a 
realistic  solution.  The  arch  is  a  physical  object  that  must 
have  a  realistic  finite  cross  sectional  area;  however,  these 
section  dimensions  must  also  remain  small  relative  to  the 
radius  of  curvature  by  definition  of  the  arch.  Thus,  the  side 
constraints  for  the  base  and  height  dimensions  are  as  follows: 

0.03in.  ^B(i)  ^e.Oin.  (2.8) 

0.03in.  <B(i)  <6.0in.  (2.9) 

In  the  future,  additional  constraints  should  be  considered 
such  as  global  buckling  and  local  crippling. 
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III.  OPTIMIZATION  ANALYSIS 


To  perform  the  computer  optimization,  the  Design 
Optimizatijn  Tools  (DOT)  software  package  is  used  due  to  its 
availability,  user  friendliness,  and  reputation.  DOT,  a 
FORTRAN  77  optimization  software  package  available  from  VMA 
Engineering,  uses  numerical  search  methods  to  seek  a  minimum 
value  of  one  function,  the  objective,  subject  to  the  limits  of 
others,  the  constraints  [Ref.  7]  .  DOT  has  two  methods  for 
iteratively  solving  constrained  optimization  problems,  the 
Modified  Method  of  Feasible  Directions  and  the  Sequential 
Linear  Programming  Method. 

A.  MODIFIED  METHOD  OF  FEASIBLE  DIRECTIONS 

Modified  Method  of  Feasible  Directions  is  a  numerical 
method  that  deals  directly  with  nonlinear  problems.  For  this 
method,  a  search  direction  vector,  is  first  found.  The 

design  point  is  then  moved  in  this  direction  to  update  the 
design  variable  vector,  X,  according  to  the  equation: 

X  +a’5  (3.1) 

— q  —q-1  — q 

where  the  scaler  quantity  a*  defines  the  distance  moved  in  the 
^  direction,  and  q  represents  the  iteration  number. 
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For  an  initial  design,  say  JC/  the  design  is  moved  in  the 
direction  of  the  steepest  descent  until  a  constraint  is 
encountered . 


Figure  3.1  Modified  Method  of  Feasible  Directions:  The 
Search  Direction 

Having  encountered  the  constraint  boundary,  a  new  search 
direction  is  found  by  solving  the  subproblem: 

Maximize ; 

MAXIMIZE:  £' X  (3.2) 
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Subject  to: 


Ax^O  (3.3) 

O'*) 


where 


Z  = 


(3.5) 


£  = 


(3.6) 


"V^g,  (x) 
V’’g2  ( x) 

V^g^  ( X) 
V^F(x) 


(3.7) 


The  search  direction,  S_,  will  follow  the  constraint  yet  allow 
the  design  to  leave  a  constraint  boundary  if  the  objective 
will  reduce  farther.  In  general,  the  form  for  inequality 
constraint  problems  is: 
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Maximize : 


-Vf(x) •  S  (3.8) 

Subject  to: 

Vg^(x)  •  s<l  jej  (3.9) 

5*5<1  (3.10) 

When  the  search  direction  is  away  from  a  currently  active 
constraint  and  the  scaler  product  of  the  gradient  of  each 
critical  constraint  with  the  ^  vector  is  less  than  zero,  the 
constraint  is  omitted  from  the  set  of  active  constraints.  If 
^  is  the  null  vector  or  numerically  small,  the  optimization 
process  is  terminated  because  the  Kuhn-Tucker  conditions  for 
optimality  have  been  met. 

B.  SEQUENTIAL  LINEAR  PROGRAMMING 

The  second  numerical  method.  Sequential  Linear  Programming 
(SLP) ,  linearizes  nonlinear  objective  and  constraint  functions 
and  then  obtains  a  solution  using  linear  programming  methods. 
Once  the  approximate  solution  is  found,  the  functions  are 
linearized  about  the  new  design  point  and  the  a  linear 
programming  problem  approximated  and  solved.  By  repeatedly 
linearizing  and  solving  the  resulting  problem,  a  precise 
solution  is  achieved. 

In  general  format,  the  nonlinear  functions  are  linearized 
via  a  first-order  Taylor  series  expansion  as  follows: 
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Minimize : 


F(x)  -F(X  )  +VF(X  )  •  5X  (3.11) 

Subject  to: 

gjU)  -gjU^)  •  8x<0  j  =  l,/n  (3.12) 

where 

8X  =  X-X^  (3.13) 

and  the  zero  subscript  identifies  the  point  about  which  this 
Taylor  series  expansion  is  performed.  At  the  initial  design, 
X^,  the  objective  and  constraints  are  linearized  to  give 
straight  line  representations  of  the  functions. 

Typically,  this  method  converges  to  the  optimum  solution 
with  fewer  iterations  than  the  previous  method  mentioned. 
However,  as  seen  in  Figure  3.2,  the  optimum  of  the 
approximated  linear  problem  is  infeasible  (i.e.,  a  design  that 
violates  some  or  all  of  the  constraints) .  Additionally, 
certain  linearizations  produce  unbounded  linear  problems. 
However,  imposing  move  limits  on  the  linear  approximation 
helps  ensure  that  the  optimum  will  eventually  be  reached. 

C.  DOT  PROGRAM  PARAMETERS 

For  both  numerical  methods,  there  are  several  parameters 
that  can  be  adjusted  within  DOT  in  order  to  'fine  tune'  the 
program  for  a  specific  problem.  Fine  tuning  is  a  process  in 
which  the  program  parameters  are  internally  adjusted  to 
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optimize  the  optimizer  performance.  With  proper  tuning,  the 
optimization  process  can  be  designed  to  remain  within 
specified  tolerances  and  operate  more  efficiently.  A  complete 
listing  of  all  the  DOT  parameters  is  contained  in  Appendix  A. 
However,  for  the  purpose  of  this  investigation,  only  the 
constraint  boundaries,  auto  scaling,  and  termination  tolerance 
parameters  were  tuned  to  enhance  optimization  performance. 
[Ref.  7] 

First,  for  constrained  optimization,  the  constraint 
boundary  must  be  established.  Mathematically  defined,  the 
constraint  is  considered  active  if  its  numerical  value  is 
between  the  value  of  CT  and  CTMIN,  and  violated  if  its 
numerical  value  is  greater  than  CTMIN.  By  using  a  narrow  band 
to  approximate  the  constraint  function,  the  optimizer  is  less 
likely  to  exceed  convergence  criteria  without  achieving  an 
optimal  design.  In  the  realm  of  design,  CTMIN  is  of  particular 
concern.  Principally,  it  is  a  small  positive  number  that 
controls  how  far  the  design  can  deviate  from  the  constraint 
boundaries  and  still  be  considered  a  feasible  design.  In 
theory,  CTMIN  can  be  reduce  to  zero  to  avoid  any  constraint 
violations,  however,  it  is  not  practical  due  to  the  large 
number  of  iterations  and  computer  expense  required. 

In  addition,  it  is  normally  considered  good  engineering 
practice  to  normalize  design  variables  and  nondimensionalize 
basic  parameters  [Ref.  1].  For  optimization,  variables  are 
scaled  to  affect  normalizing  by  evaluating  the  diagonals  of 
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the  Hessian  matrix  of  the  objective  and  constraint  functions. 
As  the  optimization  proceeds,  reevaluation  is  sometimes 
necessary  to  rescale  the  variables.  The  DOT  parameter  ISCAL 
may  be  selected  to  rescale  the  design  variables  over  an 
interval  or  eliminate  the  scaling  function  all  together. 
Unfortunately,  the  DOT  manual  indicates  that  there  is  no 
established  theory  for  scaling.  Scaling  is  therefore  a 
function  of  trial  and  error. 

Last,  the  termination  criteria  also  has  a  major  effect  on 
the  efficiency  and  reliability  of  the  optimization  process. 
Termination  criteria  is  established  so  that  the  design  process 
is  stopped  when  the  number  of  iterations  exceeds  a  specified 
limit.  DOT  parameters  ITMAX  and  JTMAX  specify  the  maximum 
number  of  iterations  allowed  for  the  Modified  Method  of 
Feasible  Directions  and  the  SLP  method  respectively.  This 
ensures  that  the  program  will  not  iterate  indefinitely. 
Furthermore,  the  progress  of  the  optimization  is  checked  for 
convergence.  Design  convergence  is  achieved  when  the  change  in 
the  value  of  the  objective  function  from  one  iteration  to  the 
next  approaches  zero.  The  DABOBJ  parameter  is  a  specified 
tolerance  for  which  the  maximum  absolute  change  in  the 
objective  function  between  iteration  is  numerically  small. 
Additionally,  ITRMOR  and  ITRMST  are  parameters  which  specify 
the  number  of  consecutive  iterations  for  which  the  design 
change  is  less  than  DABOBJ  for  Modified  Method  of  Feasible 
Directions  and  the  SLP  method  respectively. 
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IV.  STRESS  ANALYSIS 


The  objective  of  this  investigation  is  to  minimize  the 
total  weight  (volume)  of  a  load  bearing  arch  subject  to 
specified  constants.  To  obtain  an  optimal  structure,  DOT  is 
interfaced  with  an  analysis  program  which  computes  the  values 
of  the  objective  and  constraint  functions  in  terms  of  the 
design  variables,  specifically  the  cross  sectional  dimensions. 
Since  the  strength  constraint  requires  that  the  stresses  at 
any  point  do  not  exceed  the  yield  strength  of  the  arch 
material,  the  stress  distribution  over  the  domain  of  the  arch 
must  be  known.  However,  as  indicated  in  Chapter  II,  the  stress 
distribution  is  not  readily  available  in  terms  of  the  cross 
sectional  dimensions.  Therefore  the  following  stress 
development  is  pursued  for  optimization. 

A.  STRESS  DEVELOPMENT 

For  this  study,  the  strength  constraint  requires  that  the 
applied  load  will  not  cause  the  arch  to  fail  by  yielding. 
Therefore,  the  internal  stresses  developed  must  remain  below 
the  yield  strength  of  the  material.  For  this  study,  the 
stresses  considered  are  composed  of  normal  stresses  due  to 
bending  moments  and  axial  forces  where  the  total  normal  stress 
is  the  algebraic  sum  of  these  components  expressed  as  follows: 
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total  normal  stress 

the  normal  stress  due  to  bending 

the  normal  stress  due  to  axial  force 
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that  these  stresses  are  negligible.  (See  Appendix  B  for  the 
complete  justification  for  Shear  stress  omission.) 

To  compute  the  two  normal  stress  components,  the  arch  is 
sectioned  and  approximated  by  straight  frame  elements.  Thus, 
the  stresses  can  be  determined  for  each  element  endpoint  (or 
nodal  point)  in  order  to  establish  the  stress  distribution. 
Each  element  is  considered  to  behave  as  both  a  tapered  beam, 
to  calculate  the  stresses  due  to  bending,  and  a  tapered  bar, 
to  calculate  stresses  due  to  axial  forces. 

First,  for  a  straight  beam  segment,  the  maximum  normal 
stress  due  to  bending,  hereafter  referred  to  as  bending 
stresses,  is  defined  by  the  following  equation: 


where  c  is  the  distance  from  the  neutral  axis  to  the  point 
furthest  from  the  neutral  axis.  The  moment,  M,  at  a  section  is 
calculated  by: 

M  =  EIv"  (4.3) 

resulting  from  the  beam  equilibrium  equation  (1.1).  With 
substitution  and  simplification.  Equation  (4.2)  becomes: 

Q^^ECv"  (4.4) 

In  the  same  manner,  the  normal  stress  due  to  axial 
behavior  is  determined.  For  a  bar  element,  the  normal  stress 
due  to  axial  forces,  hereafter  referred  to  as  axial  stresses, 
is  defined  by  the  equation: 
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a 


(4.5) 


F 
A 

where  A  is  the  cross  section  area  and  the  axial  force,  F,  is 
calculated  by: 

F-AEu'  (4.6) 

resulting  from  the  bar  equilibrium  equation  (1.2).  Again, 
substituting  and  simplifying.  Equation  (4.5)  becomes: 

0,  =  £u'  (4.7) 

Final  substitution  into  Equation  (4.1)  results  in  an 
equation  for  total  normal  stress  as  follows: 

E(cv"  +  v')  (4.8) 

where  Young's  Modulus  of  elasticity,  E,  is  a  function  of 
material  selection,  the  distance  from  the  neutral  axis  to 
extreme  fiber,  c,  is  a  function  of  cross  section  height,  and 
u'  and  v"  are  the  first  and  second  derivatives  of  axial  and 
lateral  displacements  respectively.  Using  the  Galerkin  Finite 
Element  Method,  approximate  values  for  the  axial  and  lateral 
displacements  can  be  determined  at  element  endpoints.  From 
these  values,  the  stress  distribution  is  computed  and  the 
optimization  process  can  proceed. 

B.  THE  FINITE  ELEMENT  BEAM  EQUATION  DEVELOPMENT 

The  Galer)cin  Finite  Element  Method  (FEM)  is  an 
approximation  method  which  transforms  a  linear  differential 
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equation  into  a  system  of  linear  algebraic  equations.  Using 
the  beam  equilibrium  equation  (1.1),  approximate  lateral 
displacements  for  the  arch  can  be  determined  at  the  system 
nodal  points.  For  this  method,  a  family  of  hermite  cubic  shape 
function  which  possess  the  Kronecker  Delta  property,  are 
introduced  in  order  to  maintain  the  necessary  function  and 
slope  continuity  for  the  fourth  order  beam  equation.  An 
approximate  solution,  for  displacement,  V,  is  formed  as 
follows : 

(4.9) 

where  v  is  the  exact  solution  of  the  beam  equation  in 
continuous  space,  is  the  approximate  solution  in  discrete 
space,  Qj  is  the  transpose  of  a  column  vector  of  the  cubic 
shape  functions,  and  V  is  the  vector  of  lateral  displacements 
and  slopes. 

After  the  approximation  is  formulated,  the  next  step  in 
the  Galerkin  method  is  to  form  the  residual,  R,  in  the 
following  format: 

R  =  S£(^)  -Py(s)  (4.10) 

where  py  is  the  lateral  excitation  force  and  2  denotes  the 
differential  operator  which  in  the  case  of  the  beam 
equilibrium  equation  is  defined  by; 

Sf(v)  =  [EI(v")  ]"  (4.11) 
With  substitution,  the  residual  becomes; 
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From  the  residual,  the  Galerkin  Equations  are  formed: 


r£(i?)ds=£  (4.13) 

D 

where  0  is  the  null  vector.  Further  substitution  for  R  into 
the  Galerkin  vector  equation  results  in  : 


]''ds  -  j^yis)  ds  =£  (4.14) 

D  C 


To  solve  the  Galerkin  Equation,  integration  by  parts  is 
performed  twice  which  yields; 


^  [  El  (J2’V)  -gfEI  "  I  3 

{Q^v)"ds-^QPy(s)  ds  =  0 


(4.15) 


where  |b  denotes  evaluation  of  these  vectors  at  the  boundary 
points  of  the  structure.  Recognizing  that  the  lateral 
displacement  and  slope  vector  is  constant.  Equation  4.15  is 
rewritten  as; 

2  [  El  (£")  "  ] '  vl  3  -^EI  {QJ)  "v\  ^ 

'El  {Qj)"dsv 


^QPyis)  ds=£ 


From  the  beam  equilibrium  Equation  (1.1),  the  shear,  V,  is 
defined  by: 
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V=EIv‘ 


(4.17) 


/// 


and  Moment,  M,  by: 


M-EIv" 


(4.18) 


Thus,  the  boundary  term  load  vectors  are  defined  by: 

V  =  £[EH£'^)"  yvl  g  (4.19a) 

and 


M  =  B  (4.19b) 

Additionally,  for  convenience  a  system  stiffness  Matrix,  K®, 
is  defined  by: 


=  j^'EI  (2'^)"ds  (4.19c) 

and  a  system  Force  vector,  by: 

F'==  (ip^(s)ds  (4.19d) 

D 

Substitution  of  Equations  (4.19  a  through  d)  into  Equation 
(4.16)  results  in  the  following  system  of  linear  algebraic 
equations : 

Zl  b""I  (4.20) 

Further  simplification  is  possible  by  defining  F^*  as  the  load 
vector  of  internal  and  external  applied  lateral  loads  by: 

=  fb  +  M|  ^ - V|  ^  (4.21) 

Thus,  Equation  (4.20)  reduces  to: 
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(4.22) 


K^v  =  F^ 

where  the  global  or  system  bending  stiffness  matrix,  K®,  is 
constructed  from  the  union  of  all  the  elemental  bending 
stiffness  matrices  and  the  global  bending  force  vector, 
F*,  is  constructed  from  the  union  of  all  the  elemental  bending 
force  vectors,  . 

C.  THE  FINITE  ELEMENT  BAR  EQUATION  DEVELOPMENT 

In  a  similar  manner  to  the  beam  equation,  the  Galerkin 
Finite  Element  Method  is  applied  to  the  bar  equilibrium 
equation  (1.2)  to  approximate  the  axial  displacements  at  the 
endpoints  of  a  bar  element.  However,  the  bar  equation  is  only 
a  second  order  linear  differential  equation.  Therefore,  a 
family  of  linear  shape  functions  which  posses  the  Kronecker 
Delta  property,  are  used  in  order  to  maintain  the  necessary 
function  continuity  only.  An  approximate  solution,  0,  for 
axial  displacement,  u,  is  formed  as  follows: 

u  =  0  =  G"jj  (4 . 23) 

where  u  is  the  exact  solution  of  the  bar  equation  in 
continuous  space,  Q  is  the  approximate  solution  in  discrete 
space,  G’'  is  the  transpose  of  a  column  vector  of  the  linear 
shape  functions,  and  u  is  the  vector  of  axial  displacements. 

After  the  approximation  is  formulated,  the  next  step  in 
the  Galerkin  method  is  to  form  the  residual,  R,  in  the 
following  format: 
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R=i£{a)  *pjs) 


(4.24) 


where  is  the  axial  excitation  force  and  2  denotes  the 
differential  operator  which  in  the  case  of  the  bar  equilibrium 
equation  is  defined  by: 

«g(u)  =  [A£{u)']'  (4.25) 

With  substitution,  the  residual  becomes: 

K=  tA£(G’'u)']' +p^(s)  (4.26) 

From  the  residual,  the  Galerkin  Equation  is  formed: 


Jg  (R)  ds=^ 

D 


(4.27) 


where  £  represents  the  null  vector.  Further  substitution  into 
the  residual  equation  results  in  : 

rGtAE(G''u)']'ds+ r^^(s)  ds=£  (4.28) 

D  D 

Unlike  the  beam  equation  development,  only  single 
integration  by  parts  is  performed  to  solve  the  Galerkin 
Equation.  This  results  in: 

A£G(G’'u)'I  B-fe  [AE(G’'u)']  ds+jGp^(s)  =£  (4.29) 

D  D 

where  represents  evaluation  at  the  boundaries  of  the 

structure.  Recognizing  that  the  axial  displacement  vector  is 
constant.  Equation  (4.29)  is  rewritten  as: 
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From  the  bar  equilibrium  Equation  (1.2),  the  axial  force,  F, 
is  defined  by: 

F=AEu'  (4.31) 

Thus,  the  boundary  term  load  vectors  are  defined  by: 

P=AEG(G’')'u|  B  (4.32a) 

Additionally,  for  convenience  a  system  stiffness  Matrix,  K^, 
is  defined  by: 

[A£(G’"]  ds  (4.32b) 

D 

and  a  system  Force  vector,  F^,  by: 

£^=JjQ3*(s)  (4.32c) 

Substitution  of  Equations  (4.32  a  through  d)  into  equation 
(4.30)  results  in  the  following  system  of  linear  algebraic 
equations : 

P-K^u*F‘  =  0_  (4.33) 

Further  simplification  is  possible  by  defining  as  the  load 
vector  of  internal  and  external  applied  lateral  loads  by; 

+  p  (4.34) 

Thus,  Equation  (4.33)  reduces  to: 
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(4.35) 


K'^u=F* 

here  the  global  or  system  axial  stiffness  matrix,  K^,  is 
constructed  from  the  union  of  all  the  elemental  axial 
stiffness  matrices  }c‘^  and  the  global  axial  force  vector,  F*, 
is  constructed  from  the  union  of  all  the  elemental  axial  force 
vectors,  . 

D.  THE  ELEMEMTAL  STIFFMESS  MATRIX 

The  global  Galerkin  FEM  Equations  (4.22  and  4.35)  are 
constructed  from  the  union  of  elemental  axial  and  bending 
stiffness  matrices,  and  and  axial  and  lateral  force 

vectors,  and  .  For  the  beam  element,  the  elemental 

degrees  of  freedom  in  which  the  elemental  forces  act  are  shown 
in  Figure  (4.2)  . 


Tigur*  4.2  B«aa  ElMiMt  -  Dttgr««s  of  Troodom 
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Thus,  the  stiffness  matrix,  for  bending  results  in  a  4  X  4 
matrix  of  the  form; 


•11 

^12 

^13 

Kn 

bl 

21 

k,V 

>^23 

].bi 

K24 

.ti 

31 

^33 

j.bl 

X34 

bl 

41 

kti 

j-bi 

X43 

X44 

(4.36) 


For  the  bar  element,  the  elemental  degrees  of  freedom  in 
which  the  elemental  forces  act  are  shown  in  Figure  (4.3). 


Figur«  4.3  Bar  Blamant  Dagraas  of  Fraadom 

Thus  tn?  stiffness  matrix,  k"\  for  axial  force  results  in  a 
2X2  matrix  of  the  form: 
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Jr 

■^11 

0 

0 

■^12 

0 

0 

0 

■^11 

■^12 

0 

^13 

j,bJ 

■<14 

0 

K21 

K22 

0 

■^23 

■<24 

^21 

0 

0 

X22 

0 

0 

0 

•^31 

^32 

0 

k^i 

■<33 

1'*-' 

X34 

0 

Kii 

k^^ 

0 

X43 

■<44 

(4.38) 


The  elemental  displacements  and  forces  follow  suit  and  are 
defined  as  follows: 

The  elemental  displacements  vector,  6^',  becomes: 


(8^')-=<8r,5r,5f,8r,5r,8f  > 


(4.39) 


where  for  the  ith  element 


8/' 

8,^' 


the  axial  displacement  at  local  node  1 
the  lateral  displacement  at  local  node  1 
the  beam  slope  at  local  node  1 
the  axial  displacement  at  local  node  2 
the  lateral  displacement  at  local  node  2 
the  beam  slope  at  local  node  2 


The  elemental  force  vector,  f-*',  becomes: 


(f^')  ^  =  <  r/'r  fi',  fi',  fi',  ft'  > 


where  for  the  ith  element 

fi^’  =  the  axial  force  at  local  node  1 

f2^\  =  the  lateral  force  at  local  node  1 

=  the  moment  at  local  node  1 

f/'  =  the  axial  force  at  local  node  2 

^5^'  =  the  lateral  force  at  local  node  2 

fi^'  =  the  moment  at  local  node  2 


(4.40) 


Thus,  the  combination  of  the  Galerkin  Beam  and  Bar  Equations 


for  each  element  simplifies  to: 
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(4.41) 


where  the  elemental  stiffness  matrix,  k^‘ ,  in  terms  of  known 
quantities  becomes: 

AE/ti  0  0  -AE/ti  0  0 

0  12EI/ti^  6EI/ii  0  -12FJ/IJ  6EJ/ii 

0  €EI/ti^  ^El/ti  0  -6El/(i  2EJ/ti 

-AE/ti  0  0  AE/I^  0  0 

0  -12EI/^  -6EI/li  0  12EI/lii  -6£J/^ 

0  €EI/ii  2EI/ti  0  -6EI/ti 

It  should  be  noted  that  the  bar  and  beam  have  uncoupled 
behavior. 

E.  COORDINATE  TRANSFORMATION  OF  TBE  ELEMENTAL  SYSTEM 
OF  EQUATIONS 

For  curved  structures  such  as  the  arch,  each  element  has 
a  unique  orientation  with  respect  to  the  global  x  and  y  axes. 
Therefore,  to  solve  the  global  system  of  equations,  the 
elemental  Galerkin  Equation  (4.41)  is  transformed  from  local 
to  global  coordinates.  The  horizontal  and  vertical  axes  of  the 
arch  are  chosen  for  a  global  reference  coordinate  system. 
Figure  (4.5)  depicts  the  angle  the  i’'*’  element  makes  with  the 
horizontal  x-axis  as  a^,  and  the  compliment  angle.  Pi,  as  the 
angle  the  i^”  element  makes  with  the  vertical  y-axis. 

From  these  definitions,  the  local  displacements  and 
forces,  marked  by  a  prime  to  indicate  element  degree  of 
freedom  are  defined  in  terms  of  the  reference  coordinates  axes 
as  follows: 


(4.42) 
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Figure  4.5  Bar-Beam  Element  Degrees 
Transformation 


of 


Freedom 


=  5-"  cos  (a^)  +  8/  cos 

=  -8/  cos  +  8^^  cos  Caj 

=  5r 

=  8/  cos  (a^)  +  85^  cos  cPj 

=  -84''  cos  ^Pj  +  85^  cos 

=  5.^ 


=  f;-  COS  raj  +  fj  cos  rPj 

=  -f/  cos  rpj  +  cos  raj 

= 

=  f/  cos  raj  +  rj-*  cos  rpi-) 

=  -f/  cos  rpj  +  fs^  cos  raj 
=  f/ 


(4.45) 


(4.46) 


Accordingly,  a  transformation  matrix,  FJ  for  the  ith  element 


becomes ; 


cos(aj  cos(Pj  00  00 

-cos(pJ  cos(aj  00  00 

0  0  1  0  0  0 

0  0  0  cos  (aj  cos(Pj  0 

0  0  0  -cos(Pj  cos  (aj  0 

0  0  0  0  0  1 


(4.47) 


which  reduces  the  notation  of  Equations  (4.45)  and  (4.46)  to: 


8^'  =n8^ 


(4.48) 
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(4.49) 


1''  ' 

where 

(5^)  "  =  <5:,  bU  5^,  5i,  bi,  bi>  (4.50) 

(f ^)  ^  =  <  fi,  fi,  f-l,  fL  fi,  ft  >  (4.51) 

Thus,  the  transformed  elemental  stiffness  equation  becomes; 

(4.51) 

by  substituting  Equations  (4.48)  and  (4.49)  into  Equation 
(4.41)  . 

By  multiplying  both  sides  of  Equation  (4.52)  with  the 
inverse  of  the  transformation  matrix,  F* ,  an  orthogonal  matrix 
(i.e.,  r'‘*r^)  ,  yields: 

(P)7^i'(ri)5i«£i  (4.53) 

where  the  elemental  stiffness  matrix,  ^* ,  in  terms  of  the 
global  X  and  y  coordinates  is  defined  by: 

=  (r')^^*'(n)  (4.54) 


F .  SOLUTION 

Recall  from  the  Beam  and  Bar  FEM  development  that  the 
global  system  of  equations  result  from  the  union  of  the 
elemental  stiffness  matrices  and  force  vectors  such  that; 
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(4.55) 


ka  =  f 

where  the  global  or  system  force  vector,  F,  is  the  union  of 
the  transformed  local  force  vectors,  and  the  global  or 
system  stiffness  matrix,  K,  is  the  union  of  transformed  local 
stiffness  matrices,  .  Thus,  Equation  (4.55)  is  solved  for 
the  global  displacement  vector,  A. 

These  global  horizontal,  vertical,  and  rotational  degrees 
of  freedom  are  transformed  back  to  local  axial,  lateral,  and 
rotational  displacements  by  the  same  transformation 
: elat ionships  of  section  E  (Equations  4.45  and  4.46).  From 
these  local  displacements,  the  virtual  loads  at  the  element 
endpoints  are  computed  from  Equation  (4.41): 

=  (4.56) 

where  the  elemental  stiffness  matrix,  k'' ,  is  defined  by 
Equation  (4.42)  . 

The  node  point  virtual  loads,  f ,  equate  to  the  virtual 
axial  and  lateral  forces,  and  bending  moments  located  at  the 
endpoints  of  each  element.  From  Equation  (4.2)  and  (4.5), 
bending  and  axial  stresses  are  calculated.  For  continuity,  the 
stresses  of  internal  global  nodal  points  are  averaged  since 
physically,  local  nodal  point  2  of  the  element  is  the  same 
point  as  local  nodal  point  1  of  the  i’^*’  +  1  element.  Therefore, 
using  Equation  M.l),  the  normal  stresses  can  be  determined 
for  each  global  nodal  point. 
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V.  PROGRAM  DESCRIPTION  AMD  VALIDATION 


From  the  development  of  Chapters  II  and  IV,  a  VAX  Fortran 
77  Code  for  FEM  analysis  of  an  arch  was  written  to  interface 
with  the  DOT  software  package.  The  main  program,  ARCH_OPT . FOR, 
and  associated  common  program,  ARCH_COM.FOR,  are  contained  in 
Appendix  C.  Briefly,  ARCH_OPT.FOR  opens  and  reads  an  input 
file,  ARCH_IN.DAT,  before  it  is  divided  into  several 
subroutines  that  perform  the  FEM  analysis.  Table  5.1  lists  the 
input  data  fields  required  of  ARCH_IN.DAT  along  with  a  brief 
description  of  each. 


TABLE  5.1  ARCH  IN. DAT  FIELD  PARAMETERS 


Input  File  ANGLE,  RADIUS,  YOUNG,  YIELD, 

form:  NEL,  METHOD,  IPRINT, 

D"1BG,  DVILO,  DVIUP, 

DV2BG,  DV2LO,  DV2UP, 

CLAN,  FX,  FY,  FM,  FA, 

OPTDCS,  ITERATE,  PRCSN, 

BXl,  BYl,  BMl,  BX2,  BY2,  BM2, 
LABEL 

Parameter  Description 


Parameter 

ANGLE 


RADIUS 


YOUNG 


YIELD 


The  angle  from  0  to  359  degrees  subtended  by  the  arch 
structure . 

The  length  of  the  arch  radius  of  curvature.  (The  dimension 
is  arbitrary,  however  all  remaining  inputs  must  be 
consistent . ) 

Young'  Modulus  of  Elasticity  for  the  arch  material. 

The  yield  strength  of  the  arch  material.  If  a  factor  of 
safety  is  desired,  it  should  be  accounted  for  prior  to 
input . 

An  integer  number  of  elements,  from  1  to  32,  used  to 
approximate  the  arch  structure. 


METHOD  The  optimizer  method  to  be  used. 

METHOD  =  0  or  1 :  Modified  Method  of  Feasible  Directions 

METHOD  =  2:  Sequential  Linear  Programming 

ZPRZMT  On  screen  print  control  parameter.  Integers  from  0  to  5 
indicate  increasing  screen  printout. 

DV_BG  The  best  guess  for  design  variable  1,  the  base  dimension, 
or  2,  the  height  dimension.  Nodal  point  dimensions  are 
initialized  to  the  best  guess  value,  thus  establishes  the 
optimization  starting  point. 

Dy_LO  The  lower  limit  or  side  constraint  for  design  variable  1, 
the  base  dimension,  or  2,  the  height  dimension. 

DV_tJP  The  upper  limit  or  side  constraint  for  design  variable  1, 
the  base  dimension,  or  2,  the  height  dimension. 

CLAN  An  integer  from  1  to  NEL  +  1  that  indicates  the  node  at 

which  the  concentrated  load  is  to  be  applied. 

FX  The  magnitude  of  the  concentrated  load  in  the  horizontal 

direction  applied  at  node  CLAN. 

FY  The  magnitude  of  the  concentrated  load  in  the  vertical 

direction  applied  at  node  CLAN. 

FM  The  magnitude  of  the  concentrated  moment  applied  at  node 

CLAN. 

FA  The  magnitude  of  the  uniformly  distributed  load  in  the 

radial  direction  which  spans  the  entire  length  of  the  arch. 

OPTDCS  Optimization  option 

OPTDCS  =  1;  Optimize  the  dimensions  of  the  problem. 

OPTDCS  *  2:  Do  not  optimize  the  problem.  Based  on  the  initial 

design,  calculate  the  stress  distribution  only. 

ITERATE  The  number  of  iterations  performed.  The  resulting  optimized 
variables  are  re-entered  into  DOT  and  the  optimization 
performed  ITERATE  times  to  effect  an  iteration. 

PRCSN  Computer  precision  used  by  the  equation  solver. 


PRCSN  s  1: 

single 

precision 

PRCSN  =  2: 

double 

precision 
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BX_ 

Boundary  conditions  for  horizontal  displacement 

at  1 

! 

the 

first 

node  of  the  arch,  node  1,  or  2,  the  last  node 

of 

the 

arch. 

node  NEL  +  1 . 

BX_ 

s 

0: 

The  node  is  free  to  move  horizontally. 

BX_ 

= 

1: 

The  node  is  not  free  to  move  horizontally. 

By_ 

Boundary  conditions  for  vertical  displacement  at 

1, 

the 

first 

node  of  the  arch,  node  1,  or  2,  the  last  node 

of 

the 

arch. 

node  NEL  +  1 . 

BY_ 

0; 

The  node  is  free  to  move  vertically. 

BY_ 

= 

1; 

The  node  is  not  free  to  move  vertically. 

BM_ 

Boundary  conditions  for  the  beam  slope  at  1,  the 

first 

node 

of  the  arch,  node  1,  or  2,  the  last  node  of  the 

arch 

t 

node 

NEL  + 

1 . 

BM_ 

0: 

The  node  is  free  to  rotate. 

BM_ 

= 

1: 

The  node  is  not  free  to  rotate. 

LABEL 

A  character  string  used  to  identify  the  output. 

As  outlined  in  Figure  5.1,  the  main  program,  ARCH_OPT.FOR 
is  divided  into  subroutines.  In  general,  subroutine  Geometry 
is  called  in  order  to  generate  the  x  and  y  coordinates  of  the 
global  nodal  points  and  determine  the  orientation  of  each 
element.  Following  Geometry,  subroutine  Opt imization_tool 
establishes  the  DOT  parameters  prior  to  the  first  call  of  the 
DOT  program.  The  first  call  serves  only  to  record  the  DOT 
parameters  selected  in  DOT'S  internal  arrays.  After  DOT  is 
called,  the  Optimizat ion_tool  subroutine,  calls  Eval  to 
evaluate  the  objective  function  and  constraint  functions 
originally  outlined  in  the  problem  formulation  of  Chapter  II. 

As  detailed  in  Chapter  IV,  the  constraint  functions  are 
made  functions  of  the  design  variables  through  Finite  Element 
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Figure  5.1  Arch  Opt  Program  Structure 
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Method  analysis.  Subroutines  Form  and  Force_vector  develop  the 
global  stiffness  matrix  and  force  vector,  which  are  modified 
by  subroutine  Bndary  for  the  appropriate  boundary  conditions. 
The  equation  solver,  L2ARG,  from  the  IMSL  library  is  called  to 
solve  for  the  global  displacements,  which  in  turn  are  used  to 
calculate  the  nodal  stresses.  Once  the  constraints  are 
evaluated  for  the  initial  design,  the  problem  is  returned  to 
DOT  where  the  move  direction  is  computed  and  an  updated  design 
point  chosen.  The  objective  and  constraint  functions  are 
reevaluated  for  the  updated  design  point  before  returning  to 
DOT  for  further  iteration. 

Once  termination  criteria  for  optimization  are  reached, 
the  main  program  creates  the  output  file,  ARCH_OUT.DAT.  This 
file  contains  the  problem  parameters,  optimized  design 
variables,  and  the  resulting  objective  function  value  along 
with  a  variety  of  additional  information.  Summarizing,  for  a 
given  geometry,  loading,  and  set  of  end  conditions,  the 
program  is  capable  of  finding  the  optimum  cross  section 
dimensions  of  each  nodal  point  along  the  length  of  the  arch. 

To  validate  the  FEM  analysis,  several  non-optimum  straight 
beam  and  arch  problems  with  known  analytical  solutions  were 
solved.  A  straight  cantilever  structure,  subject  to  a 
concentrated  lateral  end  load,  axial  load,  and  end  moment;  and 
a  quarter  cantilever  arch,  subject  to  a  lateral  end  load  were 
analyzed.  These  test  problems  established  the  program  error 
for  stress  and  displacement  calculations.  Additionally,  the 
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quarter  cantilever  arch  and  a  hinged-hinged  semi-circular  arch 
structure,  subject  to  a  lateral  load  on  the  axis  of  symmetry, 
establish  trends  in  a  relationship  between  the  number  of 
elements  used  to  approximate  the  arch  and  accuracy.  The 
remainder  of  this  chapter  is  a  summary  of  the  results  and  the 
conclusions  drawn  from  each  validation  problem  studied.  The 
complete  solution  of  each  problem  is  contained  in  Appendix  D. 


A.  VALIDATION  I:  CANTILEVER  BEAM 

A  cantilever  beam  is  subject  to  a  concentrated  end  load  as 
shown  in  Figure  5.2.  [Ref.  8] 


s 

4 

3 

Z 


0.0 


HM 


a2.s 


3Z7S 


1  \  4S. 

mfnrr 

Figure  5.2  Validation  Case  #1 


£,  =  45.00  inches 
S  =  1 . 50  inches 
H-3.00  inches 

I  =  45^  =3.375  inches* 

P  =  1,000.0  ibf 

M  =  y 


) 


ARCH_OPT.FOR  was  run  for  this  beam  structure  using  an 
angle  of  45.0  X  10"’  radians  and  a  radius  of  10‘  inches  to 
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approximate  a  straight  beam  of  45  inches.  The  four  element  FEM 
solution  is  compared  to  the  analytical  solution  in  Table  5.2. 

TABLE  5.2 


NODE 

THEORETICAL 

STRESS 

FEM  ANALYTICAL 
STRESS 

%  EBROR 

1 

20000.0 

19999.7 

0.0015 

2 

15000.0 

14999.7 

0.0020 

3 

10000.0 

9999.8 

0.0020 

4 

5000.0 

4999.9 

0.0020 

5 

o 

• 

o 

0 

0.0000 

where  the  percent  error  is  defined  as: 

%  Error  =  (  Theory  -  FEM  Analysis  )/  Theory  *  100 

B.  VALIDATION  II:  PRISMATIC  BAR 

Similarly,  a  prismatic  bar  is  subject  to  an  axial  load  as 
shown  in  Figure  5.3.  [Ref.  8] 

L  =  45 . 00  inches 
B  =  1 .50  inches 
B  =  3 . 00  inches 
A  =  3  •  B  =  4 . 50  inches^ 

i 

1,  000.0  Ibf 
a  =  222.7 

Figure  5.3  Validation  Casa  #2 


Input  values  for  angle  and  radius  remained  the  same  to 
approximate  the  straight  bar.  The  four  element  FEM  solution  is 
compared  in  Table  5.3. 


TABLE  5.3 


1  NODE 

THEORETICAL 

STRESS 

FEM  ANALYTICAL 
STRESS 

%  ERROR 

bbi 

222.2 

222.2 

0.0000 

222.2 

222.2 

0.0000 

222.2 

222.2 

0.0000 

1^^ 

222.2 

222.2 

0.0000 

222.2 

222.2 

0.0000 

C.  VALIDATION  III:  CANTILEVER  BEAM 

The  cantilever  beam  is  subject  to  a  concentrated  moment  at 
the  free  end  as  shown  in  Figure  5.4.  [Ref.  8] 


L  =  45 . 00  inches 

B  =  1 . 50  inches 

H  =  3 . 00  inches 

37.76" 

1  =  4^  =  3.375  inches* 

21.S" 

E  =  30  xlO^psi 

M=10,  OOOibf 

11.25" 

Slope  S' 

0. 

.  Mx^ 

Displacement  S  - 

Figur*  5.4  Validation  Casa  #3 
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The  four  element  FEM  solution  for  both  slope  and 


displacement  is  compared  in  Table  5.4. 


TABLE  5.4 


NODE 

THEORETICAL 

SLOPE 

FEM  ANALYTICAL 
SLOPE 

%  ERROR 

1 

0.00000000 

0.00000000 

0.0000 

2 

0.00111111 

0.00111109 

0.0019 

3 

0.00222222 

0.00222218 

0.0019 

4 

0.00333333 

0.00333328 

0.0016 

5 

0.00444444 

0.00444438 

0.0014 

NODE 

THEORETICAL 

DISPLACEMENT 

FEM  ANALYTICAL 
DISPLACEMENT 

%  ERROR 

1 

0.00000000 

0.00000000 

0.0000 

2 

0.00625000 

0.00624985 

0.0024 

3 

0.02500000 

0.02499940 

0.0024 

4 

0.05625000 

0.05624880 

0.0021 

5 

0.10000000 

0.09999790 

0.0021 

D.  DATION  IV:  CANTILEVER  QUARTER  ARCH 

A  cantilever  quarter  arch  is  subject  to  a  lateral  load  as 
shown  in  Figure  5.5.  [Ref.  8] 


Figure  5.5  Validation  Casa  #4 


L  =  45 . 00  inches 
B  =  1 . 50  inches 
H  =  3 . 00  inches 

1  =  4^  =  3.375  inches^ 

E  =  30  X 10*  psi 
P  =  l,  000.01bf 

s  _ 


ARCH_0PT.F0R  was  run  for  this  structure  using  an  angle  of 
90.0  degrees  and  a  radius  of  45  inches.  To  approximate  the 
arch  a  four,  six,  eight,  ten  and  12  element  FEM  solution  is 
solved  and  compared  to  the  analytical  solution  presented  in 
Table  5.5. 


TABLE  5.5 


NODE 

%  ERROR 

4 

0.450000 

0.446951 

0.677556 

1  ^ 

0.450000 

0.448382 

0.359556  1 

1  ^ 

0.450000 

0.448854 

0.254667 

1 

0.450000 

0.448790 

0.268889 

I  12 

0.450000 

0.449100 

0.200000 
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E.  VALIDATION  V:  HINGED-HINGED  SEHI -CIRCULAR  ARCH 


A  hinged-hinged  semi-circular  arch  structure  is  subject  to 
a  lateral  load  along  the  axis  of  symmetry  as  shown  in  Figure 
5.6.  [Ref.  8 ] 


JR  =  32 . 0  inches 
S  =  1 . 50  inches 
H  =  3.00  inches 


P=  10,  000.0  Ibf 

(l-cosG)  {BmB) 


Results  are  tabulated  in  Table  5.6  for  comparison  of  the 
four,  six,  eight,  ten,  12,  14,  and  16  element  FEM  solutions  to 
the  analytical  solution.  It  should  be  noted  that  by  using 
symmetry,  the  arch  structure  is  approximated  by  twice  the 
number  of  elements  shown  in  the  calculations. 


51 


TABLE  5.6 


NODE 

e 

THEORETICAL 

STRESS 

FEM  ANALYTICAL 
STRESS 

%  ERROR 

4 

0 

0.0 

0 

0.000000 

22.5 

11911.4 

12338 

3.581870 

45 

11183.3 

11971.7 

7.049945 

67.5 

2073.4 

1043.2 

49.685474 

90 

25840.4 

24725.3 

4.315231 

6 

0 

0.0 

0 

0.000000 

15 

9293.9 

9411.2 

1.262367 

30 

13108.3 

13334.9 

1.728775 

45 

11183.3 

11503.9 

2.866917 

60 

3650.1 

4042.8 

10.759982 

75 

8978.0 

8539.9 

4.879970 

90 

25840.4 

25386.8 

1.755283 

8 

0 

0.0 

0 

0.000000 

11.25 

7465.5 

7509.3 

0.586658 

22.5 

11911.4 

11997.3 

0.721573 

33.75 

13166.7 

13291.6 

0.948666 

45 

11183.3 

11342.3 

1.421903 

56.25 

6037.3 

6224.5 

3.099877 

67.5 

2073.4 

1865.4 

10.029988 

78.75 

12837.1 

12616.3 

1.720372 

90 

25840.4 

25615.1 

0.87178 

10 

0 

0.0 

0 

0.000000 

9 

6206.4 

6225.1 

0.301171 

18 

10509.0 

10546.3 

0.352057 

27 

12801.8 

12856.2 

0.424639 

36 

13028.5 

13098.8 

0.539901 

45 

11183.3 

11268.0 

0.757518 

54 

7311.7 

7408.7 

1.325965 

63 

1509.2 

1616.0 

7.077882 

72 

6081.5 

5967.4 

1.876789 

81 

15273.5 

15155.0 

0.775860 

90 

25840.4 

25720.4 

0.464280 
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NODE 


e 

THEORETICAL 

STRESS 

FEM  ANALYTICAL 
STRESS 

%  ERROR 

0 

.0 

0 

0.000000 

15 

9293.9 

9309.8 

0.171326 

30 

13108.3 

13139.3 

0.236589 

45 

11183.3 

11227.5 

0.395371 

60 

3650.1 

3704.7 

1.497107 

75 

8978.0 

8917.1 

0.678600 

90 

25840.4 

25777.3 

0.244082 

0 

0.0 

0 

0.000000 

6.428571 

4621.6 

4624.5 

0.063120 

12.85714 

8290.8 

8296.5 

0.06S954 

19.28571 

10961.5 

10970.1 

0.078835 

25.71429 

12600.0 

12611.5 

0.091081 

32.14286 

13185.9 

13200.1 

0.107896 

38.557143 

12711.6 

12728.6 

0.133432 

45 

11183.3 

11202.7 

0.173611 

51.42857 

8620.0 

8641.9 

0.253703 

57.85714 

5054.1 

5078.2 

0.476611 

64.28571 

530.4 

556.4 

4.907803 

70.71429 

4894.3 

4866.7 

0.563998 

77.14286 

11151.7 

111123.1 

0.256399 

83.57143 

18163.1 

18133.9 

0.160800 

90 

25840.4 

25810.9 

0.114053 

0 

0.0 

0 

0.000000 

11.25 

7465.5 

7465 

0.006738 

22.5 

11911.4 

11911.1 

0.002106 

33.75 

13166.7 

13167.7 

0.007655 

45 

11183.3 

11186 

0.024281 

56.25 

6037.3 

6042.2 

0.080341  1 

67.5 

2073.4 

2066.5 

0.330744  1 

78.78 

12837.1 

12828.8 

0.065020  1 

90 

25840.4 

25831.4 

0.034720  1 
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F .  CONCLUSIONS 


The  four  element  approximation  for  a  straight  cantilever 
structure  produced  an  error  no  greater  than  0.016%.  The 
cantilever  quarter  arch  produced  an  error  less  than  0.70%  for 
the  four  element  model  which  reduced  to  less  than  0.20%  with 
12  elements.  The  results  of  the  hinged-hinged  arch  indicate, 
as  expected,  that  the  more  elements  used  the  better  the 
solution.  Considering  only  meaningful  stresses,  stresses  in 
excess  of  10,000  psi,  the  error  is  less  than  2%  for  eight 
elements  and  less  than  0.8%  for  12  elements. 

In  general,  the  percent  error  recorded  for  the  first  three 
validation  cases  seemed  insignificant.  Four  element 
approximations  sufficed  to  solve  the  stresses,  slopes,  and 
displacements  for  straight  structures.  Therefore,  it  was 
concluded  that  the  program  was  producing  accurate  results  for 
analysis  of  straight  beams. 

Unfortunately,  for  the  arch  structures,  the  error  of  the 
four  element  model  was  significant  (greater  than  45%)  . 
However,  the  error  reduced  significantly  when  more  elements 
were  used  to  approximate  the  structure.  Grid  independence,  (2% 
error) ,  was  not  achieved  for  the  hinged-hinged  arch  until  at 
least  eight  elements  are  used  to  approximate  the  structure. 
This  indicates  that  an  element  cannot  be  used  to  span  more 
than  11.25  degrees  of  arch.  The  resulting  trend,  as  expected, 
confirms  that  the  more  elements  used,  the  better  the  model. 
However,  computer  time  and  computer  error  increase  with 
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increase  in  the  number  of  elements  and  models  of  more  than 
eight  elements  were  not  used. 
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VI.  CASE  STUDIES 


Results  are  presented  for  a  number  of  cises  with  regard  to 
optimization  scheme  and  stress  analysis.  The  case  studies 
range  from  the  simple  cantilever  beam  to  complex  arch  struc¬ 
tures.  In  addition,  for  many  cases  one  parameter  of  the  same 
structure  was  modified  and  the  problem  was  reoptimized  to 
establish  a  comparison.  The  straight  beam  is  examined  first, 
followed  by  five  cases  studying  the  quarter  cantilever  arch 
with  varied  loadings.  Cases  #7  and  #8  are  symmetric 
semicircular  arches  comparing  simply  supported  arch  structures 
with  fixed-end  arch  structures.  The  remaining  cases  are 
asymmetric  semicircular  arch  structures.  Cases  #9  through  #11 
investigate  various  end  conditions  and  Cases  #12  through  #14 
various  combined  loadings.  The  cases  conclude  with  Case  #15 
which  combines  a  concentrated  lateral  load,  applied  moment, 
and  distributed  load  across  the  arch  structure. 

For  each  case,  interpretations  of  the  results  are 
accompanied  by  a  schematic  drawing  of  the  structure  modeled, 
a  plot  of  the  cross  section  dimensions  and  area  as  functions 
of  nodal  points,  and  a  plot  of  the  axial  and  bending  stresses 
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as  functions  of  nodal  points.  Eight  elements  were  selected  to 
model  the  arch  structures  and  the  material  properties  were 
selected  such  that  the  yield  strength  was  imputed  as  52,000 
psi  and  Young's  Modulus  as  30,000,000  psi.  For  reference,  the 
Modified  Method  of  Feasible  Directions  will  be  referred  to  as 
Method  1  and  the  Sequential  Linear  Programming  Method  will  be 
referred  to  as  Method  2.  Additionally,  each  endpoint,  unless 
geometrically  restricted  by  imposed  boundary  conditions,  can 
have  three  'means  of  displacement,'  MOD.  An  endpoint  can 
rotate  about  the  z-a.xi3,  displace  in  the  x  direction,  and 
displace  in  the  y  direction.  For  reference,  an  endpoint  will 
be  described  by  a  number  from  zero  to  three  reflecting  the 
means  of  diplacen.ent .  As  an  example,  a  fixed  end  is  considered 
to  have  zero  means  of  displacement  because  it  cannot  rotate  or 
displace  in  either  the  x  or  y  direction.  A  free  end  which  can 
rotate  and  displace  in  both  the  x  and  y  direction  is 
considered  to  have  three  means  of  displacement.  A  hinge  which 
can  only  rotate  has  one  MOD.  The  complete  computer  data 
printout  is  presented  in  Appendix  D. 
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A.  CASE  #1:  CANTILEVER  BEAM  WITH  LATERAL  LOADING 

The  cantilever  beam  was  optimized  first  in  order  to 
provide  guidance  for  adjusting  the  various  parameters 
discussed  in  Chapter  III.  Satisfactory  results  were  produced 
by  turning  the  auto  scaling  function  off,  reducing  CT  and 
CTMIN,  and  establishing  the  termination  criteria.  Using  the 
Modified  Method  of  Feasible  Directions,  henceforth  referred  to 
as  Method  1,  the  cross  section  dimensions  and  stresses  were 
plotted.  As  expected,  the  dimensions  form  a  parabolic  function 
over  the  length  of  the  beam.  Furthermore,  the  beam  exhibits 
only  stress  due  to  bending  moments.  The  normal  stresses  are 
virtually  nonexistent  which  likewise  is  as  expected. 


Loads 
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B.  CASE  #1A:  CANTILEVER  BEAM  WITH  LATERAL  LOADING 

For  comparison  with  Case  #1,  the  same  cantilever  beam  was 
optimized  using  the  Sequential  Linear  Programming  Method, 
henceforth  referred  to  as  Method  2.  The  results  are  quite 
similar.  In  total  structure  volume,  the  difference  is  less 
than  0.07%.  The  only  significant  difference  appears  at  nodal 
point  9,  the  free  end.  In  theory,  the  free  end  of  a  beam  can 
support  no  bending  stresses.  For  this  case,  nodal  point  9  has 
no  stresses  unlike  the  previous  case  which  had  relatively 
small  bending  stresses  at  nodal  point  9.  However,  from  this 
result  alone  it  is  not  conclusive  that  Method  2  is  superior  to 
Method  1 . 


Loads 


Lateral 

=  2,000  lbs 

Axi  al 

=  C  lbs 

Moment 

=  0  in-lbs 

End  conditions 

Node  L 

0  MOO 

Node  9 

3  MOO 

Dimensions 

Radius 

=  not  applicable 

Theta 

=  not  applicable 

Total  volume 

Volume 

33.15  in’ 
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C.  CASE  #2:  CANTILEVER  ARCH  WITH  LATERAL  LOADING 

Since  Case  #1  did  not  strongly  suggest  a  preferential 
method,  the  cantilever  arch  was  optimized  with  Method  1.  At 
most  nodal  points,  the  total  stresses  were  well  below  the 
yield  stress  which  indicates  that  this  design  is  far  from  an 
optimum  structure.  Additionally,  the  height  and  base 
dimensions  hovered  around  the  initial  starting  point  of  2 
inches  by  2  inches  and  produced  a  structure  only  7.42%  less  in 
volume  than  that  of  the  initial  structure.  It  appears  that  the 
optimizer  failed  to  achieve  an  optimum  solution  using  this 
Method . 


Loads 


Lateral 

=  2,000  lbs 

Axial 

=  C  lbs 

Komer.t 

=  0  in-lbs 

End  conditions 

Node  1 

0  MOD 

Node  9 

1  KC2 

Dimensions 

Radius 

=  32  in 

Tr.eta 

=  90  degree: 

Total  volume 

VoiuTe 

=  186.15  in' 
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D.  CASE  #2A:  CANTILEVER  ARCH  WITH  LATERAL  LOADING 

For  comparison,  the  same  arch  structure  was  reoptimized 
with  Method  2.  Each  element  of  the  structure  now  supports 
stresses  equal  to  the  yield  stress  producing  an  efficient 
structure.  The  total  volume  was  reduced  from  the  initial 
starting  point  by  61.32%.  For  this  structure.  Method  2  also 
produced  results  with  fewer  iterations  than  Method  1.  With 
these  observations  in  mind.  Method  2  was  selected  as  the 
preferred  method  for  quarter  arches.  Additionally,  it  is 
interesting  to  note  that  the  axial  stresses  only  remotely 
effect  the  stress  total  for  the  first  5  nodal  points,  hence 
the  first  45  degrees  of  arch.  After  node  5,  the  height  reduces 
significantly,  however  the  area  remains  roughly  the  same. 


Loads 


Lateral 

2,0C0  lbs 

Axial 

=  0  lbs 

Moment 

=  0  in-lbs 

End  conditions 

Node  1 

0  MOD 

Node  9 

3  MOD 

Dimensions 

Radius 

=  32  in 

Theta 

=  90  degrees 

Total  volume 

Volume 

=  77.78  in’ 
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E.  CASE  #3:  CANTILEVER  ARCH  WITH  AXIAL  LOADING 

This  case  presents  a  quarter  arch  structure  subject  to  an 
axial  load  vice  the  lateral  load  of  Case  #2.  Unlike  Case  #2 
the  axial  stresses  increase  significantly  along  the  length  of 
the  arch  and  the  bending  stresses  decrease.  The  net  result  is 
an  arch  structure  of  27,15%  less  material.  This  seems  to 
indicate  a  dominant  relationship  between  area  and  bending 
stress.  Additionally,  this  case  exemplifies  the  difficulty 
experienced  by  approximating  an  arch  of  90  degrees  with  eight 
straight  segments.  The  plots  appear  very  disjointed,  hence  the 
data  points  seem  circumspect.  However,  the  effect  can  be 
minimized  as  presented  in  Case  #3a. 


ZOOO^ 


Loacs 


Lateral 

=  0  lbs 

Axial 

=  2, COO  lbs 

Mc.-nenc 

=  0  in-lbs 

End  conditions 

Node  1 

0  KCD 

Node  9 

3  MOO 

Dimensions 

Radius 

=  32  in 

Theta 

=  90  degrees 

Total  volume 

Volume 

=  56.66  in’ 
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F.  CASE  #3A:  CANTILEVER  ARCH  NITH  LATERAL  LOADING 

Thus  far,  each  case  has  started  with  an  initial  design  of 
2  inches  by  2  inches.  For  comparison,  the  arch  structure  of 
Case  #3  was  optimized  a  second  timeusing  the  results  of  Case 
#3  as  the  initial  design.  Reoptimizing  had  the  desired  effect 
of  smoothing  the  results  and  in  graphical  form,  both  the  area 
and  stress  curves  take  on  a  fairer  shape.  In  terms  of  total 
structure  volume,  the  reoptimized  arch  was  27.03%  smaller  than 
that  of  Case  #3.  In  all  subsequent  cases  this  strategy  of 
reoptimization  will  be  referred  to  as  a  two-stage  optimization 
strategy. 


Lateral 

Loads 

0  lbs 

Axial 

s 

2, CCD  lbs 

Women t 

0  in-lbs 

End  conditions 

Node  1 

0  MOD 

Node  9 

3  KOD 

Radius 

Dimensions 

32  in 

Theta 

9C  degree: 

Total  volume 


If 


Voi'wroe 


41.34  in’ 


Case  #3a 


Node 


^  Base 


Height 


Area  (in  ^2) 


6.  CASE  #4;  CANTILEVER  ARCH  NITH  LATERAL  LOADING  AND  MCX4ENT 

For  this  case,  a  lateral  load  and  concentrated  moment  were 
applied  at  nodal  point  9.  The  shape  of  the  dimension  plot 
curves  are  very  similar  to  those  of  the  cantilever  beam, 
parabolic  in  form.  In  comparison  with  the  same  structure 
subject  only  to  the  lateral  load.  Case  #2,  the  total  structure 
volume  is  reduced  by  18.67%.  The  concentrated  end  moment 
negates  the  effect  of  lateral  load  on  the  extreme  fibers  by 
producing  compressive  stresses  on  the  outer  fibers  and  tensile 
stresses  on  the  inner  fibers  of  the  arch.  Thus,  the  cross 
sectional  dimensions  necessary  to  withstand  the  total  normal 
stress  is  reduced  thereby  reducing  the  total  structure  volume. 


Lateral 

Axial 

Moment 


Nooe  1 
Noae  9 


Radius 

Theta 


Volume 


2, CCS  lbs 
C  lbs 

l,CuO  in-lbs 


C  MCL 
3  MOL 

Dimensions 


32  in 

9C  degrees 


“Ota.  volum.e 


63.25  in’ 
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H.  CASE  #4A:  CANTILEVER  ARCH  WITH  LATERAL  LOADING  AND  MC»<ENT 


To  further  emphasize  the  effect  of  the  concentrated  end 
moment,  the  structure  of  Case  #4  was  subject  to  the  same 
lateral  load  while  the  moment  at  the  end  point  was  increased 
by  a  factor  of  10.  By  increasing  the  applied  moment,  the 
effect  of  the  lateral  load  on  the  extreme  fibers  is  negated 
further  which  reduces  the  cross  sectional  area  necessary  to 
withstand  the  total  stresses.  Expectedly,  the  volume  reduced 
from  Case  #4  by  15.88%  for  a  total  reduction  from  Case  #2  of 
31.60%.  It  is  interesting  to  note  that  the  shape  of  the 
dimension  curves  still  remain  parabolic  in  form. 


Loads 
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I.  CASE  #5:  CANTILEVER  ARCH  WITH  DISTRIBUTED  LOADING 

This  case  is  presented  to  display  some  of  the  versatility 
of  the  program.  A  load  acting  radially  inward  is  distributed 
along  the  length  of  the  arch.  The  cross  section  dimensions  and 
area  curves  appear  to  be  almost  linear  and  the  bending 
stresses  dominate  the  total  stresses.  Since  the  bending  stress 
is  a  function  of  height  squared,  the  optimizer  tried  to 
maximize  the  height  dimension  until  the  geometric  constraint 
was  violated.  At  each  nodal  point,  the  height  is  10  times  the 
size  of  the  base  except  at  the  end  point  for  which  both 
dimensions  reach  the  minimum  side  constraint.  Had  the  arch 
structure  not  been  optimized,  the  volume  necessary  to  support 
the  distributed  load  would  increase  by  225%. 


Loads 


Lateral 

Axial 

Moment 

Distrib. 

=  0  lbs 

=  0  lbs 

=  0  in-lbs 

=  1,000  Ibs/in 

End  conditions 

Node  1 
Node  9 

0  MOO 

3  MOD 

Dimensions 

Radius 

Theta 

=  32  in 

=  90  degrees 

Total  volume 

Volume 

=  55.70  in’ 
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Case  #5 


J.  CASE  #6:  CANTILEVER  ARCH  WITH  LATERAL  AND 
DISTRIBUTED  LOADING 

To  build  on  Case  #5,  a  lateral  load  was  applied  at  the  end 
point  in  addition  to  the  distributed  load.  In  comparison,  the 
volume  required  to  withstand  the  lateral  load  only  is  77.78 


in^  (Case 

#2)  . 

The 

volume  required 

to 

withstand 

the 

distributed 

load 

only 

is  55.70  in^. 

Yet 

the  volume 

to 

withstand  both  the  lateral  load  and  the  distributed  load 
presented  in  this  case  is  97.47  in^.  By  combining  loads  which 
produce  opposing  bending  moments,  the  volume  of  the  resultant 
optimized  arch  is  not  equal  to  the  sum  of  the  volume  of  arches 
optimized  subject  to  the  individual  loads.  Therefore,  it  is 
possible  to  achieve  a  more  efficient  structure  through 
resourceful  combination  loadings. 


Loads 


Lateral 

Axial 

Mo.'^ent 

Di str lb . 

=  2,000  lbs 

=  0  lbs 

0  in-lbs 

=  100  Ibs/in 

End  conditions 

Node  1 
Node  9 

3  HOD 

Dimensions 

Radius 

Theta 

=  32  in 

=  90  degrees 

Total  volume 

Volume 

=  97.  <5  7  in’ 
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K.  CASE  #6A:  CANTILEVER  ARCH  WITH  LATERAL  AND 

DISTRIBUTED  LOADING 

For  comparison,  the  same  structure  (Case  #6)  was  optimized 
with  the  DOT  auto  scaling  function  switched  on.  Changing  this 
parameter  seemed  to  have  little  effect  on  the  overall  volume 
indicated  by  an  increase  by  only  4.40%.  However,  the 
computation  effort  judged  by  total  computer  time  nearly 
doubled  and  both  the  dimension  and  stress  curves  have 
unexpected  behavior  near  the  endpoint.  This  comparison 
confirmed  that  better  results  were  achieved  by  switching  the 
auto  scaling  function  off  for  these  structures. 


Lateral 
Axial 
MoT.er.t 
Di strib. 


Node  1 
\'ode  9 


Radius 
"net  a 


Vclure 


Loads 

=  2,000  lbs 

=  0  lbs 

=  0  in-ibs 

=  ICO  Ibs/in 

End  conditions 

C  MOD 
3  yOD 


DiPiCnsi  ons 


=  32  in 

=  90  degrees 

Total  volu-e 

=  101. -76  in’ 
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L.  CASE  #7:  HINGED-HINGED  ARCH  WITH  LATERAL  LOADING 


For  the  remaining  cases,  it  was  observed  that  the  only 
reliable  and  consistent  results  were  obtained  by  using  Method 
1  for  optimization.  It  is  theorized  that  restricting 
displacements  at  both  endpoints  may  have  caused  Method  2  to 
become  mathematically  unstable  and  therefore  unsuitable  to 
solve  such  problem.  For  this  particular  case,  it  is 
interesting  to  note  that  at  the  base,  node  1,  and  56.25 
degrees  from  the  base,  node  6,  the  axial  stress  completely 
dominates  the  total  stresses  because  there  is  virtually  no 
binding  force.  At  these  points,  the  dimensions  of  the  cross 
section,  dictated  strictly  by  the  axial  stress,  form  a  square 
to  produce  the  minimum  area. 

Leads 

Lateral  =  16,000  lbs 

Axtal  =  0  lbs 
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M.  CASE  #7A:  HIMGE-BIMGED  ARCH  WITH  LATERAL  LOADING 

Similar  to  Case  #3a,  the  arch  structure  of  Case  #7  was 
reoptimized  using  the  results  achieved  as  the  initial  design 
in  order  to  apply  the  two-stage  optimization  strategy.  Again, 
reoptimizing  had  the  desired  effect  of  smoothing  the  results, 
however,  this  effect  was  not  as  dramatic  for  Method  1  as  for 
Method  2.  The  two-stage  optimization  strategy  only  reduced  the 
volume  by  4.47%  using  Method  1  as  opposed  to  the  27,03% 
reduction  using  Method  2.  Additionally,  at  node  6,  the  total 
stresses  exceeded  the  yield  stress  by  2.38%.  Fortunately,  this 
occurrence  did  not  repeat  in  any  other  cases  due  to 
reoptimization.  Therefore,  the  two-stage  optimization  strategy 
was  applied  for  the  remaining  cases. 


Loads 


Lateral 

=  16, 

WWW  1 Os 

Axial 

=  0  : 

DS 

Moment 

=  0  in-lbs 

End  conditions 

Node  1 

1  yoD 

Ncae  9 

1  MOD 

Dimensions 

Radius 

=  32 

in 

Theta 

18S 

degrees 

Total  volume 

Volume 

=  123 

.35  in* 
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N.  CASE  #8:  FIXED-FIXED  ARCH  WITH  LATERAL  LOADING 

For  this  case,  the  same  loading  of  Case  #7  was  applied  to 
a  semicircular  arch  with  fixed  end  points.  This  produces  a 
statically  indeterminate  structure  with  zero  irieans  of 
displacement  at  both  the  boundaries.  As  a  result,  the  peaks  of 
the  axial  stress  curve  are  dampened  and  shifted  towards  the 
center  by  approximately  15  degrees.  A  larger  bending  moment  is 
produceu  at  the  base  since  it  is  no  longer  free  to  rotate. 
However,  the  net  results  is  that  the  total  structure  volume  of 
Case  #7  is  reduced  by  14.08%  by  changing  the  end  conditions 
from  simply-supported  to  fixed.  From  this,  as  expected,  a 
structure  more  statically  indeterminate  results  in  a  more 
efficient  structure. 


Loads 


Lateral 

-■T 

16,CCC  lbs 

Axial 

0  lbs 

Kor.er.z 

= 

0  in-lbs 

Er.d  cor.diticr 

,£ 

Node  1 

C  MCL 

Noce  9 

c  y.cL 

Dir'.ensior.s 

Radi  us 

=. 

32  in 

Tneta 

= 

180  degrees 

Total  volume 

Voiure 

= 

1  05.98  in’ 
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O.  CASE  #8A:  FIXED-FIXED  ARCH  WITH  LATERAL  LOADING 

In  Chapter  I,  thesis  research  performed  by  Scott  McDavid 
was  mentioned  as  the  predecessor  for  this  investigation.  His 
results  optimized  a  fixed-fixed  arch  subject  to  the  same 
lateral  load  with  respect  to  the  base  dimension  only.  By 
holding  the  height  dimension  constant,  the  structure  must  have 
twice  the  volume  in  order  to  withstand  the  loading.  For  most 
arch  structures,  the  bending  stress  is  usually  the  more 
dominate  stress.  Therefore,  the  height  dimension  has  more 
effect  on  the  total  volume  than  the  base  dimension  because  the 
bending  stress  is  a  function  of  base  times  height  squared. 
When  the  bending  stress  dominates,  the  optimizer  will  seek  to 
maximize  the  height. 


Loads 


Lateral 

=  16,000  lbs 

Axial 

=  0  lbs 

Moment 

=  0  in-lbs 

End  conditions 

Node  1 

0  MOD 

Node  9 

0  MOD 

Dimensions 

Radius 

=  32  in 

Theta 

180  degrees 

Total  volume 

Volume 

=  215.17  in* 
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P.  CASE  #9:  FIXED-HINGED  ROLLER  ARCH  WITH  AXIAL  LOADING 

Unable  to  invoke  symmetry  on  the  remaining  cases,  the 
elements  used  to  model  an  asymmetric  semicircular  arch  must 
span  the  full  180  degrees.  The  largest  number  of  elements  used 
to  model  the  structure  and  produce  consistent  results  remained 
only  eight.  It  is  again  suggested  that  restricting 
displacements  at  both  endpoints  cause  the  problem  to  become 
mathematically  unstable.  For  this  particular  case,  the  arch 
has  zero  means  of  displacement  at  node  1  and  two  means  of 
displacement  at  node  9.  As  expected,  the  arch  is  quite  large 
at  node  1  to  support  the  resultant  moment.  At  node  3,  45 
degrees  up  from  node  1,  the  axial  stress  dominates  the  total 
stress  and  the  size  decreases.  For  reference,  this  arch  is 
more  than  twice  the  volume  of  the  arches  in  Case  #7  and  #8. 


Loads 


Lateral 

=  0  lbs 

Axial 

=  16,000  lbs 

Moment 

=  0  in-lbs 

End  conditions 

Node  1 

0  MOD 

Node  9 

2  MOD 

Dimensions 

Radius 

=  32  in 

Theta 

=  180  degrees 

Total  volume 

Volume 

287.15  in’ 
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Dimension  (in) 


Q.  CASE  #9A:  HINGED-HINGED  ROLLER  ARCH  WITH  AXIAL  LOADING 

To  emphasis  the  conclusion  drawn  from  case  #9  about  the 
endpoints,  the  arch  structure  and  loading  studied  for  Case  #9 
was  modified  by  adding  an  additional  means  of  displacement  at 
node  1.  Allowing  nodal  point  1  to  rotate  freely,  the  dimension 
and  stress  curves  alter  drastically.  The  total  structure 
volume  increased  by  19.92%,  yet  the  structure  cannot  withstand 
the  stresses.  The  total  stresses  exceed  the  yield  stresses  by 
54.81%  resulting  in  an  infeasible  design.  It  appears  that  the 
optimizer  failed  to  achieve  an  optimal  solution  for  this  arch 
structure  due  to  the  additional  means  of  displacement. 


Loads 


Lateral 

=  2,000  lbs 

Axial 

=  16,  000  lbs 

Moment 

=  0  in-lbs 

End  conditions 

Node  1 

1  MOO 

Node  9 

2  MOO 

Dimensions 

Radius 

=  32  in 

Theta 

=  180  degrees 

Total  volume 

Volume 

=  34<.34  in’ 
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R.  CASE  #10:  FIXED-FIXED  ROLLER  ARCH  WITH  AXIAL  LOADING 

For  this  case,  the  same  arch  structure  and  loading  studied 
in  Case  #9  was  modified  by  reducing  one  means  of  displacement 
at  node  point  9.  This  produces  a  more  redundant  structure  with 
a  resultant  decrease  in  total  volume  of  10.64%.  In  comparison. 
Case  #8  with  one  less  degree  of  freedom  than  Case  #7  at  both 
node  1  and  node  9,  had  a  reduction  in  total  volume  of  14.08%. 
Again,  it  is  suggested  that  a  structure  more  statically 
indeterminate  results  in  a  more  efficient  structure. 
Additionally,  it  is  noted  that  when  the  axial  stress  dominates 
the  total  stresses,  the  area  is  reduced  significantly  and  the 
cross  section  dimensions  reduce  to  form  a  square. 

Loads 

=  0  lbs 

=  16,  000  lbs 

=  0  in-lbs 

End  conditions 

0  K02 
1  KOO 

Dimensions 

=  32  in 

=  180  degrees 

Total  volume 

=  256.61  in’ 


1“ 

' - ►J 


t4O00^ 


Lateral 

Axia. 

Moment 


Node  1 
Node  9 


Radius 

Theta 


Volume 
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S.  CASE  #11:  HINGED-FIXED  ARCH  WITH  LATERAL  LOADING 
AND  MOMENT 

To  investigate  the  possibility  that  dominant  axial 
stresses  result  in  volume  reduction,  a  semicircular  arch  with 
one  degree  of  freedom  at  node  1  and  zero  degrees  of  freedom  at 
node  9  was  subjected  to  a  lateral  load  and  applied  bending 
moment.  From  this,  it  appears  that  the  cross  sectional  area  is 
inversely  proportional  to  the  axial  stresses.  Additionally,  it 
appears  that  the  dimension  and  stress  curves  of  the  left  half 
of  the  structure  behaves  exactly  as  those  of  Case  #7  which  has 
the  identical  end  conditions.  Similarly,  the  curves  of  the 
right  half  of  the  arch  behaves  exactly  as  those  of  Case  #8. 
This  suggests  that  the  boundary  conditions  do  not  effect  the 
structure  past  the  midpoint. 

Leads 

Lateral  «  12,0CC  lbs 


94 


Dimension  (in) 


Case  #1 1 


35 


3-  l-X 

I  \ 


¥ 

0- - - - - 1 

123456789 


Node 


'  Base  — ^  Height  Area  (in ''  2) 


100000 :r 


123456789 

Node 
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T.  CASE  #11A:  HINGED-FIXED  ARCH  WITH  LATERAL  LOADING 

AND  MCXfENT 

In  order  to  test  the  possibility  that  dominant  axial 
stresses  might  reduce  the  cross  section  area  and  hence  reduce 
the  total  structure  volume,  the  structure  of  Case  #11  was 
subject  to  the  same  bending  moment  while  the  lateral  load  was 
increased  by  a  factor  of  2.  As  a  result,  the  axial  stresses 
increased  overall.  Again,  it  appears  that  the  cross  sectional 
area  is  inversely  proportional  to  the  axial  stresses.  The 
dimension  and  stress  curves  displayed  the  same  shape  as  noted 
before  but  the  total  volume  increased  by  57.95%.  Therefore,  it 
was  concluded  that  increasing  axial  stresses  may  reduce  the 
cross  sectional  area  at  specific  nodes  but  the  overall 
structure  volume  is  not  reduced. 


l.oaas 

Lateral 

=  2<,C0C  lbs 

Axial 

=  C  lbs 

Moment 

=  1,000  in-ibs 

End  conditions 

Node  1 

1  MOO 

Node  9 

0  MOO 

Dimensions 

Radius 

=  32  in 

Theta 

=  180  degrees 

Total  volume 

Volume 

=  241.08  in’ 
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Case  #11a 


U.  CASE  #12:  FIXED-FIXED  ARCH  WITH  MULTIPLE  LOADING 

To  demonstrate  further  versatility  of  this  program,  a 
fixed-fixed  arch  was  subjected  to  a  combination  load  applied 
at  an  angle  45  degrees  up  from  node  9.  The  load  consisted  of 
a  concentrated  lateral  and  axial  load  and  an  applied  bending 
moment.  As  anticipated,  there  is  a  jump  in  the  dimension 
curves  at  node  7  were  the  load  was  applied.  Interestingly, 
22.5  degrees  from  each  endpoint,  the  axial  stress  dominates 
and  accordingly,  the  cross  sectional  area  reduces 
significantly.  Additionally,  at  node  5,  the  midpoint  of  the 
arch  structure,  the  cross  sectional  area  is  significantly 
smaller  as  a  result  of  an  increase  in  the  axial  stress. 


Loads 


Lateral 

17,000  lbs 

Axial 

9,  000  Ibr 

Hoff.er.t 

=  1,000  in-  bs 

End  conditions 

Node  1 

0  MOO 

Node  9 

0  MOD 

Dimensions 

Radius 

=  32  in 

Theta 

=  180  degrees 

Total  volume 

Volume 

=  156.55  in’ 
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Case  #12 


V.  CASE  #13:  FIXED-FIXED  ROLLER  ARCH  WITH  AXIAL  LOADING 
AND  MOMENT 

For  this  case,  the  same  structure  of  Case  #10  is  subjected 
to  an  equivalent  axial  load  with  an  additional  bending  moment 
applied  at  nodal  point  6,  67.5  degrees  up  from  node  9. 
Shifting  the  load  by  22.5  degrees  and  adding  the  applied 
bending  moment  appeared  to  have  little  effect  on  the  overall 
design.  In  fact,  the  volume  is  increased  from  Case  #10  by  only 
3.65%  and  the  dimension  curves  exhibit  very  similar 
characteristics.  However,  the  dip  observed  previously  in  the 
axial  stress  curve  at  node  6  of  Case  #10  is  not  present  in  the 
axial  stress  curve  of  this  case. 

Loads 

=  0  lbs 

=  16,000  lbs 

=  1,000  in-lbs 

End  conditions 

0  KOD 
1  MOD 

Dimensions 

=  32  in 

=  180  degrees 

Total  volume 

=  265.96  in’ 


L, 


I  KM 


Lateral 

Axial 

Moment 


Node  1 
Node  9 


Radius 

Theta 


Volume 
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W.  CASE  #13A:  FIXED-FIXED  ROLLER  ARCH  WITH  AXIAL  LOADING 

AND  MC»1ENT 

For  comparison,  the  same  structure  and  loading  were 
reoptimized  starting  from  a  different  initial  design  point. 
Previously,  each  optimization  began  from  an  initial  design  of 
2  inches  by  2  inches  at  each  node.  For  this  case,  the  base 
dimension  at  each  nodal  point  was  0.5  inches,  and  the  height 
dimension  at  each  nodal  point  was  3.5  inches.  Incredibly,  the 
volume  of  the  resultant  structure  is  33.96%  less  than  the 
volume  of  the  structure  optimized  in  Case  #13.  Obviously, 
optimization  can  be  a  function  of  the  starting  point. 
Fortunately  in  the  previous  cases,  various  initial  design 
points  were  tested  and  this  occurence  did  not  repeat  itself. 


Loads 


Lateral 

0  lbs 

Axial 

-  16,000  lbs 

Moment 

=  1,000  in-lbs 

End  conditions 

Node  1 

0  MOD 

Node  9 

1  MOD 

Dimensions 

Radius 

=  32  in 

Theta 

=  180  degrees 

Total  volume 

Volume 

=  175.65  in’ 
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Case  #13a 


X.  CASE  #14:  FIXED-HINGED  ARCH  WITH  MULTIPLE  LOADING 

In  contrast  to  Case  #13,  this  asynunetric  arch  was  loaded 
at  an  angle  on  the  side  of  the  arch  with  zero  means  of 
displacement  at  the  endpoint.  The  behavior  exhibited  by  the 
dimension  and  stress  curves  was  similar  to  that  of  Case  #12 
which  had  zero  means  of  displacement  at  both  endpoints. 
Allowing  for  the  difference  in  the  magnitude  and  direction  of 
the  load,  the  only  significant  difference  between  Case  #12  and 
this  case  appears  at  node  8  and  9.  It  is  presumed  that  the 
added  means  of  displacement  at  node  9  caused  such  a 
difference.  To  ensure  that  a  true  optimum  had  been  reached, 
attempts  were  made  to  optimize  this  structure  for  several 
different  initial  starting  points.  Consistent  results  were  not 
obtained. 


Loads 


Lateral 

=  5,000 

lbs 

Axial 

=  9,000 

lbs 

Moment 

=  1,000 

in-lbs 

End  conditions 

Node  1 

0  MOD 

Node  9 

1  MOD 

Dimensions 

Radius 

=  32  in 

Theta 

=  180  degrees 

Total  volume 

Volume 

121.28 

in* 
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y.  CASE  #15:  FIXED-HINGED  ROLLER  ARCH  WITH  MULTIPLE  LOADING 

For  the  last  case  studied,  a  concentrated  lateral  load  and 
bending  moment  are  applied  at  the  midpoint  in  combination  with 
a  load  acting  radially  outward  distributed  along  the  length  of 
the  arch.  The  cross  sectional  area  behaved  as  anticipated  from 
Case  #11,  inversely  proportional  to  the  axial  stress.  Again, 
to  ensure  that  a  true  optimum  had  been  reached,  attempts  were 
made  to  optimize  this  structure  for  several  different  initial 
starting  points.  Consistent  results  were  not  obtained. 

Loads 

Lateral  =  16,000  lbs 

Axial  =  0  lbs 

Moment  =  1,000  in-lbs 

Distrib.  =  100  Ibs/in. 

End  conditions 

Node  1  0  MOD 

Node  9  2MOD 

Dimensions 

Radius  =  32  in 

Theta  =  180  degrees 

Total  volume 

Volume  =  285.75  in’ 


106 


>-  Base 


Height 


Z.  CASE  #15A:  FIXED-HINGED  ROLLER  ARCH  WITH  MULTIPLE  LOADING 

For  comparison,  the  structure  of  Case  #15  was  subject  to 
the  same  bending  moment  and  distributed  load  while  the  lateral 
load  was  doubled  in  value.  As  demonstrated  by  Case  #lla,  an 
overall  increase  in  the  axial  stresses  results  does  not  effect 
the  shape  of  the  dimension  and  stress  curves.  However,  the 
volume  from  Case  #15  is  increased  by  79.65%.  In  comparison, 
doubling  the  lateral  load  for  Case  #11  resulted  in  an  increase 
in  volume  of  57.95%.  Of  interest,  it  appears  that  the  majority 
of  the  volume  increase  is  centered  around  the  midpoint  were 
the  increased  load  was  applied. 


Loads 


108 


VII.  CONCLDSIOKS 


The  conclusions  of  this  study  are  as  follows: 


•  The  bar-beam  model  for  stress  analysis  yielded  results 
which  deviated  from  known  analytical  solutions  with  an 
error  of  less  than  2%.  Therefore,  the  technique  of 
modeling  arch  structures  with  bar-beam  elements  is  deemed 
a  viable  approximation.  (Chapter  V) 

•  From  the  specific  cases  studied,  the  Sequential  Linear 
Programming  method,  (Method  2) ,  best  performed  the 
optimization  for  cantilever  arch  structures.  The  Modified 
Method  of  Feasible  Directions  (Method  1)  best  performed 
the  optimization  for  arch  structures  with  restrictive 
boundary  conditions  at  both  endpoints.  (Case  #2,  2a,  1) 

•  Reoptimization  of  an  optimal  solution  has  the  effect  of 
smoothing  the  results  and  reducing  the  volume  of  the 
structure.  The  effect  of  this  two-stage  optimization 
strategy  was  more  significant  for  Method  2  than  Method  1. 
(Case  #3a,  7a) 

•  The  DOT  auto  scaling  function  inhibited  the  optimizer 
performance.  (Case  #6a) 

•  Applying  moments  that  produce  stresses  that  oppose  the 
stresses  produced  by  a  concentrated  load  reduce  the  total 
structure  volume  required  to  withstand  the  combined  load. 
Through  prestressing  one-way  loaded  structures,  more 
efficient  structures  can  be  achieved.  (Case  #4,  4a,  6) 

•  The  cross  sectional  shape  is  dependant  on  the  type  of 
stress  experienced.  When  bending  stresses  dominate,  the 
optimal  cross  section  forms  a  tall  rectangle  limited  only 
by  the  geometric  constraint.  When  axial  stresses  dominate, 
the  optimal  cross  section  dimensions  form,  a  square.  (Case 
#5,  7,  10) 

•  Structures  which  are  more  statically  indeterminate  are 
more  efficient  under  identical  loading  than  less  redundant 
structures.  (Case  #8,  10) 

•  Asymmetric  structures  are  more  likely  to  produce  erroneous 
results  due  to  the  limit  of  the  number  of  elements  used  to 
obtain  results.  (Case  #9) 
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•  The  boundary  conditions  act  as  an  excitation  which  follow 
the  St.  Venant  principle.  The  information  from  the 
boundary  condition  diminishes  such  that  for  a  semicircular 
'  arch,  the  boundary  conditions  do  not  effect  the  cross 

sectional  shape  past  the  arch  midpoint.  (Case  #11,  14) 

T  •  Optimization  is  a  function  of  the  initial  design  starting 

point.  (Case  #13a) 


From  this  investigation,  the  following  is  suggested  as  a 
possibility  for  future  research  in  the  realm  of  weight  optimum 
arch  structures: 


•  Continue  to  record  results  for  a  comprehensive  study  of 
all  combinations  of  parameters,  loadings,  and  end  condi¬ 
tions 

•  Optimize  the  arch  structure  using  varied  cross  sections 
such  as  a  C,  L,  or  I  beam,  a  box  beam,  or  a  circular  beam. 

•  Remove  the  assumption  that  the  arch  maintains  a  constant 
radius  of  curvature  and  optimize  the  arch  shape. 

•  Apply  additional  constraints  such  as  global  buckling  in 
order  to  present  a  more  accurate  model. 


Ill 


APPENDIX  A 

DOT  PROGRAM  PARAMETERS 


The  information  in  the  following  tables  is  taken  from 
[Ref.  7] 


SCALAR  PARAMETERS  STORED  IN  RPRM 


Ill  1. 

LOCATION 

NAME 

DEFAULT  VALUE 

RPRM ( 1 ) 

CT 

-0.05 

RPRM (2) 

CTMIN 

0.003 

RPRM (3) 

DABOBJ 

MAX[0. 001*ABS (FO) ,0.0001] 

RPRM ( 4 ) 

DELOBJ 

0.001 

RPRM (5) 

DOBJl 

0.1 

RPRM (6) 

DOBBJ2 

0.2*ABS  (FO) 

RPRM ( 7 ) 

DXl 

0.01 

RPRM ( 8 ) 

DX2 

0.2*,AX[X(1)  ] 

RPRM  (9) 

FDCH 

0.001 

RPRM (10) 

FDCHM 

0.0001 

RPRM (11) 

RMVLMZ 

0.4 

RPRM(12) 

DABSTR 

MAX[0.001*ABS (FO) ,0.00001] 

RPRM (13) 

DELSTR 

0.001 

1  RPRM (14) -RPRM (20)  RESERVED  FOR  INTERNAL  USE 

NOTE:  FO  =  The  value  of  the  objective  function  at  the  start  of 
optimization  (for  the  initial  values  of  X) . 


LOC. 


PARAM 


CT 


CTMIN 


DABOBJ 


DELOBJ 


DOBJl 


DOBJ2 


DXl 


DX2 


FDCH 
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FDCHM 


11 


RMVLMZ 


12 


DABSTR 


DEFINITIONS  OF  PARAMETERS  CONTAINED 
IN  THE  RPRM  ARRAY 


DEFINITION 

A  constraint  is  active  if  its  numerical  value  is 
more  positive  than  CT.  CT  is  a  small  negative  niimber 

A  constraint  is  violated  if  its  numerical  value  is 
more  positive  than  CTMIN 

Maximum  absolute  change  in  the  objective  between 

ITRMOP  consecutive  iterations  to  indicate 
convergence  in  optimization 

Maximum  relative  change  in  the  objective  between 

ITRMOP  consecutive  iterations  to  indicate 
convergence  in  optimization 

Relative  change  in  the  objective  function  attempted 
on  the  first  optimization  iteration.  Used  to 
estimate  initial  move  in  the  one-dimensional  search. 
Updated  as  the  optimization  progresses. 

Absolute  change  in  the  objective  function  attempted 
on  the  first  optimization  iteration 

Maximum  relative  change  in  a  design  variable 
attempted  on  the  first  optimization  iteration.  Used 
to  estimate  the  initial  move  in  the  one-dimensional 
search.  Updated  as  the  optimization  progresses 

Maximum  absolute  change  in  a  design  variable 
attempted  on  the  first  optimization  iteration.  Used 
to  estimate  the  initial  move  inthe  one-dimensional 
search.  Updated  as  the  optimization  progresses. 

Relative  finite  difference  step  when  calculating 
gradients 

Minimum  absolute  value  of  the  finite  difference  step 
when  calculating  gradients.  This  prevents  too  small 
a  step  when  X(l)  is  near  zero 

Maximum  relative  change  in  design  variable  during 
the  first  approximate  subproblem  in  the  Sequential 
Linear  Programming  Method.  This  is,  each  design 
variable  is  initially  allowed  to  change  by  ±40%. 

This  move  limit  is  reduced  as  the  optimization 
progresses . 

Maximum  absolute  change  in  the  objective  between 
itrmst  consecutive  iterations  of  the  Sequential 

Linear  Programming  method  to  indicate  convergence  to 
the  optimum 

Maximum  relative  change  in  the  objective  between 

ITRMST  consecutive  iterations  of  the  Sequental 

Linear  Programming  method  to  indicate  convergence  to 
the  optimum 
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DELSTR 


PARAMETERS  IN  THE  IPRM  ARRAY 


LOCATION 

NAME 

DEFAULT  VALUE 

IPRMd) 

I  GRAD 

0 

IPRM(2) 

IS  CAL 

NDV 

IPRM(3) 

ITMAX 

40 

IPRM(4) 

ITRMOP 

2 

IPRM(5) 

IWRITE 

6 

IPRM(6) 

NCOLA 

NCON+NDV,  but  at  least  2*NDV  and 
not  more  than  10*NDV 

IPRM(7) 

IGMAX 

0 

IPRM{8) 

JTMAX 

20 

IPRM(9) 

ITRMST 

2 

IPRM(IO) 

JPRINT 

0 

IPRM(ll) 

IPRNTl 

0 

IPRM(12) 

IPRNT2 

0 

IPRM (13 

JWRITE 

0 

1  IPRM (14) -IPRM (18) 

RESERVED  FOR  FUTURE  USE 

IPRM(19) 

NEWITR 

INTERNALLY  DEFINED 

IPRM(20) 

NGT 

INTERNALLY  DEFINED 
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DEFINITIONS  OF  PARAMETERS  CONTAINED  IN  THE  IPRM  ARRAY 


PARAX . 


ic-ra:: 


:  TMAX 


ITRMC? 


NCOLA 


IGKAX 


:?rnt: 


:?RN':2 


PFINiriON 


Specifies  wnecr.er  'Che  gradients  a 
(IGRA3=G)  or  by  the  user  (IGRAD=i 


iateo  by  DO. 


Design  variables  are  rescaled  every  ISCAL  iterations.  Set 
rSCAL=-l  to  turn  off  scaling 


Maximur  number  of  iterations  allowed  at  the  ootimize  level 


The  numiber  of  consecutive  iterations  for  which  the  absolute 
or  relative  convergence  criteria  must  be  met  to  indicate 
convergence  at  the  optimizer  level 


If  IGKAX=C,  only  gradients  of  active  and  violated 
constraints  are  calucated.  If  IGMAX>C,  up  to  NCOLA 
gradients  are  calculated,  including  active,  violated,  and 
near  active  constraints 


Maximum  num.ber  of  iterations  allowed  for  the  Sequential 
Linear  Program.m.ing  method.  This  is  the  numbier  of  linearized 
super obiem.s  solved. 


iT.e  num.ber  of  consecutive 
or  relative  convergence  cr 
convercence  in  the  Seouent 


solute 

ate 

ca 


Sequential  Linear  Programming  subproblem  print.  If 
JPRINT>C,  IPRINT  is  turned  on  during  approxim.ate  linear 
subproblem..  This  is  for  debugging  only 


:prnti=i. 


sea 

g  f ac 

r  the  X  vector 


If  I?RNT2  =  1,  print  m.iscel  ianeous  search  inf  orm.ation .  If 
I?RNT2=2,  turn  on  print  during  one-dimensional  search 
process.  This  is  for  debugging  only 


File  num.ber  to  write  iteration  history  inform, ation  to.  This 
is  useful  for  using  postprocessing  program,  to  plot  the 
Iteration  process.  This  is  only  used  if  JWRITF>C 


Normally  =-l.  Set  =n  at  the  start  of  a  new  iteration,  where 
n  is  the  nu.mlaer  of  the  iteration  just  co.mpleted .  If 
XET.4CD=C,1,  this  is  after  each  one-dimensional  seaarch.  If 
KET.nCD=2,  this  is  after  eac.h  approxim.ate  optim.i zation  If 
JWRITE>0,  the  optimization  inform.ation  will  have  just  been 
written  to  that  file.  If  you  with  to  stop  after  each 
iternation  (or  after  a  particular  iteration)  and  then  re¬ 
start  later,  NEWITR  is  a  flag  to  do  this.  NEWITR  is  defined 
internally  by  DOT 


Tne  num.ber  of  constraint  gradients  neeaed.  If  the  user 
supplies  gradients  to  DOT,  this  will  be  needed.  The 
constraint  numbers  for  which  gradients  are  needed  are 
contained  in  positiooon  1-NGT  of  the  IWK  array.  NGT  is 
defined  internally  by  DOT 


APPENDIX  B 

JUSTIFICATION  FOR  <»(ITTIN6  SHEAR  STRESSES 

(The  following  Appendix  is  taken  from  [Ref.  9]) 

The  shear  stress  distribution  through  a  beam  of 
rectangular  cross-section  has  a  parabolic  distribution  along 
the  height  of  the  member.  The  maximum  shear  stress,  located  at 
the  neutral  axis  of  the  beam,  is 

T,,,  =  1.5V/ A  (B.l) 

where  is  the  maximum  shear  stress,  V  is  the  shear  force, 

and  A  is  the  cross-sectional  area  of  the  beam.  [Ref.  8] 

The  normal  stress  due  to  bending  is  given  by  the  equation 

=  Mc/i  (B.2) 

where  a.  is  the  maximum  normal  stress,  M  is  the  bending 
moment,  and  I  is  the  cross-sectional  moment  of  inertia  which 
for  this  case  is  bh^/12  where  b  and  h  are  the  width  and  height 
respectively  of  the  cross-section. 

Redefining  the  normal  stress  in  terms  of  the  cross- 
sectional  dimensions  yields 

=  M(h/2)  /  (bhVl2) 
or 

O,.  =  6M/hA  (B.3) 

The  ratio  of  the  maximum  shear  stress  to  the  normal  stress 
due  to  bending,  is  denoted  by  r  and  given  by  the  expression: 

r  =  (B.4) 
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Substituting  Equations  (B.l)  and  (B.3)  into  Equation  (B.4) 
yields 

r  =  (1.5V/A)  /  (6M/hA) 


or 


r  =  Vh/4M  (B.5) 

For  the  cases  investigated  in  this  study,  the  maximum  value  r 
can  attain  is  when  the  loading  is  that  of  a  uniformly 
distributed  load,  Py.  Then,  where: 

V  =  PyL  (B.6) 

M  =  PyLV2  (B.7) 

which  upon  substitution  into  Equation  (B.8)  yields 

r  =  (PyL)h/4  (PyLV2) 

which  simplifies  to 

r  =  h/2L  (B.8) 

The  use  of  the  beam  equation  requires  the  length  of  the 
beam  to  be  at  a  minimum  ten  times  the  height,  that  is: 

L  >  lOh  (B.9) 

To  maximize  the  value  of  r,  let  L  equal  lOh,  the  minimum 
allowable  length.  Substitutiing  this  value  of  L  into  Equation 
(B.8)  yields 

r  ^  h/2(10h) 

or  simply 


r  <  1/20  (B.IO) 
Hence,  the  maximum  shear  stress  accounts  for  less  than  5%  of 
the  bending  stress  developed  in  the  structure.  Five  percent  is 
high  considering  this  analysis  over-assumed  the  value  of  the 
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shear  stress  by  assigning  the  maximum  shear  stress  to  the 
entire  cross-section  of  the  beam.  Moreover,  at  the  outermost 
fibers  where  o.  is  a  maximum,  the  shear  stress  is  zero. 
Therefore,  under  the  circumstances  of  this  study,  the  addition 
of  shear  stresses  was  deemed  to  be  unwarranted. 
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APPENDIX  C 

archjoptimization  computer  code 

PROGRAM  ARCH  OPTIMIZATION 


ARCH  OPTIMIZATION  ANALYSIS  CODE  * 


ALPHA - TRANSFORMATION  ANGLE  OF  ELEMENT  (ANGLE  TO  X-AXIS) 

ANGLE. .. .TOTAL  ANGLE  OF  ARCH  (IN  DEGREES) 


SAVE . THE  AVERAGE  BASE  DIMENSION  ACROSS  AN  ELEMENT 

BASE . DOT  ARRAY  CONTAINING  THE  ELEMENTAL  BASE  DIMENSIONS 


BASEL.... DOT  ARRAY  CONTAINING  THE  ELEMENTAL  BASE  DIMENSIONS  LOWER 
SIDE  CONSTRAINT 

BASEU....DOT  ARRAY  CONTAINING  THE  ELEMENTAL  BASE  DIMENSIONS  UPPER 
SIDE  CONSTRAINT 

BETA _ TRANSFORMATION  ANGLE  OF  ELEMENT  (ANGLE  TO  Y-AXIS) 


B  1 . BOUNDARY  TERMS  APPLIED  AT  END  ”1" 

B~2 . BOUNDARY  TERMS  APPLIED  AT  END  "2" 

CT,  .. ,C5. CONSTANTS  RELATED  TO  ELEMENT  STIFFNESS  COEFFICIENTS 

CLAN . CONCENTRATED  LOAD  APPLICATION  NODE  (THE  NODE  FX.FY.FM  ARE 

APPLIED 

COUNT ....  COUNTS  THE  NUMBER  OF  ITERATIONS  COMPLETED 

DOF . DEGREE  OF  FREEDOMS  (UNKNOWN  DISPLACEMENTS  &  SLOPES) 

DSN . DESIGN  VARIABLE  FOR  EACH  ELEMENT 

DESIGN... DOT  ARRAY  CONTAINING  THE  ELEMENTAL  BASE  AND  HEIGHT  DIMENSIONS 


OESIGNL..OOT  ARRAY  CONTAINING  THE  ELEMENTAL  BASE  AND  HEIGHT  DIMENSIONS 
LOWER  SIDE  CONTRAINT 

DESIGNU..DOT  ARRAY  CONTAINING  THE  ELEMENTAL  BASE  AMD  HEIGHT  DIMENSIONS 
UPPER  SIDE  CONSTRAINT 

DVIBG. .. .DESIGN  VARIABLE  *1  (BASE  DIMENSION)  INITIAL  ESTIMATE 
DVILO. .. .DESIGN  VARIABLE  tl  (BASE  DIMENSION)  LOWER  SIDE  CONSTRAINT 
DVIUP. .. .DESIGN  VARIABLE  #1  (BASE  DIMENSION)  UPPER  SIDE  CONSTRAINT 
OV2BG. .. .DESIGN  VARIABLE  12  (HEIGHT  DIMENSION)  INITIAL  ESTIMATE 
OV2LO. .. .DESIGN  VARIABLE  #2  (HEIGHT  DIMENSION)  LOWER  SIDE  CONSTRAINT 
OV2UP. .. .DESIGN  VARIABLE  #2  (HEIGHT  DIMENSION)  UPPER  SIDE  CONSTRAINT 


EK . 6X6  ELEMENT  STIFFNESS  MATRIX  IN  LOCAL  X,Y  COORDINATES 

EKPR . 6X6  ELEMENT  STIFFNESS  MATRIX  IN  ELEMENT  LOCAL  COORDINATES 

ELEN . LENGTH  OF  ELEMENT 

FA . CONSTANT  DISTRIBUTED  LOAD  OUTWARD  FROM  END  TO  END 

FH . CONCENTRATED  MOMENT  AT  FREE  END 

FX . CONCENTRATED  LOAD  IN  X  DIRECTION  AT  FREE  END 

FY . CONCENTRATED  LOAD  IN  Y  DIRECTION  AT  FREE  END 

G . THE  ARRAY  OF  CONSTRAINT  FUNCTIONS 

GAMMA.... 6X6  ELEMENT  TRANSFORMATION  MATRIX 

GR . (NDOF)X(NOOF)  GLOBAL  STIFFNESS  MATRIX 

RAVE . THE  AVERAGE  HEIGHT  DIN.  ACROSS  THE  ELEMENT 

HGT . DOT  ARRAY  CONTAINING  THE  ELEMENTAL  HEIGHT  DIMENSIONS 

HGTL . DOT  ARRAY  CONTAINING  THE  ELEMENTAL  HEIGHT  DIMENSIONS 

LOWER  SIDE  CONSTRAINT 

HGTU . DOT  ARRAY  CONATINING  THE  ELEMENTAL  HEIGHT  DIMENSIONS 

UPPER  SIDE  CONSTAINT 

INOSN. .. .INITIAL  (UNIFORM)  DESIGN  DIMENSION 

INFO . DOT  PARAMETER  USED  TO  SIGNAL  THAT  THE  OPT  IS  COMPLETE 


IPRINT...OOT  PARAMETER  USED  SELECT  THE  DATA  OUTPUT  FORMAT 
IPRM . DOT  SELECTABLE  INTEGER  PARAMETERS 

ITERATE.. THE  NUMBER  OF  TIMES  DOT  IS  TO  BE  RELOADED  WITH  THE 
PRECEEDING  DATA 

IWK . DOT  INTERNAL  WORK  SPACE  ARRAY 
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oono  noon  no  on 


METHOD... DOT  PARAMETER  USED  TO  DEFINE  THE  OPTIMIZATION  METHOD 
MINMAX...DOT  PARAMETER  USED  TO  MINIMIZE/MAXIMIZE  THE  PROBLEM 

NOON . NUMBER  OF  DESIGN  CONSTRAINTS 

NOOF . NUMBER  OP  DEGREES  OF  FREEDOM 

NDV . NUMBER  OF  DESIGN  VARIABLES 

NEL . NUMBER  OP  ELEMENTS 

NRIWK....DOT  INTERNAL  WORK  SPACE  ARRAY  DIMENSION 

NRWK . DOT  INTERNAL  WORK  SPACE  ARRAY  DIMENSION 

NSNP . NUMBER  OF  SYSTEM  NODAL  POINTS 

OBJ . THE  OBJECTIVE  FUNCTION  OF  THE  OPTIMIZATION 

OPTDCS. . .OPTIMIZATION  DECISION  TO  OPTIMIZE  THE  PROBLEM  OR  NOT 
PI. . .PS. .PARAMETER  DIMENSION  CORRESPONDING  TO  THE  NEL,  NSNP,  NCON, 
NDOF,  AND  NDV  RESPECTIVELY 

PHI . SUBTENDED  ELELENT  ANGLE  (ALSO,  PHIANG  IN  DEGREES) 

PRCSN....THE  PRECISION  DESIRED  TO  SOLVE  THE  PEN  SYSTEM  OP  EQUATIONS 
RADIUS. . .ARCH  RADIUS 

RPRM . DOT  SELECTABLE  REAL  PARAMETERS 

SIGMA  B..THE  ELEMENTAL  NORMAL  STRESS  DUE  TO  BENDING 
SIGNA~N..THE  ELEMENTAL  NORMAL  STRESS  DUE  TO  AXIAL  FORCES 
SIGMA  T..THE  MAXIMUM  TOTAL  STRESS  IN  EACH  ELEMENT 

U _ T...THE  "DISPLACEMENT"  VECTOR  OF  THE  SYSTEM  OF  LINEAR  EQUATIONS 

WK . DOT  INTERNAL  WORK  AREA 

X . GLOBAL  HORIZONTAL  COORDINATE 

Y . GLOBAL  VERTICAL  COORDINATE 

YIELD _ YIELD  STRENGTH  OF  THE  ARCH  MATERIAL 

YOUNG ....  YOUNG ' S  MODULUS  OF  THE  ARCH  MATERIAL 


C  ....declare  the  variables . 

INCLUDE  'ARCH^COM.FOR' 

....read  the  input  parameters . 

0PEN(8,  FILE-'ARCH_IN.DAT' ,  STATUS- ' OLD ' ) 

READ ( 8 , • )  ANGLE , RADIUS , YOUNG , YI ELD , NEL , METHOD , I PRINT , DVl BG , 

S  DVILO , DVIUP , DV2BG , DV2LO , DV2UP , CLAN , PX , FY , FM , FA, OPTDCS , 

S  ITERATE , PRCSN , BXl , B Yl , BMl , BX2 , BY2 , BH2 , LABEL 

....define  constants . 

NSNP  •  NEL  *  1 
NDOF  -  3*NSNP 
NCON  -  3 "NSNP 
NDV  -  2*NSNP 

....determine  the  system  nodal  coord  and  element  orientation.... 
CALL  GEOMETRY  (NEL, NSNP, ANGLE, RADIUS, X,Y, ALPHA, BETA, ELEN) 

....define  the  size  of  the  work  arrays  for  DOT . 

NRWK  -  38800 
NRIWK  -  1000 

....optimize  the  problem . 

CALL  OPTIMIZATION_TOOL 

....compile  and  format  the  output . 

CALL  ARCH  OUTPUT 


C 

END 


SUBROUTINE  GEOMETRY  (NEL, NSNP. ANGLE, RADIUS, X,Y, ALPHA, BETA, ELEN) 
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I  This  routine  is  used  by  main  ARCH_OPTiniZATZON  to  generate  | 
I  the  X-,  y-coordinates  of  each  system  node,  to  determine  j 
I  the  orientation  of  each  element,  and  to  calculate  the  | 
I  length  of  each  element.  | 


....declare  the  variables . 

INTEGER  NEL,NSNP,P1,P2 

PARAMETER  ( Pl-32 , P2-33 ) 

REAL  ANGLE, RADIUS, ELEN, X(P2) ,Y(P2) ,ALPHA(Pl) ,BETA(Pl) , 
6  PI,PHI,ANG,yNUM,XOEN 

PARAMETER  (PI-3.141593) 


....determine  the  geometric  constants 
PHI  -  (ANGLE/NEL)*( PI/18 0.0) 


X(l)  -  0.0 
Yd)  -  0.0 

ANG  -  0.0 


100 


DO  100  i-l, 
ANG 
X{i*l) 
Y(i+1) 
YNUM 
XOEN 
ALPHA( i) 
BETA(i) 
CONTINUE 


NEL 

-  ANG  4-  PHI 

-  RADIUS  *  (1.0  -  COS(ANG)) 

-  RADIUS  •  SIN(ANG) 

-  (Y(ifl)  -  Y(i)) 

•  (X(i-t>l)  -  X(i) ) 

-  ATAN2( YNUM.XDEN) 

-  (PI/2.0)  -  ALPHA(i) 


....determine  the  length  of  each  element 
ELEN  -  SQRT(X(2)**2.0  +  Y(2)**2.0) 

RETURN 

END 


* 

SUBROUTINE  OPTIMIZATION  TOOL 


This  subroutine  directs  the  program  flow  optimization  decision 
l.e.,  optimize  the  problem  or  not.  It  also  serves  to  set  up  fi 
execute  the  DOT  optimization  software. 


....declare  the  variables . 

INCLUDE  'ARCH_COH.rOR' 

INTEGER  i 

....zero  out  the  RPRM  and  IPRM  arrays 
DO  100  i*l,20 
RPRM(i)  -  0.0 
IPRM(i)  -  0 
100  CONTINUE 

_ initialize  COUNT . 

COUNT  -  1 
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....refine  the  constraint  tolerence 
RPRM(2)  -  0.0001 
RPRMO)  -  0.001 


....turn  off  DOT'S  auto  scaling . 

IPRM(2)  «  -1 

....increase  DOT'S  default  nuotber  of  iterations . 

IPRrt( 3 )  •  1000 
IPRM(8)  •  1000 

....increase  DOT'S  nuniber  of  consecutive  convergence  criteria. 
IPRM(4)  -  3 
IPR«(9)  -  3 

....define  KZNHAX*-!  to  miniisize  the  objective  function . 

niNHAX  -  -1 

....initialize  the  design  variable  linits  and  best  guess . 

00  200  i-l,NSNP 

BASC<i)  -  DVIBG 
BASEL(i)  -  DVILO 
BASEO(i)  -  DVIOP 
HGT(i)  -  OV2BG 
HGTL(i)  -  DV2LO 
HGTO(i)  •  DV2UP 
200  CONTINUE 

....conbine  base  and  HGT  arrays  into  design  array . 

00  250  i-l,NSNP 
j-NSNP+i 

DESIGN(i)  •  BASE(i) 

DESIGNL( i)  -  BASEL( i ) 

0ESIGNU( i )  -  BASEU( i ) 

DeSIGS(j)  -  HGT(i) 

OESIGNL(  j  )  -  RGTl.(  i  ) 

DESIGNU( j )  -  HGTU( i ) 

250  CONTINUE 

....■ake  optimization  decision . 

IP  (OPTDCS  .NE.  1)  THEN 
CALL  EVAL 
RETURN 
END  IP 

....ready  to  optimize . 

INPO  -  0 

300  CALL  DOT  ( INPO, NETHOO, IPRINT, NOV, NCON. DESIGN. DESIGNL.DESIGNU, 

4  OBJ , NINNAX. G, RPRH , I PRN , WK , NRHR , INK , NRINK ) 

....evaluate  the  objective  function  and  constraints . 

IP  (INPO  .GT.  0)  THEN 
CALL  EVAL 
GOTO  300 
END  IP 
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....refine  the  solution  vector  by  reoptimising 
IF  (COUNT  .LT.  ITERATE)  THEN 
INFO  -  0 
COUNT  -  COUNT+1 
GOTO  300 
END  IF 

RETURN 

END 


SUBROUTINE  EVAL 


This  subroutine  is  used  to  evaluate  the  Objective  function, 
constraint  functions,  and  side  constraints  of  the  optimization 
problem. 


....declare  the  variables . 

INCLUDE  'ARCH_COM.FOR' 

INTEGER  i,j 

....separate  the  design  array  into  base  and  HGT  arrays 
DO  50  i-l,NSNP 

j-NSNP+i 

BASE(i)  «  DESIGN(i) 

HGT(i)  -  DESIGN! j) 

50  CONTINUE 
C 

C  . . . .calculate  the  objective  function . 

OBJ  -  0.0 
C 

DO  100  1-1, NEL 

BAVE(i)  -  (BASE!  DEBASE  (i-M)  )/2.0 
HAVE(i)  -  (HGT(  i)+HGT(i-*.l)  )/2.0 

OBJ  -  OBJ  *  BAVE( i) •HAVE! i) •ELEN 
100  CONTINUE 
C 

C  ....initialize  the  design  constraint  vector . 

DO  200  i-l,NCON 
G(i)  -  0.0 
200  CONTINUE 
C 

C  ....determine  the  de >ign  constraints . 

CALL  ARCH_STRESS 
C 

DO  210  i-l,NSNP 

j-i+NSNP 

k-i*( 2*NSNP) 

G(i)  -  (SIGKA  T(i)/YIELD)  -  1.0 
G(j)  -  BASE(iT-(3.0^HGT(l)) 

G(k)  -  HGT(i)-(10.0^BASE(i)) 

210  CONTINUE 
C 

RETURN 


END 


SUBROUTINE  ARCH  STRESS 
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This  subroutine  is  used  perform  the  Finite  Element  analysis 
of  the  stresses  developed  in  an  arch  or  beam  for  a  given  load¬ 
ing. 


....declare  the  variables . 

INCLUDE  ’ARCH_COM.FOR' 

INTEGER  I PVT (99) 

REAL  F(P4) 

REAL*8  BK(P4,P4) ,BF(P4),BU(P4),PAC( 9801),WORK( 99) 


....form  the  element  and  system  matrices . 

CALL  FORM 

....form  the  Force  vector,  F . 

CALL  rORCE_VECTOR  ( NEL , NDOF , ELEN, ALPHA, BETA, FA, F ) 

....set  the  boundary  conditiona  and  loads . . 

CALL  BNDARY  ( NDOF , GK , CLAN, FX , FY, FH, F, BXl , BYl , BNl , BX2 , BY2 4BM2 ) 

....solve  the  system  of  equations . . .' . 

IF  (PRCSN  .EQ.  2)  THEN 

....change  GR  and  F  arrays  to  double  precision . 

CALL  UPSCALE  ( NDOF , GK , F ,BK, BF ) 

....solve  the  system  of  equations . 

CALL  0L2ARG  (NDOF, BK,P4,Br,l.BU, FAC, IPVT, WORK) 

....change  BU  array  to  single  presicion . 

CALL  DOWNSCALE  (NDOF,BU,U) 

ELSE 

....solve  the  system  of  equations . 

CALL  L2ARG  ( NDOF , GK , P4 , F, 1 , U, FAC, IPVT, WORK ) 

END  IF 

....determine  the  stress  distribution . 

CALL  STRESS 

RETURN 

END 


SUBROUTINE  TORN 


This  subroutine  is  used  to  construct  the  global  stiffness  mat¬ 
rix  for  the  arch  problem. 


....declare  the  variables . 

INCLUDE  'ARCH_COM.FOR' 

INTEGER  IEL,I,J.K,II, JJ,KK,III,JJJ 

REAL  Cl,C2,C3,C4,C5,CA,CB,EK(Pl,6,6) ,GArtnA(6,6 ) ,EKGA(6,6 ) , 
a  GAEKGA(6,6) ,BH,BH3 


....define  the  constants  Cx 
Cl  -  YOUNG/ELEN 
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DO  220  I  -  1,6 
DO  215  J  -  1,6 
DO  210  K  -  1,6 

EKGA(I,J)  -  EKGA(I,J)  EKPR(  lEL ,  I ,  K )  *GAnnA(  K ,  J  ) 
210  CONTINUE 

215  CONTINUE 

220  CONTINUE 

....determine  the  GAEKGA  array . . . 

DO  240  I  -  1,6 
DO  235  J  -  1,6 
DO  230  K  -  1,6 

GAEKGA(I,J)  -  GAEKGAd,  J)-t-GAHnA(K,I)*ERGA(K,  J) 
230  CONTINUE 

235  CONTINUE 

240  CONTINUE 

....copy  the  GAEKGA  array  into  the  EK  array . 

DO  260  I  -  1,6 
DO  250  J  -  1,6 

EK(IEL,I,J)  -  GAEKGAd, J) 

250  CONTINUE 

260  CONTINUE 

120  CONTINUE 

....initialize  the  GK  array . 

DO  150  I  -  1,  NDOF 
DO  140  J  -  1,  NDOF 
GKd,J)  -  0.0 
140  CONTINUE 

150  CONTINUE 

....construct  the  GR  matrix . 

DO  300  lEL  -  1,  NEL 

II  -  3»dEL-l) 

DO  290  J  ■  1,  6 

JJ  -  II  J 
DO  280  R  -  1,  6 

KK  -  II  +  K 

GK(JJ,KK)  -  GK( JJ,KK)-)-EKdEL, J,K) 

280  CONTINUE 

290  CONTINUE 

300  CONTINUE 

RETURN 

END 


C 

c 

c 


SUBROUTINE  PORCE_VECTOR  ( NEL, NDOF, ELEN, ALPHA, BETA, FA, F ) 

This  subroutine  is  used  to  construct  the  force  vector  for  the 
FEM  problem  specified. 

....declare  the  variables . 

INTEGER  NEL, NDOF, i, 1 1,1 2, 1 3, PI, P4 

PARAMETER ( PI- 3 2 , P4-99  ) 

REAL  ELEN, ALPHA (PI) , BETA( PI ) , FA,  F( P4 ) 
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C2  -  (1.0/ELEN) **2.0 
C3  -  (1.0)/(2.0*ELEN) 
C4  -  (1.0)/3.0 
C5  -  (1.0)/6.0 


C 

C 


90 

100 

C 

C 


C 

c 


c 

c 


....initialize  the  work  arrays.. 

DO  120  lEL  -1,NEL 

DO  100  1-1,6 
DO  90  J-  1,6 

EKPR(1EL,I,J)  -  0.0 
GAMKA(I,J)  -  0.0 
EKGA(I,J)  -  0.0 
GAEKGA(I,J)  -  0.0 
EK(IEL,I,J)  -  0.0 
CONTINUE 
CONTINUE 


....calculate  the  area  and  inertia  terms 
BH  -  BAVE(IEL)*MAVE(IEL) 

BH3  -  BAVE(IEL)*(HAVE(IEL)**3.0) 


. . . .determine 
EKPR( ICL,l,i; 
ERPR( IEL,1,4) 
EKPR(IEL,2,2) 
EKPR(IBL,2,3) 
ERPR( ZEL,2,S} 
EKPR( IEL,2,6) 
EKPR(IEL,3,2) 
EKPR(IEL,3,3) 
BKPR(IEL,3,5) 
EKPR(IEL,3,6) 
EKPR(IEL,4,1) 
EKPR<ZEL,4,4) 
BKPR(IEL,5,2) 
EKPR( IEL,5, 3) 
BKPR( IEL,5,5) 
BKPR( IEL,S,6) 
EKPR( IEL,6,2) 
EKPR(IEL,6,3) 
EKPR(IEL,6,5) 
BXPR( IEL,6,6) 


the  EKPR  matrix 

-  C1*BB 

-  -C1*BH 

-  Cl*C2*BH3 

-  C1*C3*BH3 

-  -C1*C2*BH3 

-  C1*C3*BH3 

-  C1*C3*BH3 

-  C1*C4*BH3 

-  -C1*C3*BH3 

-  C1*C5*BH3 

-  -C1*BH 

-  C1*BH 

-  -C1*C2*BH3 

-  -C1*C3*BH3 

-  Cl*C2*BH3 

-  -C1*C3*BH3 

-  C1*C3*BH3 

-  C1*C5*BH3 

-  -C1*C3*BH3 

-  C1*C4*BH3 


....determine  the  GAHNA  matrix 
CA  -  COS (ALPHA (I EL ) ) 

CB  -  COS(BETA(IEL) } 


GAMMA(1,1) 

CA 

GAMKA(1,2) 

CB 

GAMHA(2,1) 

-CB 

GAnnA{2,2) 

CA 

GAKKA( 3,3) 

1.0 

GAI1KA(  4,4) 

CA 

GAW1A(  4,5) 

CB 

GAMHA(5,4) 

-CB 

GAW(A(S,5) 

CA 

GAIinA(6,6} 

1.0 

C 

C  ....determine  the  EKGA  array 
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c 


....form  the  F-vector . 

F(l)  -  (ELEN/2.0)  *  (-C0S(BETA(1) ) ) 

F(2)  -  (ELEN/2.0)  *  ( COS ( ALPaA( I ) ) ) 
r(3)  -  0.0 

c 

DO  100  i«2,NEL 

11  -  (i-l)*3  +  1 

12  -  (i-l)*3  +  2 

13  -  (i-l)*3  +  3 

C 

F(Il)  -  (ELEN/2. 0)*(-COS(BETA(i) ) ) 

&  +(ELEN/2.0)*(-COS(BETA(i-l) ) ) 

F(I2)  -  ( ELEN/2. 0) *(COS(ALPHA(i) ) ) 

C  +(ELEN/2.0)*(COS(ALPHA( i-1) ) ) 

F(I3)  -  0.0 
100  CONTINUE 
C 

F(NDOF-2)  -  (ELEN/2.0)*(-COS(BETA(NBL) ) ) 
F(NDOF-l)  -  ( ELEN/2. 0)*( COS (ALPHA( NED ) ) 
F ( NDOF )  -0.0 

....scale  the  F-vector  by  FA . 

DO  200  i-l,NDOF 
F(i)  -  FA*F(i) 

200  CONTINUE 

BETURN 
END 


SUBROUTINE  BNDARY  ( NDOP,GK,CLAN, PX, PY , PK, P,BX1 , BYl , BNl ,BX2 , 

£  BY2,BH2) 


This  subroutine  is  used  to  impose  the  boundary  conditions  upon 
the  global  stiffness  matrix  and  force  vector. 


....declare  the  variables . 

INTEGER  NDOP , BXl , BYl , BNl , BX2 , BY2 , Bn2 , CLAN ,i,N,Il,I2,l3,P4 

PARAnETER(P4-99) 

REAL  GK(F4,P4) ,PX.PY,Pn,P(P4) 

....invoke  the  essential  boundary  conditions . 

IP  (BXl  .EQ.  1)  THEN 

CALL  IHPOSE  BC  (NOOP,GR,l,P) 

ENDIP 
C 

IP  (BYl  .EO.  1)  THEN 

CALL  inPOSE_BC  (NOOP,GR,2,P} 

ENDIP 

C 

IP  (BNl  .EQ.  1)  THEN 

CALL  IHPOSE_BC  ( NDOP,GR, 3 , P) 

ENDIP 

C 

IP  (BX2  .EQ.  1}  THEN 
N-NDOP-2 

CALL  IHPOSE_BC  (NDOP,GK,N.r) 

END  IF 
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C 

IF  (BY2  .EQ.  1)  THEN 
N-NDOF-1 

CALL  1MP0SE_BC  (NDOF.GK.N, F) 

END  IF 
C 

IF  (BM2  .EQ.  1)  THEN 

CALL  IMPOSE  BC  ( NDOF.GK.NDOF, F) 

ENDIF 

....add  the  concentrated  load  to  the  force  vector . 

1 1- (  CLAN-1 

12-  (CLAN-1) *3+2 

13- ( CLAN-1) *3+3 

F{I1)-F(I1)+FX 
F( I2)-F( I2)+FY 
F(13)-F(I3)+FM 

RETURN 
END 

*•**••**+*••******•***********•************************* t 

SUBROUTINE  IMPOSE_BC  ( NDOF ,GK,N, F ) 


This  subroutine  is  used  to  do  the  redundant  leg  wor)c  of  impos¬ 
ing  the  boundary  conditions. 


....declare  the  variables . 

INTEGER  N00r,N,i,P4 

PARAHETER(P4-99) 

REAL  GK(P4,P4) ,F(P4) 

....impose  the  boundary  condition  on  the  GK  and  F  arrays 
DO  100  i-l,N00F 
GK(N,i)  -  0.0 
100  CONTINUE 

GK(N,N)  -  1.0 


RETURN 

END 


•  * 
SUBROUTINE  UPSCALE ( NDOF , GR , F , BK , BF ) 

C  This  subroutine  is  used  to  change  the  stiffness  matrix  (  force 

C  vector  from  single  precision  to  double  precision  in  order  to 

C  solve  the  linear  system  of  equations  in  double  precision. 

C 

C  ....declare  the  variables . 

INTEGER  NDOF,i,j,P4 

PARAMETER  (P4-99) 

REAL  GK(P4,P4) ,F(P4) 

REAL*8  BK(P4,P4) ,BF(P4) 
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....generate  the  doubleprecision  compliments  of  GK  and  r 
DO  110  i-l,NDOr 
DO  100  j-l,NDOr 

BK(i,j)  -  GK(i,j) 


100 

CONTINUE 

BF(i)  - 

110 

CONTINUE 

RETURN 

END 

* 

SUBROUTINE  DOWNSCALE (NOOF,BU,U) 


This  subroutine  is  used  to  do  dovn  scale  the  double  precision 
solution  of  the  linear  system  of  equations  back  to  single  pre¬ 
cision.  DOT  could  have  problems  with  double  precision  numbers! 


....declare  the  variables . 

INTEGER  NOOF,i,F4 

PARAMETER  (P4-99) 

REAL  U(P4) 

REALMS  BU(P4) 

....generate  the  doubleprecision  compliments  of  GK  and  F 
DO  100  i-l,NDOF 
U(i)  -  BU(i) 

100  CONTINUE 

RETURN 

END 


SUBROUTINE  STRESS 


This  subroutine  computes  the  stress  at  each  nodal  point. 

C  .............................................................. 

C  ....declarations . 

INCLUDE  'ARCH_COM.FOR' 

INTEGER  11,12,13,14,15,16 

REAL  CA1,CB1,K1,K2,FPR(P4,6),UPR(6) ,NORNl,NORM2, 

&  BEN01,BEND2 

C 

C  ....determine  local  forces  from  stiffness  and  displacement.... 

DO  100  i-1,  NEL 

11- (1-1)*3-^1 

12- (i-l)*3+2 

13- (i-l)*3+3 

14- (i)*3+l 

15- (i)*3+2 

16- < i)*3>3 
C 

CBl-  COS(BETA(i) ) 
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CM-  COS(ALPHA(  i  )  ) 

UPR(l)-  U(11)*CA1  +  U(I2)*CB1 

UPR(2)-  -U(I1)*CB1  U(I2)*CA1 

UPR(3)-  U(I3) 

UPR(4)-  U(I4)*CA1  +  U(I5)*CB1 

UPR(5)-  -U(I4)*CB1  +  U(I5)*CA1 

UPR(6)-  U(I6) 

DO  250  L-1,6 

rPR(i,L)-  0.0 
CONTINUE 

DO  300  J-1,6 
DO  350  K-1,6 

FPR(i,J)-  FPR(i,J)  +  EKPR(i,J,K)*UPR(K) 

CONTINUE 
CONTINUE 
CONTINUE 

....determine  the  bending  and  normal  streaaes . 

SIGMA_N(1)  -  ABS(FPR(1,1)*(1.0/(BASE(1)*HGT(1) ) ) ) 

SIGMA  B(l)  -  ABS(FPR(1,3)*(6.0/(BASE(1)*(HGT(1)**2.0)))). 
SIGKA”T(1)  -  SIGMA_B(1)  +  SIGMA_N(1) 

DO  400  i-2,NEL 

Kl  -  1.0/(BASE(i)*HGT(i)) 

K2  -  6.0/(BASE(i)*(HGT(i)**2.0)) 

NORMl  -  ABS(FPR(i,l)*Kl) 

NORM2  -  ABS(FPR(i-l,4)*Kl) 

BENDl  -  ABS(FPR(i,3)*K2) 

BEND2  -  ABS(FPR(i-l,6)*K2) 

SIGMA_N(i)  -  (NORM1+NORM2)/2.0 
SIGNA_B(i)  -  (BENOl-»-BENO2)/2.0 

SIGHA_T(i)  -  SIGMA_B(i)  +  SIGMA_N(i) 

400  CONTINUE 

SIGMA  N(NSNP}  -  ABS ( FPR( NEL, 4 ) * ( 1 . 0/( BASE< NSNP ) *HGT( NSNP ) ) ) ) 
SIGMA~B(NSNP)  -  ABS( FPR( NEL,6 ) * 

(6.0/(BASE(NSNP)*(HGT(N5NP)**2.0) ) ) ) 
SIGHA_T(NSNP)  -  SIGMA_B(NSNP)  SIGMA_N( NSNF ) 

RETURN 
END 


SUBROUTINE  ARCH  OUTPUT 


This  subroutine  formats  the  final  results  and  output  of  the 
optimization  problem  and  stores  it  in  a  file  named  ARCH_OUT.DAT 


C  ....declare  variables. 

INCLUDE  'ARCH  CON. FOR' 
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REAL  VOL, VOLUME 


....open  output  file  and  write  header . 

OPEN ( 9 ,  PILE-' ARCH_OUT . DAT ' ,  STATUS- ' UNKNOWN ' ) 

write(9,ioo)  label 

WRITE(9,100)  '  OPTIMIZATION  SOLUTION' 

WRITE(9,105)  '  - 

« - ' 

100  rORMAT(/5X,A) 

105  PORMAT(5X,A) 

_ .section  "A" . 

NRITE(9,100)  '  A)  Problem  Parameters:' 

WRITE(9,110)  '  Arch  Angle  ANGLE,  '  youngs  Modulus :' .YOUNG 

WR1TE(9,110)  '  Arch  Radius:',  RADIUS.  '  Yield  Strength; ' .YIELD 

WRITE(9,115)  '  No  of  Design  Var:',  NDV,  '  No  of  Elements :' ,NEL 
110  PORMAT(ax,A,P12.3,T38,A,Fl2.1) 

115  PORMAT(8X,A,I7,T38,A,I10) 

....section  "B" . . . . . 

WRITE(9,100)  '  B1  Derived  Constants:' 

WR1TE(9,120)  '  No  of  System  Nodal  Points. ..' ,NSNP 
HRITE(9,120i  '  No  of  Degrees  of  Freedom. ...' ,ND0P 


WRITE(9,125)  '  Length  per  Element . '.ELEN 

WRITB(9,125)  '  Phi  Angle  per  Element . '.PHIANG 

WRITE(9,120)  '  Number  of  Iterations . '.ITERATE 


120  FORnAT(8X.A,l6) 

125  FORnAT(8X,A.F12.4) 


....  section  "C* . 

WRITE(9,100)  '  C)  Structure  Loading:' 

WRITE(9,125)  'FX . 

WRITE(9,125)  'FY . 

WRITE(9,12S)  'FM . 

WHITE(9,125)  'FA . 


'  ,PX 
'  ,FY 
'  ,FM 
'  ,FA 


....section  "D" . 

WRITE(9.100)  '  D)  Elemental  Dimensions  and  Stress  Distribution:' 
WRITE(9,210)  'Node' , 'Height' , 'Base' , 'Length' , 'Area' 

210  fQRNAT(6X.A,T21,A,T36,A,T49,A,T62.A) 

220  FORMAT(8X,14,T17,F10.5,T32.F10.5,T48,F8.5,T60,F8.5) 

VOLUME  -0.0 


DO  300  i-l,NSNP 

AREA  -  BASE(i)*BGT(i) 

WRITE! 9, 220)  i,RGT( 1) .BASE! 1) .ELEN, AREA 
300  CONTINUE 


. . . .section  "B" . 

WRITB(9,100)  '  E)  Objective  Function;' 

C 

WRITE(9,310)  '  Total  structure  Volume: ' .OBJ 
310  FORMAT(/12X,A,F12.6/) 

C 

WRITE(9,330)  'Node' , 'Normal  Stress' , 'Bending  Stress 'Total ' 
DO  320  i-l,NSNP 
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WRITE(9,340)  i, SIGMA  N( i ) ,SIGMA_B( i ) , SIGMA_T( i ) 
320  CONTINUE 

33C  F0RMAT(8X,A,T18,A,T35,A,T57.A) 

340  rORMAT(8X,I4,T15,Fl4.1,T32,ri4.1,T49,F14.1) 

....  section  ”F" . 

WRITE(9,100)  '  F)  Boundary  Conditions:' 
WRITE{9,410)  'Node' , 'X-Displ' , 'Y-Displ' , 'Slope' 
WRITE(9,430)  1 , BXl , BYl , BMl 
WRITE(9,430)  NEL-:-!  ,  BX2  ,  BY2  ,  BM2 

. . . .section  "G" . 

WRITE(9,100)  '  G)  Solution  Vector:' 

WRITE(9,410)  'Node' , 'X-Displ' , 'Y-Oispl' , 'Slope' 

DO  400  i-l,NSNP 

11- (i-l)*3+l 

12- (i-l)*3-^2 

13- (i-l)*3+3 

WRITE(9,420)  i  ,U(  II) ,U< I2) ,U( I3) 

400  CONTINUE 

410  F0RMAT(T9,A,T17,A,T31,A,T46,A) 

420  F0RHAT(7X,I5,3E14.6) 

430  rORMAT(7X,15,T20,I4,T34,I4,T48,I4) 

RETURN 

END 
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ARCH  COMMON 


.  . .  .definitions . 

Pi . The  maximum  number  of  elements 

P2 . The  maximum  number  of  global  nodal  points 

P3 . The  maximum  number  of  design  constraints 

P4 . The  maximum  number  of  degrees  of  freedom 

p5 . The  maximum  number  of  design  variables 

....declare  the  variables . 

INTEGER  NEI. . NCON , NSNP , NOOF , NOV , METHOD, HINMAX , INFO , I PRINT , 

«  I WK ( 1 0  0  0 ) , NRWK , NRI WK , I PRM (20), COUNT , OPTDCS , I TERATE , 

S  PRCSN, CLAN , BXl , BYl , BHl , BX2 ,By2 , BM2 , PI , P2 . P3 , P4 , PS 

PARAMETER( Pl-32 , P2-3 3 , P3-96 , P4-99 , P5-64 ) 

REAL  ANGLE , RADIUS , ELEN, X( P2 ) ,¥( P2 ) , ALPHA( PI ) , BETA( PI ) , 
a  YOUNG, YI ELD, WK( 38800) ,RPRM( 20) ,OBJ,G( P3) , 

«  DV1BG,0V1L0,DV1UP,BASE(P1) ,BASEL(P1) ,BASEU(P1) , 

4  DV2BG,DV2LO,DV2UP,HGT(Pl) ,HGTL(P1) ,HGTU(P1) , 

4  DESIGN(P5) ,OESIGNL(PS) ,DESIGNU(P5) , 

4  FA, FX, FY, FM, U(P4) , SIGMA  T( P4 ) , SIGMA_N( P4 ) , 

4  SIGMA  B(P4) ,BAVE( PI) , 

4  HAVE(?1) ,GK(P4,P4) ,EKPR(P4,6,6) 

CHARACTER*30  LABEL 

....malce  in  common . 

COMMON  NEL , NCON , NSNP , NDOF , NDV , METHOD , M I NMAX ,INFO,IPRINT,IWK, 
4  NRWK , NRI WK , I PRM , COUNT , OPTDCS , I TERATE , PRCSN , CLAN , 

4  BX1,BY1,BN1,BX2,BY2,BH2, 

4  ANGLE, RADIUS, ELEN, X,Y, ALPHA, BETA, YOUNG, YIELD, 

4  WK , RPRM , OB J , G , DVl BG , DVl LO , OVl UP , BAS  E , BASEL , BAS  EU , 

4  OV2BG , OV2  LO , DV2UP , HGT , HGTL , HGTU , 

4  DESIGN, OESIGNL.DESIGNU, 

4  FA, FX,FY,FH,U,SIGHA_T, LABEL, SIGMA  N, SIGMA  B, 

4  BAVE,HAVE,GK,EKPR 
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APPENDIX  D 
VALIDATION  CASES 


validation  tl 
OPTIMIZATION  SOLUTION 


A)  Problea  Paraaatars: 

Arch  Angle  :  0.003 

Arch  Radius:  1000000.000 
No  of  Design  Var:  10 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus: 

30000000.0 

Yield  strength: 

52000.0 

No  of  Elements: 

4 

5 

15 

11.2500 

1 

C)  structure  Loading: 


FX .  1000.0000 

FY .  0.0000 

FM .  0.0000 

FA .  0.0000 


0)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

1 

3.00000 

1.50000 

11.24996 

2 

3.00000 

1.50000 

11.24996 

3 

3.00000 

1.50000 

11.24996 

4 

3.00000 

1.50000 

11.24996 

5 

3.00000 

1.50000 

11.24996 

Area 

4.50000 

4.50000 
4.50000 
4.50000 
4.50000 


E)  Objective  Function: 

Total  structure  Volume:  202.499207 


Node 

Normal  Stress 

Bending  Stress 

1 

0.0 

19999.7 

2 

0.0 

14999.7 

3 

0.0 

9999.8 

4 

0.0 

4999.9 

5 

0.0 

0.0 

Total 

19999.7 

14999.7 
9999. S 
4999.9 

0.0 


F)  Boundary  Conditions: 
Node  X-Oispl 

1  1 

5  0 


Y-Displ 

1 

0 


Slope 

1 

0 


G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOE-fOO 

2  0.257807B>01 

3  0.937478E-01 

4  0.189839E-K00 

5  0.2999938-^00 


Y-Displ 

O.OOOOOOE-i-00 

0.112327E-08 

0.409056E-08 

0.828721E-08 

0.130985E-07 


Slope 

O.OOOOOOE-fOO 

-0.437492E-02 

-0.749985E-02 

-0.937481E-02 

-0.999979E-02 
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Validation  #2 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  0.003 

Arch  Radius:  1000000.000 
No  of  Design  Var:  10 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

4 

5 

IS 

11.2500 

1 

C)  Structure  Loading: 


FX .  0.0000 

FY .  1000.0000 

FM .  0.0000 

FA .  0.0000 


D)  Elemental 

Dimensions 

and  Stress  : 

Distribution: 

Node 

Height 

Base 

Length 

1 

3.00000 

1.50000 

11.24996 

2 

3.00000 

1.50000 

11.24996 

3 

3.00000 

1.50000 

11.24996 

4 

3.00000 

1.50000 

11.24996 

5 

3.00000 

1.50000 

11.24996 

Area 

4.50000 

4.50000 

4.50000 

4.50000 

4.50000 


E)  Objective  Function: 

Total  structure  Volume:  202.499207 


Node 

Normal  Stress 

Bending  Stress 

1 

222.2 

0.0 

2 

222.2 

0.0 

3 

222.2 

0.0 

4 

222.2 

0.0 

5 

222.2 

0.0 

Total 

222.2 

22^.2 

222.2 

222.2 

222.2 


F)  Boundary  Conditions; 
Node  X-Displ 

1  1 

5  0 


Y-Displ 

1 

0 


Slope 

1 

0 


G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOEfOO 

2  0.112327E-06 

3  0.409056E-0e 

4  0.S28721E-08 

5  0.130985E-07 


Y-Displ 

O.OOOOOOE-fOO 

0.833330E-04 

0.166666E-03 

0.249999E-03 

0.333332E-03 


Slope 

O.O0OOOOE4-0O 

-0.191234E-09 

-0.327829E-09 

-0.409786E-09 

-0.437105E-09 
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Validation  13 
OPTIMIZATION  SOLUTION 


A)  Pcoblam  Parameters: 

Arch  Angle  :  0.003 

Arch  Radius:  1000000.000 
No  o£  Design  Var:  10 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 

rx . 

FY . 

FM . 

FA . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

4 

5 

15 

11.2500 

1 

0.0000 

0.0000 

10000.0000 

0.0000 

0)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

1 

3.00000 

1.50000 

11.24996 

2 

3.00000 

1.50000 

11.24996 

3 

3.00000 

1.50000 

11.24996 

4 

3.00000 

1.50000 

11.24996 

5 

3.00000 

1.50000 

11.24996 

Area 

4.50000 

4.50000 

4.50000 

4.50000 

4.50000 


E)  Objective  Function: 

Total  structure  Volume:  202.499207 
Node  Normal  Stress  Bending  Stress 


1 

0.0 

4444.4 

2 

0.0 

4444.4 

3 

0.0 

4444.4 

4 

0.0 

4444.4 

5 

0.0 

4444.4 

Total 

4444.4 

4444.4 

4444.4 

4444.4 

4444.4 


F)  Boundary  Conditions: 

Node  X-Oispl  Y-Displ 

11  1 
5  0  0 


Slope 

1 

0 


G)  Solution  Vector: 
Node  X-Displ 

1  0.000000E4-00 

2  -0.6249e5E-02 

3  -0.249994E-01 

4  -0.562488E-01 

5  -0.999979E-01 


Y-Displ 

O.OOOOOOE-t-00 

-0.273190E-09 

-0.109276E-08 

-0.245871E-08 

-0.437105E-08 


Slope 

O.OOOOOOE*00 

0.111109E-02 

0.22221SE-02 

0.333328E-02 

0.444438E-02 
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Validation  14 
OPTIMIZATION  SOLUTION 


A)  Problam  Pacametars: 


Arch  Angla  :  90.000 

Youngs  Modulus: 

30000000.0 

Arch  Radius:  45.000 

Yield  Strength: 

52000.0 

NO  of  Dasign  Var:  10 

No  < 

of  Elements: 

4 

B)  Oarivad  Constants: 

No  of  Systam  Nodal  Points 

•  a  a 

S 

No  of  Oagraas  of  Fraadom. 

a  a  a 

IS 

Length  par  Elanant . 

a  a  a 

17.5581 

Number  of  Iterations . 

a  a  a 

1 

C)  Structure  Loading: 

FX . 

0.0000 

FY . 

1000.0000 

FM . 

*  , 

0.0000 

FA . 

0.0000 

- 

D)  Elemental  Dimensions  and 

Stress 

Distribution: 

■* 

Node  Height 

Base 

Length 

Area 

1  3.00000 

1.50000 

17.55813 

4.50000 

2  3.00000 

1.50000 

17.55813 

4.50000 

3  3.00000 

l.SOOOO 

17.55813 

4.50000 

4  3.00000 

1.50000 

17.55813 

4.50000 

5  3.00000 

1.50000 

17.55813 

4.50000 

B>  Objaetiva  Function: 

Total  structuta  Voluna:  316.0463S6 


Node 

Normal  Stress 

Bending  Stress 

1 

217.9 

19996.8 

2 

201.4 

18475.5 

3 

154.1 

14141.1 

4 

83.4 

7653.6 

5 

43.3 

0.0 

Total 

20214.8 

18676.8 
14295.2 

7737.0 

43.3 


F)  Boundary  Conditions: 


Noda 
1 
5 


X-Displ 
1 
0 


Y-Displ 

1 

0 


Slopa 

1 

0 


G)  Solution  Vector: 

Node  X'Displ 

Y-Displ 

1  O.OOOOOOEaOO 

O.OOOOOOEaOO 

2  >0.654S28E>01 

0.131494E-01 

3  >0.223473Ea00 

0.118865Ea00 

4  -0.381495EaOO 

0.3S5493Ea00 

5  -0.446951E-^00 

0.684692Ea00 

Slopa 

O.OOOOOOE-^00 

0,750557E-02 

0.138688E-01 

0.181207E-01 

0.196139E>01 
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Validation  #4 
OPTIMIZATION  SOLUTION 


A)  Problam  ParaBCtara; 

Arch  Angle  :  90.000 
Arch  Radius:  45.000 
No  of  Design  Vac:  14 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  6 


B)  Derived  Constants: 

No  of  System  Nodal  Points...  7 

No  of  Degrees  of  Freedom....  21 

Length  per  Element .  11.7474 

Number  of  Iterations .  1 


C)  Structure  Loading: 

rx . 

FY . 

FM . 

FA . 


0.0000 

1000.0000 

0.0000 

0.0000 


D)  Elemental 

Dimensions 

and  Stress  Distribution: 

Node 

Height 

Base 

Length 

1 

3.00000 

1.50000 

11.74736 

2 

3.00000 

1.50000 

11.74736 

3 

3.00000 

1.50000 

11.74736 

4 

3.00000 

1.50000 

11.74736 

5 

3.00000 

1.50000 

11.74736 

6 

3.00000 

1.50000 

11.74736 

7 

3.00000 

1.50000 

11.74736 

Area 

4.50000 

4.50C00 

4.50000 

4.50000 

4.50000 

4.50000 

4.50000 


E)  Objective  Function: 

Total  structure  Volume:  317.178711 


Node 

Normal  Stress 

Bending  Stress 

1 

220.3 

19988.4 

2 

212.8 

19310.1 

3 

190.9 

17315.2 

4 

155.9 

14139.8 

5 

110.3 

9999.2 

6 

57.1 

5176.4 

7 

29.0 

0.2 

Total 

20208.7 
19522.9 
17506.0 

14295.7 
10109.4 

5233.4 

29.2 


F)  Boundary  Conditions: 
Node  X-Oispl 

1  1 

7  0 


Y-Dlspl 

1 

0 


Slope 

1 

0 


G)  Solution  Vector: 

Node  X-Displ 

Y-Displ 

1  O.OOOOOOE-t-00 

0.000000E4-00 

2  -0.300317E-01 

0.404077E-02 

3  -0.11208SE-*-00 

0.3811S4E-01 

4  >0.224178E-f00 

0.124215E4.00 

5  -0.336278E^00 

0.270393E4-00 

6  -0.418343E-t-00 

0.468604E4-00 

7  -0.448382B4.00 

0.696858E4-00 

Slope 

O.OOOOOOE-^00 

0.512948E-02 

0.991003E-02 

0.140157E-01 

0.171665E-01 

0.191473E-01 

0.198230E-01 
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Validation  #4 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000 
Arch  Radius:  45.000 
No  of  Design  Var:  16 


Youngs  Modulus:  30000000.0 
Yield  Strength;  52000.0 
No  of  Elements:  8 


B)  Derived  Constants: 

No  of  System  Nodal  Points...  9 

No  of  Degrees  of  Freedom....  27 

Length  per  Element .  8.8215 

Number  of  Iterations .  1 


Structure 

FX . 

FY . 

FM . . 

Loading: 

. .  0.0000 

. .  1000.0000 

.  0.0000 

FA . 

0.0000 

- 

Elemental 

Dimensions 

and  Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1.50000 

8.82154 

4.50000 

2 

3.00000 

1.50000 

8.82154 

4.50000 

3 

3.00000 

1.50000 

8.82154 

4.50000 

4 

3.00000 

1.50000 

8.82154 

4.50000 

5 

3.00000 

1.50000 

8.82154 

4.50000 

6 

3.00000 

1.50000 

8.82154 

4.50000 

7 

3.00000 

1.50000 

8.82154 

4.50000 

a 

3.00000 

1.50000 

8.82154 

4.50000 

9 

3.00000 

1.50000 

8.82154 

4.50000 

E) 


Objective  Function: 

Total  structure  Volume:  317.575592 


Node 

Normal  Stress 

Bending  Stress 

Total 

1 

221.2 

19983.1 

20204.3 

2 

217.0 

19602.4 

19819.4 

3 

204.4 

18467.7 

18672.1 

4 

184.0 

16622.6 

16806.6 

5 

156.6 

14138.1 

14294.6 

6 

123.0 

11109.3 

11232.4 

7 

84.7 

7653.0 

7737.7 

8 

43.2 

3901.8 

3944.9 

9 

21.9 

0.1 

22.0 

G)  Solution  Vector: 

Node  X-Oispl  Y-Displ 

1  O.OOOOOOE-fOO  O.OOOOOOE>00 

2  -0.170799E-01  0.17475aE-02 

3  -0.657226E-01  0.1656S6E-01 

4  -0.138527E+00  0.55548SE-01 


Slope 

O.OOOOOOE-^00 

0.388006E-02 

0.761158E-02 

0.110510E-01 


5  -0.224410E-f00  0 . 126097EfOO 

6  -0.310296E-^00  0 . 230817E-*-00 

7  -0.383113B+00  0.367110E+00 

8  >0.431768E>00  0 .  S27570E-«>00 

9  -•0.448854E'f00  0 . 701109E-f00 


0.140661E-01 

0.16S408E-01 

0.183798E-01 

0.19S124E-01 

0.198948E-01 


1  O.OOOOOOE-t-OO 

2  -0.109773E-01 

3  -0.4283868-01 

4  -0.924699E-01 

5  -0.155017E+00 

6  -0.224357E+00 

7  -0.293705E+00 

8  -0.356268E-t-00 

9  -0.405923E4-00 

10  -0,4378048*00 

11  -0.448790E*00 


O.OOOOOOE*00 

0.916237E-03 

0.861779E-02 

0.292281E-01 

0.676091E-01 

0.126884E*00 

0.208132E*00 

0.310279E+00 

0.430208E*00 

0,563056E*00 

0.702697E*00 


O.000000E*0O 

0.311363E-02 

0.615142E-02 

0.903846E-02 

0.117035E-01 

0.140808E-01 

0.161115E-01 

0.177456E-01 

0.169426E-01 

0.196729E-01 

0.199184E-01 
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validation  #4 
OPTIMIZATION  SOLUTION 


A)  Problam  Paraaeters: 

Arch  Angle  :  90.000 

Arch  Radius:  45.000 

No  of  Design  Var:  26 

h)  Derived  Constants: 

No  of  Systea  Nodal  Points... 
No  of  Degrees  of  Preedoa. . . . 

Length  per  Eleaent . 

Nuaber  of  iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Eleaents: 

12 

13 

39 

5.8863 

1 

C)  Structure  Loading: 


PX .  0.0000 

PY .  1000.0000 

PM .  0.0000 

PA .  0.0000 


Eleaental 

Oiaensio.^s 

and  Stress  1 

Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1.50000 

5.88628 

4.50000 

2 

3.00000 

1.50000 

5.88628 

4.50000 

3 

3.00000 

1.50000 

5.88628 

4.50000 

4 

3,00000 

1.50000 

5.88628 

4.50000 

5 

3.00000 

1.50000 

5.88628 

4.50000 

6 

3.00000 

1.50000 

5.88628 

4.50000 

7 

3.00000 

1.50000 

5.88628 

4.50000 

8 

3.00000 

1.50000 

5.88628 

4.50000 

9 

3.00000 

1.50000 

5.88628 

4,50000 

10 

3.00000 

1.50000 

5.88628 

4.50000 

11 

3.00000 

1.50000 

5.88628 

4.50000 

12 

3.00000 

1.50000 

5.88628 

4.50000 

13 

3.00000 

1.50000 

5.88628 

4.50000 

Objective  function; 

Total  structure  Voluae:  317.859253 
Node  Noraal  Stress  Bending  Stress 

Total 

1 

221.7 

19973.5 

20195.2 

2 

219.8 

19805.9 

20025.8 

3 

214.2 

19299.1 

19513,3 

4 

204.9 

18461,6 

18666.5 

5 

192.1 

17308.0 

17500,2 

6 

176.0 

15857.4 

16033.4 

7 

156.9 

14135.3 

14292.3 

8 

135.2 

12171.0 

12306.2 

9 

lll.l 

9998.4 

10109.5 

10 

85.0 

7653.1 

7738.1 

11 

57.4 

5176.2 

5233.7 

12 

29,1 

2610.7 

2639.8 

13 

14.7 

O.l 

14,7 
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r)  Boundary  conditions: 

Node  X-Displ  V-Displ  Slope 

11  1  1 
13  0  0  0 


G)  Solution  Vector:' 

Node  X-Oispl  y-Displ  Slope 

1  O.OOOOOOE-^OO  O.OOOOOOE-fOO  O.OOOOOOE-fOO 

2  -0.764e66E-02  0.544914E-03  0.260170E-02 

3  -0.3007S0E-01  0.504941E-02  0.51S929E-02 

4  •0.6S7S29E-01  0.172040E-01  0.762896E-02 

5  -0.112252E4-00  0.401787E-01  0.996841E-02 

6  >0.166406E'«-00  0.764064E-01  0.121375E-01 

7  -0.224523E't-00  0 . 127417E-f 00  0.140991E-01 

8  -0.282643E-)-00  0 . 19373SE-^00  0.158197E-01 

9  -0.336806E-t-00  0.2748388-^00  0.172696E-01 

10  -O.lSSllBE-i-OO  0.369200E-t-00  0.184241E-01 

11  -0.419010Ef00  0.4743878-^00  0. 1926328-01 

12  -0.441447e-»-00  0 .  S87230E-t-00  0.1977258-01 

13  -0.449100E-t-00  0 . 704036E-»'00  0.1994328-01 
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Validation  t5 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000  Youngs  Modulus:  30000000.0 


Arch  Radius:  32. 

000 

Yield  Strength: 

52000.0 

No  of  Design  Var: 

10 

No 

of  Elements: 

4 

B)  Derived  Constants: 

No  of  System  Nodal  Points 

•  •  • 

5 

No  of  Degrees  of  Freedom. 

•  •  • 

15 

Length  per 

Element. . . 

•  •  •  • 

•  •  • 

12.4858 

Number  of 

Iterations. 

.  -  .  . 

1 

C)  Structure 

Loading: 

FX . 

•  •  •  • 

*  , 

0.0000 

FY . 

•  •  •  • 

..  -5000.0000 

FM . 

•  •  «  • 

0.0000 

FA . 

.  .  . . 

0.0000 

- 

D)  Elemental 

Dimensions 

and 

Stress 

D  stribution: 

“ 

Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1.50000 

12.48578 

4.50000 

2 

3.00000 

1.50000 

12.48578 

4.50000 

3 

3.00000 

1.50000 

12.48578 

4.50000 

4 

3.00000 

1.50000 

12.48578 

4.50000 

5 

3.00000 

1.50000 

12.48578 

4.50000 

E)  Objective 
Total 


Function: 
structure  Volume: 


224.744095 


Node  Normal  Stress 

1  1231.2 

2  127S.8 

3  1273.2 

4  1073.8 

5  927.6 


Bending  Stress 
0.0 
12338.0 
11971.7 
1043.2 
24725.3 


Total 

1231.2 

1?516.9 

13245.0 

2116.9 

25652.9 


F)  Boundary  Conditions: 
Node  X-oispl 

1  1 

5  1 


Y-Oispl 

1 

0 


Slope 

0 

1 


G)  Solution  Vector: 

Node 

X-Displ 

Y-Displ 

1 

0.000000E4-00 

O.OOOOOOE-^00 

2 

-0.301621E-01 

0.S47718E-02 

3 

-0.265139E-01 

0.237552E-02 

4 

-0.530910E-02 

-0.302736E-01 

5 

o.ooooooE-*-oo 

-0.589432E-01 

Slope 

0.302544E-02 

0.131377E-02 

-0.205875E-02 

-0.357488E-02 

O.OOOOOOE-^00 
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Validation  15 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000 
Arch  Radius:  32.000 
No  of  Design  var:  14 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  6 


B)  Derived  Constants: 

No  of  System  Nodal  Points...  7 

No  of  Degrees  of  freedom....  21 

Length  per  Element .  8.3537 

Number  of  Iterations .  1 


C)  Structure  Loading: 


PX .  0.0000 

PY .  -5000.0000 

PM .  0.0000 

PA .  0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1.50000 

8.35368 

4.50000 

2 

3.00000 

1.50000 

8.35368 

4.50000 

3 

3.00000 

1.50000 

8.35368 

4.50000 

4 

3.00000 

1.50000 

8.35368 

4.50000 

5 

3.00000 

1.50000 

8.35368 

4.50000 

6 

3.00000 

1.50000 

8.35368 

4.50000 

7 

3.00000 

1.50000 

8.35368 

4.50000 

E)  Objective  Punction: 

Total  structure  Volume;  225.549301 


Node 

1 

2 

3 

4 

5 

6 
7 


Normal  Stress 
1194.9 
1247.4 

1308.2 
1279.8 

1164.2 

969.3 

853.4 


Bending  Stress 
0.0 
9411.2 

13334.9 

11503.9 

4042.8 

8539.9 
25336.8 


Total 

1194.9 

10658.6 

14643.1 

12783.7 

5207.1 

9509.2 

26240.2 


P)  Boundary  Conditions: 
Node  X-Displ 

1  1 

7  1 


y-Displ 

1 

0 


Slope 

0 

1 


G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOE-^00 

2  -0.238606E-01 

3  -0.340200E-01 

4  -0.276224E-01 

5  -0.128425E-01 

6  -0.176538E-02 

7  O.OOOOOOE-fOO 


Y-Displ 

O.OOOOOOE-^00 

0.280572E-02 

0.662210E-02 

0.12S100E-02 

-0.185793E-01 

-0.46111SE-01 

-0.61341SE-01 


Slope 

0.316688E-02 

0.229335E-02 

0.182085E-03 

-0.212342E-02 

-0.356644E-02 

-0.314903E-02 

O.OOOOOOE-fOO 
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Validation  IS 
OPTIMIZATION  SOLUTION 


A)  Pcoblam  Paramtters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
Mo  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus:  30000000.0 
Yield  strength:  52000.0 
No  of  Elements:  8 


9 

27 

6.2731 

1 


C)  Structure  Loading: 


M .  0.0000 

PY .  -5000,0000 

PM .  0.0000 

PA .  0.0000 


Elemental 

Dimensions 

and  Stress  ; 

Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1.50000 

6.27310 

4.50000 

2 

3.00000 

1.50000 

6.27310 

4.50000 

3 

3.00000 

l.SOOOO 

6.27310 

4.50000 

4 

3.00000 

1.50000 

6.27310 

4.50000 

5 

3.00000 

l.SOOOO 

6.27310 

4.50000 

6 

3.00000 

1.50000 

6.27310 

4.50000 

7 

3.00000 

1.50000 

6.27310 

4.50000 

8 

3.00000 

1.50000 

6.27310 

4.50000 

9 

3.00000 

1.50000 

6.27310 

4.50000 

B)  Objective  Function: 

Total  structure  Volume:  225.831528 


Node 

Normal  Stress 

Bending  Stress 

Total 

1 

1175.4 

0.0 

1175.4 

2 

1222.5 

7509.3 

8731.8 

3 

1292.3 

11997.3 

13289.6 

4 

1312.4 

13291.6 

14604.0 

5 

1282.1 

11342.3 

12624.5 

6 

1202.5 

6224.5 

7427.0 

7 

1076.8 

1865.4 

2942.1 

8 

909.6 

12616.3 

13525.8 

9 

816.4 

25615.1 

26431.5 

P)  Boundary  Conditions: 
Node  X-Displ 

1  1 

9  1 


Y-Displ 

1 

0 


Slope 

0 

1 


G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOB-»-00 

2  -0.190172E-01 

3  -0.314948B-01 

4  -0.342107E-01 


Y-Displ 

O.OOOOOOE4-00 

0.162605E-02 

0.513368E-02 

0.627359E-02 


Slope 

0.321682E-02 

0.269341E-02 

0.133378B-02 

-0.4288828-03 


5  -0.280319E-01 

6  -0.170219E-01 

7  -0.6540S1E-02 

8  -0.742209E-03 

9  O.O0OOOOE4-0O 


0.848441E-03 

-0.129808B-01 

-0.331005B-01 

-0.529373E-01 

-0.622147B-01 


-0.214589E-02 

-0.337032E-02 

-0.367416E-02 

-0.266477E-02 

O.OOOOOOE4-00 
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Validation  #5 
OPTIMIZATION  SOLUTION 


A)  Problaa  Paraaatars: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  oC  Design  Var:  22 

B)  Derived  Constants: 

No  of  Systea  Nodal  Points... 
No  of  Degrees  of  Preedoa.... 

Length  per  Bleaent . 

Number  of  Iterations . 

C)  Structure  Loading: 

PX . 

PY . 

PM . 

PA . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Bleaents: 

10 

11 

33 

5.0214 

1 

0.0000 

-5000.0000 

0.0000 

0.0000 

D)  Eleaental  Diaensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1.50000 

5.02138 

4.50000 

2 

3.00000 

1.50000 

5.02138 

4.50000 

3 

3.00000 

1.50000 

5.02138 

4.50000 

4 

3.00000 

1.50000 

5.02138 

4.50000 

5 

3.00000 

1.50000 

5.02138 

4.50000 

6 

3.00000 

1.50000 

5.02138 

4.50000 

7 

3.00000 

1.50000 

5.02138 

4.50000 

8 

3.00000 

1.50000 

5.02138 

4.S0000 

9 

3.00000 

1.50000 

5.02138 

4.50000 

10 

3.00000 

1.50000 

5.D2138 

4.50000 

11 

3.00000 

1.50000 

5.02138 

4.50000 

Objective 

Punction: 

Total 

structure  Voluae:  225.962219 

Node 

Normal  Stress 

Bending  Stress 

Total 

1 

1163.3 

0.0 

1163.3 

2 

1204.7 

6225.1 

7429.8 

3 

1272.0 

10546.0 

11817.9 

4 

1307.9 

12856.2 

14164.1 

5 

1311.7 

13098.8 

14410.6 

6 

1283.2 

11268.0 

12551.2 

7 

1223.1 

7408.7 

8631.8 

a 

1132.9 

1616.0 

2748.9 

9 

1014.8 

5967.4 

6982.2 

10 

871.6 

15155.0 

16026.7 

11 

794.2 

25720.4 

26514.7 

P)  Boundary  Conditions: 

Node  X-Displ  Y-Displ  Slope 

11  1  0 
11  1  0  1 
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G)  Solution  Vector: 

Nod*  X-DispI  Y-Displ 


1  O.OOOOOOE-^00 

2  -0.15655IE-01 

3  -0.277397E-01 

4  -0.339727E-01 

5  -0.337863E-01 

6  -0.282253E-01 

7  -0.195378E-01 

8  -0.105138B>01 

9  -0.366515E>02 

10  -0.355416E-03 

11  0.000000E4-00 


O.OOOOOOEt-OO 
0.103676E-02 
0.372354E-02 
0.607021E-02 
0.569723E-02  - 
0.660388E-03  - 
-0.9e3642E-02  - 
>0.249420E-01  - 
-0.419486E-01  - 
-0.564151E-01  - 
-0.62625SE-01 


Slope 

0.324004E-02 

0.289272E-02 

0.195701E-02 

0.6S1332E-03 

0.796782E-03 

0.215629E-02 

0.319832B-02 

0.370184E-02 

0.345906E-02 

0.228057B-02 

O.OOOOOOE-t-00 
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Validation  »5 
OPTIMIZATION  SOLUTION 


A)  Problaa  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Var:  26 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Clements: 

12 

13 

39 

4.18SS 

1 

C)  Structure  Loading: 


PX .  0.0000 

PY .  -5000.0000 

PM .  0.0000 

PA .  0.0000 


Elemental 

Dimensions  and 

Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1.50000 

4.18560 

4.50000 

2 

3.00000 

1.50000 

4.18580 

4.50000 

3 

3.00000 

1.50000 

4.16580 

4.50000 

4 

3.00000 

1.50000 

4.18580 

4.50000 

5 

3.00000 

1.50000 

4.18580 

4.50000 

6 

3.00000 

1.50000 

4.18580 

4.50000 

7 

3.00000 

1.50000 

4.18580 

4.50000 

8 

3.00000 

1.50000 

4.18580 

4.50000 

9 

3.00000 

1.50000 

4.18580 

4.50000 

10 

3.00000 

1.50000 

4.18580 

4.50000 

11 

3.00000 

1.50000 

4.18580 

4.50000 

12 

3.00000 

1.50000 

4.18560 

4.50000 

13 

3.00000 

1.50000 

4.18580 

4.50000 

Objective 

Function: 

Total 

structure  Volume:  226.033264 

Node 

Normal  Stress 

Bending  Stress 

Total 

1 

1155.0 

0.0 

1155.1 

2 

1191.5 

5308.7 

6500.2 

3 

1253.9 

9309.8 

10563.7 

4 

1294.8 

11935.0 

13229.8 

5 

1313.6 

13139.3 

14452.9 

6 

1309.9 

12902.2 

14212.1 

7 

1283.8 

11227.5 

12511.3 

8 

1235.7 

8144.1 

9379.9 

9 

1166.5 

3704.7 

4871.2 

10 

1077.3 

2015.0 

3092.3 

11 

969.7 

8917.1 

9886.8 

12 

845.5 

16883.3 

17728.8 

13 

779.5 

25777.3 

26556.9 
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F)  Boundary  Conditions: 

Node  X-Displ  Y-Displ  Slope 

11  1  0 
13  1  0  1 


G)  Solution  Vector: 

Node  X-Oispl  Y-Oispl  Slope 

1  O.OOOOOOC-fOO  O.OOOOOOE-fOO  0.325259E-02 

2  -0.132523E-01  0.707096E-03  0.300569E-02 

3  -0.243569E-01  0.274126E-02  0.232580E-02 

4  -0.317555E-01  0.506415E-02  0.133772E-02 

5  -0.347043E-01  0.631460E-02  0.171S46E-03 

6  -0.332895E-01  0.S14815E-02  -0 . 103961E-02 

7  -0.283308E-01  0.557844E-03  -0 . 216186E-02 

8  -0.211878E-01  -0 . 785492E-02  -0 . 306281E-02 

9  -0.134844E-01  -0 . 19«669E-01  -0 . 361389E-02 

10  -0.678007E-02  -0 . 336374E-01  -0 . 369247E-02 

11  -0.221779E-02  -0 . 475236E-01  -0 . 318404E-02 

12  -0.181929E-03  -0 . 584105E-01  -0 . 198409E-02 


13  0.000000E-«-00  -0.628492E-01  0 .  OOOOOOE-«-00 
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validation  *5 
OPTIMIZATION  SOLUTION 


A)  Problam  Paraaeters: 

Acch  Angle  :  90.000 

Arch  Radius:  32.000 

NO  of  Design  Var:  30 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Preedoa.... 

Length  per  Element . . . 

Nuoiber  of  Iterations . 

C)  Structure  Loading: 

FX . 

FY . 

FM . 

FA . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

14 

15 

45 

3.5885 

1 

0.0000 

-5000.0000 

0.0000 

0.0000 

Elemental 

Dimensions  and 

Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1.50000 

3.58851 

4.50000 

2 

3.00000 

1.50000 

3.58851 

4.50000 

3 

3.00000 

1.50000 

3.58851 

4.50000 

4 

3.00000 

1.50000 

3.58851 

4.50000 

5 

3.00000 

1.50000 

3.58851 

4.50000 

6 

3.00000 

1.50000 

3.58851 

4.50000 

7 

3.00000 

1.50000 

3.58851 

4.50000 

6 

3.00000 

1.50000 

3.58851 

4,50000 

9 

3.00000 

1.50000 

3.58851 

4.50000 

10 

3.00000 

1.50000 

3.58851 

4.50000 

11 

3.00000 

1.50000 

3.58851 

4.50000 

12 

3.00000 

1.50000 

3.58851 

4.50000 

13 

3.00000 

1.50000 

3.58851 

4.50000 

14 

3.00000 

1.50000 

3.58851 

4.50000 

15 

3.00000 

1.50000 

3.58851 

4.50000 

Objective 

Function: 

Total 

structure  Volume:  226.076035 

Node 

Normal  Stress 

Bending  Stress 

Total 

1 

1149.0 

0.0 

1149.0 

2 

1181.5 

4624.5 

5805.9 

3 

1238.8 

8296.5 

9535.3 

4 

1280.5 

10970.1 

12250.6 

5 

1306.1 

12611.5 

13917.5 

6 

1315.3 

13200.1 

14515.4 

7 

1307.9 

12728.6 

14036.5 

8 

1284.1 

11202.7 

12486.8 

9 

1244.2 

8641.9 

9886.1 

10 

1188.6 

5078.2 

6266.8 

11 

1118.0 

556.4 

1674.5 

12 

1033.5 

4866.7 

5900.2 

13 

935.9 

11123.1 

12058.9 

14 

826.5 

18133.9 

18960.3 

15 

769.0 

25810.9 

26579.9 
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F)  Boundary  Conditions: 

Node  X-Displ  y-Displ  Slope 

11  1  0 
15  1  0  1 


G)  Solution  Vector: 

Node  X-Displ  Y-Displ  Slope 

1  0 .  OOOOOOE-t-00  O.OOOOOOE-fOO  0.326008E-02 

2  -0.114679E-01  0.506369E-03  0.307569E-02 

3  -0.215484E-01  0.207183E-02  0.2S6050E-02 

4  -0.291563E-01  0.410633E-02  0.179229E-02 

5  -0. 3363S3E-01  0.579362E-02  0.852037E-03 

6  -0.347835E-01  0.624852E-02  -0 . 177131E-03 

7  -0.328342E-01  0.467242E-02  -0 . 121096E-02 

8  -0.283939E-01  0.49S743E-03  -0 . 216516E-02 

9  -0.223430E-01  -0 . 650i87E-02  -0 . 295641E-02 

10  -0.157056E-01  -0.161095E-01  -0 . 350346E-02 

11  -0.949819E-02  -0 . 276272E-01  . 372813E-02 

12  -0.456627E-02  -0 . 398712E-01  -0 , 355627E-02 

13  -0.142229E-02  -0 . 512110E-01  -0 , 291872E-02 

14  -0.946979E-04  -0 . 596559E-01  -0 . 175218E-02 


IS  Q.OQ0O0Qt*00  -0.629827E-01  0.  OOOOOOE-t-00 
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Validation  tS 
OPTIMIZATION  SOLUTIviN 


A)  Problaa  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Var:  34 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 

FX . 

FY . 

FM . 

FA . 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000,0 
No  of  Elements:  16 


17 

51 

3.1403 

1 


0.0000 

-5000.0000 

0.0000 

0.0000 


D)  Elemental 

Dimensions  and 

Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

3.00000 

1,50000 

3.14033 

4.50000 

2 

3.00000 

1.50000 

3.14033 

4.50000 

3 

3.00000 

1.50000 

3.14033 

4,50000 

4 

3.00000 

1,50000 

3.14033 

4.50000 

5 

3.00000 

1,50000 

3.14033 

4.50000 

6 

3.00000 

1.50000 

3.14033 

4.50000 

7 

3.00000 

1.50000 

3.14033 

4.50000 

e 

3.00000 

1.50000 

3.14033 

4.50000 

9 

3.00000 

1.50000 

3.14033 

4.50000 

10 

3,00000 

1.50000 

3.14033 

4.50000 

11 

3.00000 

1,50000 

3.14033 

4.50000 

12 

3.00000 

1.50000 

3.14033 

4.50000 

13 

3.00000 

1.50000 

3.14033 

4.50000 

14 

3.00000 

1.50000 

3.14033 

4.50000 

15 

3.00000 

1.50000 

3.14033 

4.50000 

16 

3.00000 

1.50000 

3.14033 

4.50000 

17 

3.00000 

1.50000 

3.14033 

4.50000 

t.;  Objective 

Function: 

Total 

structure  Volume:  226.103821 

Node 

Normal  Stress 

Bending  Stress 

Total 

1 

1144.4 

0.0 

1144.4 

2 

1173.5 

4094.6 

5268.1 

3 

1226.1 

7465.0 

8691.2 

4 

1266.9 

10078.9 

11345.9 

5 

1295.6 

11911.1 

13206.6 

6 

1311.7 

12943.9 

14255.6 

7 

1315.2 

13167.7 

14482.9 

8 

1306,0 

12579.8 

13885.7 

9 

1284.2 

11186.0 

12470.3 

10 

1250.1 

8999.8 

10250.0 

11 

1204.0 

6042.2 

7246.1 

12 

1146.2 

2341.6 

3487.8 

13 

1077.5 

2066.5 

3143.9 

14 

998.3 

7139,5 

8137.7 

15 

909.5 

12828.8 

13738.4 

16 

812.0 

19079.4 

19891,4 

17 

761.1 

25831.4 

26592.6 
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D  Boundary  Conditions: 
Node  X-Displ 

1  1 

17  1 


Y-Displ 

1 

0 


Slope 

0 

1 


G) 


Solution  Vector: 


Node 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


X-Displ 
O.OOOOOOE+00 
-0.100961E-01 
-0.192461E-01 
-0.266726E-01 
-0.3ia426E-01 
-0.344882E-01 
-0.346060E-01 
-0.324379E-01 
-0.284326E-01 
-0.231915E-01 
-0.174003E-01 
-0.117521E-01 
-0.686387E-02 
-0.319243E-02 
-0.953191E-03 
-0.472050E-Q4 
Q .OOOOOOE-OO 


Y-Displ 
O.OOOOOOE-fOO 
0.376056E-03 
0.160606E-02 
0.333145E-02 
0.503854E-02 
0.613847E-02 
0.604840E-02 
0. 426921E-02 
0.455661E-03 
-0.552506E-02 
-0,135497E-01 
-0.232130E-01 
-0.338212E-01 
-0,4440S4E-01 
-0.537569E-01 
-0.604801E-01 
-0.630648E-01 


Slope 

0.326459E-02 

0.312172E-02 

0.271838E-02 

0.210622E-02 

0.133893E-02 

0.471678E-03 

-0.439423E-03 

-0.133782E-02 

-0.216707E-02 

-0.287140E-02 

-0.339626E-02 

-0.368879E-02 

-0.369838E-02 

-0.337717E-02 

-0.268042E-02 

-0.156706E-02 

O.OOOOOOE+00 
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APPENDIX  E 
CASE  STX7DIES 


OPTIMIZATION  *1 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  0.002 

Arch  Radius:  1000000.000 
No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom. . . . 

Length  per  Element . 

Number  of  Iterations . . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

4.0000 

1 

C)  Structure  Loading: 


FX .  2000.0000 

FY .  0.0000 

FM .  0.0000 

FA .  0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length  Area 

1 

4.19530 

0.41953 

4.00000  1.76005 

2 

4.01266 

0.40127 

4.00000  1.61015 

3 

3.81169 

0.38117 

4.00000  1.45290 

4 

3.58695 

0.35869 

4.00000  1.28662 

5 

3.32982 

0.33298 

4.00000  1.10877 

6 

3.02540 

0.30254 

4.00000  0.91530 

7 

2.64292 

0.26429 

4.00000  0.69850 

8 

2.09772 

0.20977 

4.00000  0.44004 

9 

0.03041 

0.03000 

4.00000  0.00091 

Objective  Function: 

Total  structure 

Volume:  33.126362 

Node 

Normal  Stress  Bending  Stress 

Total 

1 

0.0 

52002.4 

52002.4 

2 

0.0 

52002.4 

52002.4 

3 

0.0 

52002.5 

52002.5 

4 

0.0 

52002.4 

52002.4 

5 

0.0 

52003.1 

52003.1 

6 

0.0 

52000.9 

52000.9 

7 

0.0 

52001.6 

52001.6 

8 

0.0 

51999.4 

51999.4 

9 

0.1 

844.7 

844.8 

Boundary  Conditions: 

Node 

X-Displ 

Y-Displ  Slope 

1 

1 

1 

1 

9 

0 

0 

0 

154 


G)  Solution  Vector: 

Node 

X-Oispl 

y-Displ 

1 

O.OOOOOOE+00 

O.OOOOOOE+00 

2 

0.691840E-02 

0.295492E-09  - 

3 

0.277399E-01 

0-119801E-06  - 

4 

0.6322838-01 

0.274074E-08  - 

5 

0.114358E+00 

0.496595E-08  - 

6 

0.182449E'«-00 

0.793075E-08  - 

7 

0.269490Ef00 

0.117209E-07  - 

8 

0.379305E+00 

0.165004E-07  - 

9 

0.61809SE^00 

0-268468E-07  - 

slope 

O.OOOOOOE-^00 

0.338401E-02 

0.693567E-02 

0.106947E-01 

0.147211E-01 

0.191152E-01 

0.240748B-01 

0.301571E-01 

0.744676E-01 
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OPTIMIZATION  lla 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  0.002 

Arch  Radius:  1000000.000 
No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 

C-)  Structure  Loading: 

FX . 

FY . .  . . 

FM . 

FA . 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  8 


9 

27 

4.0000 

1 


2000.0000 

0.0000 

0.0000 

0.0000 


0)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length  Area 

1 

4.19499 

0.41961 

4.00000  1.76025 

2 

4.01244 

0.40133 

4.00000  1.61030 

3 

3.81144 

0.38124 

4.00000  1.45308 

4 

3.58705 

0.35870 

4.00000  1.28669 

5 

3.32942 

0.33306 

4.00000  1.10891 

6 

3.02541 

0.30254 

4.00000  0.91531 

7 

2.64171 

0.26456 

4.00000  0.69888 

8 

2.09811 

0.20981 

4.00000  0.44021 

9 

0.10080 

0.03000 

4.00000  0.00302 

Objective  Function: 

Total  structure 

Volume:  33.148262 

Node 

Normal  Stress  Bending  Stress 

Total 

1 

0.0 

52001.2 

52001.2 

2 

0.0 

52001.2 

52001.2 

3 

0.0 

52000.0 

52000.0 

4 

0.0 

51998.6 

51998.6 

5 

0.0 

52003.4 

52003.4 

6 

0.0 

52000.0 

52000.0 

7 

0.0 

51997.4 

51997.4 

8 

0.0 

51970.5 

51970.5 

9 

0.0 

0.0 

0.0 

Boundary  Conditions 

; 

Node 

X-Displ 

Y-Oispl  slope 

1 

1 

1 

1 

9 

0 

0 

0 

G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOEfOO 

2  0.691870E-02 

3  0.277409E-01 

4  0.632300E-01 

5  0.114360E-«-00 

6  0.182452E4-00 

7  0.269496E-*-00 

8  0.379315E-*-00 

9  0.607042Ct-00 


Y-Displ 

0.000000E-*-00 

0.295505E-09 

0.119806E-08 

0.274082E-08 

0.496604E-0e 

0.793089E-06 

0.117212E-07 

0.16S008E-07 

0.263666E-07 


Slope 

O.OOOOOOE-^00 

-0.338415C-02 

-0.693591E-02 

-0.106947E-01 

-0.147212E-01 

-C.19I155E-01 

-0.240760E-01 

-0.301576E-01 

-0.703189E-01 
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OPTIMIZATION  12 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

NO  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

1 

C)  Structure  Loading: 


fX .  0.0000 

py .  -2000.0000 

FM .  0.0000 

FA .  0.0000 


Elemental 

Dimensions  and 

Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

1.95274 

1.95611 

6.27310 

3.81977 

2 

1.94555 

1.93241 

6,27310 

3.75960 

3 

1.89934 

1,90956 

6.27310 

3.62691 

4 

1.92039 

1.92025 

6.27310 

3.68763 

5 

1.92103 

1.91927 

6.27310 

3.68698 

6 

1.91825 

1,92011 

6.27310 

3.68326 

7 

1.91849 

1.92476 

6.27310 

3.69264 

8 

1.92622 

1.91786 

6.27310 

3.69422 

9 

1.96640 

1.96000 

6.27310 

3.85413 

Objective 

Function: 

Total 

structure  Volume:  186.151276 

Node 

Normal  Stress 

Bending  Stress 

Total 

1 

521.1 

51491.0 

52012.2 

2 

519.3 

51498.4 

52017.7 

3 

507.0 

51507.8 

52014.8 

4 

448.8 

45091.9 

45540.7 

5 

381,7 

38340.4 

38722,1 

6 

300.0 

30197.7 

30497.7 

7 

205.9 

20743.2 

20949.1 

8 

104.9 

10527.8 

10632.7 

9 

50.9 

0.3 

51,2 

F)  Boundary  Conditions: 
Node  X-Displ 

1  1 

9  0 


T-Displ 

1 

0 


Slope 

1 

0 


G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOC-oOO 

2  0.345895E-01 

3  0.134858C-^00 

4  0.287683E-»-00 

5  0.469409E-^00 

6  0.651576E-^00 

7  0.806243E-f00 

8  0.909657E4-00 

9  0.945957E+00 


Y-Displ 

O.OOOOOOE-^00 

-0.351714E-02 

-0.340464E-01 

-0.115848E-^00 

-0.265100E-^00 

-0.487184E-*^00 

-0.7766S9E+00 

-0.111768E+01 

-0.148635e-»-01 


Slope 

O.OOOOOOE-^00 

-0.1I0489E-01 

-0.222S60E-01 

-0.32e233E-01 

-0.419063E-01 

-0.493731E-01 

-0.549243E-01 

-0.583234E-01 

-0.594191E-01 
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OPTIMIZATION  »2a 
OPTIMIZATION  SOLUTION 


A)  problem  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

so  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom. . . . 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 

. . 

. . 

. . 

FA . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

1 

0.0000 

>2000.0000 

0.0000 

0.0000 

D)  Elemental 
Node 
1 
2 

3 

4 

5 

6 

7 

8 
9 


Dimensions  and  Stress  Distribution: 


Height 

Base 

Length 

Area 

3.91339 

0,49197 

6.27310 

1.92527 

4.10780 

0.43837 

6.27310 

1.80073 

4,11506 

0.41151 

6.27310 

1.69337 

3.97207 

0.39721 

6.27310 

1.57773 

3.69253 

0.39033 

6.27310 

1.44130 

2.42308 

0.70760 

6.27310 

1.71458 

2.31581 

0.53330 

6.27310 

1.23503 

1.30739 

0.84851 

6.27310 

1.10933 

0.80670 

1.49956 

6.27310 

1.20969 

C)  Objective  Function: 

Total  structure  Volume:  77.775108 


Node 

Normal  Stress 

Bending  Stress 

Total 

1 

1033.8 

50968.7 

52002.6 

2 

1084.1 

50916.9 

52000.9 

3 

1085.9 

50913.5 

51999.4 

4 

1048.9 

50949.6 

51998.5 

5 

976.5 

51020.7 

51997.2 

6 

645.0 

51350.6 

51995.6 

7 

616.7 

51379.1 

51995,9 

8 

350.0 

51653.2 

52003.2 

9 

162.0 

0.2 

162.2 

F)  Boundary  Conditions: 
Node  X'Displ 

Y-Displ  Slope 

1 

1 

1  1 

9 

0 

0  0 

G)  solution  Vector: 
Node  X-Displ 

1  0.000000E>00 

2  0.165765E-01 

3  0.640268E-01 

4  0.136735E-»'00 

5  0.22686lE-*-00 

6  0.324976E-:-00 

7  0.421348E+00 

8  0.499870E-f00 

9  0.534109Ef00 


Y-Displ 

O.OOOOOOE-fOO 

-0.185682E-02 

-0.164901E-01 

-0.556094E-01 

-0.129851E4-00 

-0.249654E-^00 

-0.430238E-^00 

-0.689424E4-00 

-0.103740E+01 


Slope 

O.OOOOOOEf-OO 

-0.530039E-02 

-0.104791E-01 

-0.157515E-01 

-0.213260E-01 

-0.277857E-01 

-0.369135E-01 

-0.482073E-01 

-0.595043E-01 
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OPTIMIZATION  »3 
OPTIMIZATION  SOLUTION 


A)  Pcoblem  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Var:  IS 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 

rx . 

py . 

PM . 

PA . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

1 

2000.0000 

0.0000 

0.0000 

0.0000 

D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

3.66465 

0.55100 

6.27310 

2.01923 

2 

3.83344 

0.40649 

6.27310 

1.55825 

3 

3.34200 

0.41258 

6.27310 

1.37883 

4 

2.46546 

0.54860 

6.27310 

1.35256 

5 

1.99742 

0.55572 

6.27310 

1.11001 

6 

1.48940 

0.58247 

6.27310 

0.86752 

7 

1.21374 

0.66160 

6.27310 

0.80301 

8 

0.58386 

0.48072 

6.27310 

0.28068 

9 

0.77051 

1.86386 

6.27310 

1.43612 

Objective 

function: 

Total 

structure 

Volume:  56.657707 

Nodr 

Normal  Stress  Bending  Stress 

Total 

1 

97.1 

51895.5 

51992.6 

2 

249.2 

51745.2 

51994.4 

3 

552.4 

51445.1 

51997.4 

4 

817.5 

51180.4 

51997.9 

5 

1267.9 

50730.6 

51998.5 

6 

1907.7 

50087.6 

51995.3 

7 

2290.0 

29991.8 

32281.8 

8 

6955.2 

45025.1 

51980.3 

9 

1386.0 

0.0 

1386.0 

1  Boundary 

Conditions 

; 

Node 

X- 

Displ 

Y-Oispl  Slope 

1 

1 

1 

1 

9 

0 

0 

0 

G)  Solution  Vector: 
Node  X-Displ 

1  0.000000E•^00 

2  0.ieS814E-01 

3  0.720148E-01 

4  0.15840SE-f00 

5  0.27582lE-^00 

6  0.416853E-»-00 

7  0.559677E+00 

8  0.678467E-f00 

9  0.725130E-^00 


Y-Displ 
O.OOOOOOE4.00 
-0.180717E-02 
-0.179297E-01 
-0.639459E-01 
-0.160028E-»'00 
-0.33136384-00 
-0.597638E-400 
-0.986551E4-00 
-0 . 145490E4-01 


Slope 

0.0O0000E4-00 

-0.574511E-02 

-0.117837E-01 

-0.190296E-01 

-0.286801E-01 

-0.409667E-01 

-0.537583E-01 

-0.722173E-01 

-0.764554E-01 


159 


OPTIMIZATION  i3a 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

2 

C)  Structure  Loading: 


FX .  2000.0000 

FY .  0.0000 

FM .  0.0000 

FA .  0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

4.19827 

0.41983 

6.27310 

1.76255 

2 

3.90887 

0.39089 

6.27310 

1.52793 

3 

3.58581 

0.35858 

6.27310 

1.20581 

4 

3.22371 

0.32237 

6.27310 

1.03923 

5 

2.81849 

0.28185 

6.27310 

0.79439 

6 

2.36310 

0.23631 

6.27310 

0.55843 

7 

1.84434 

0.18443 

6.27310 

0.34016 

8 

0.99509 

0.18101 

6.27310 

0.18012 

9 

0.60305 

0.06348 

6.27310 

0.03828 

Objective  Function: 

Total  structure 

Volume:  41.341122 

Node 

Normal  Stress  Bending  Stress 

Total 

1 

111.2 

51891.6 

52002.9 

2 

254.1 

51749.1 

52003.2 

3 

592.4 

51411.0 

52003.3 

4 

1064.0 

50938.9 

52002.9 

5 

1771.6 

50231.5 

52003.1 

6 

2963.6 

49040  2 

52003.8 

7 

5406.0 

46591.3 

51997.2 

8 

10838.0 

41165.8 

52003.8 

9 

51995.7 

0.0 

51995.7 

Boundary  Conditions 

Node 

X-Displ 

Y-Displ  Slope 

1 

1 

1 

1 

9 

0 

0 

0 

G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOE-i-00 

2  0.173605E-01 

3  0.677523E-01 

4  0.146066E-^00 

5  0.251532E-I-00 

6  0.367672E>00 

7  0.482192Ef00 

8  0.577701E+00 

9  0.627437B+00 


Y-Displ 

O.OOOOOOE4-0O 

-0.168480E-02 

-0.168e06E-01 

-0.596161E-01 

-0.1441S3E-)-00 

-0.284910E>00 

-0.497395E+00 

-0.806931Ef00 

-0.126845E+01 


Slope 

O.OOOOOOEfOO 

-0.536752E-02 

-0.111584B-01 

-0.175022E-01 

-0.246090E>01 

-0.328323E-01 

-0.428618E-01 

-0.575020E-01 

-0.822389E-01 
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OPTIMIZATION  #4 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants; 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom. . . . 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

2 

C)  Structure  Loading: 


rx .  0.0000 

ry .  -2000.0000 

PM .  1000.0000 

FA .  0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution; 


Node 

Height 

Base 

Length 

Area 

1 

4.20406 

0.42041 

6,27310 

1.76741 

2 

4.17639 

0.41764 

6,27310 

1.74422 

3 

4.09209 

0.40921 

6.27310 

1.67452 

4 

3.94749 

0.39475 

6.27310 

1.55827 

5 

3.73436 

0.37344 

6.27310 

1.39455 

6 

3,43725 

0.34373 

6.27310 

1.18147 

7 

3.02033 

0. 302^3 

6.27310 

0.91224 

8 

1.95373 

0.35108 

6.27310 

0.68592 

9 

0.79249 

0.21196 

6.27310 

0.16798 

Objective  Function: 

Total  structure 

Volume:  63.252686 

Node 

Normal  Stress  Bending  Stress 

Total 

1 

1126.2 

50868.9 

51995.1 

2 

1119.2 

50874.8 

51994.1 

3 

1098.2 

50896.4 

51994.6 

4 

1062.0 

50929,1 

51991.1 

5 

1009.3 

50986.3 

51995.6 

6 

936.0 

51055.6 

51991.6 

7 

835.0 

51156.8 

51991.8 

8 

566.1 

51424.8 

51990.9 

9 

1166.5 

45071.7 

46238.2 

1  Boundary  Conditions 

* 

Node 

X-Displ 

Y-Displ 

Slope 

1 

1 

1 

1 

9 

0 

0 

0 

G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOE-^OO 

2  0.158777E-01 

3  0.619007E-01 

4  0.133298E1-00 

5  0.222312E-»^00 

6  0.3ie467E-«-00 

7  0.408886E-f00 

8  0.479103Ef00 

9  0. 512886E-K00 


Y-Dlspl 

0.0O000OE-)>0O 

-0.180201E-02 

-0.160076E-01 

-0.54428eE-01 

-0.127764E4-00 

-0.24S2S6E4-00 

-0.414819E4-00 

-0.848809E-»-00 

-0.990889E4.00 


Slope 

0.000000E4-00 

-0.507744E-02 

-0.102264E-01 

-0.155283E-0I 

-0.2I0892E-01 

-0.270698E-01 

-0.337708E-01 

-0.425065E-01 

-0.605568E-01 
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OPTIMIZATION  «4a 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Vac:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 


youngs  Modulus: 

30000000.0 

yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

2 

C)  Structure  Loading: 


PX .  0.0000 

py .  -2000.0000 

PM .  10000.0000 

PA .  0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

3.99646 

0,39965 

6.27310 

1.59717 

2 

3.96579 

0.39658 

6.27310 

1.57275 

3 

3.87202 

0.38720 

6.27310 

1.49925 

4 

3.70940 

0.37094 

6.27310 

1.37596 

5 

3.46539 

0.34654 

6.27310 

1.20089 

6 

3.11238 

0.31124 

6.27310 

0.96869 

7 

2.57625 

0.25763 

6.27310 

0.66371 

8 

1.43829 

0.14383 

6.27310 

0.20687 

9 

2.25936 

0.22769 

6.27310 

0.51444 

Objective 

Pune t ion: 

Total 

structure 

Volume:  53.206200 

Node 

Normal  Stress  Bending  Stress 

Total 

1 

1246.2 

50757.2 

52003.4 

2 

1241.2 

50761.6 

52002.8 

3 

1226.5 

50776.2 

52002.8 

4 

1202.8 

50799.4 

52002.1 

5 

1172.0 

50828.7 

52000.7 

6 

1141.6 

50859.8 

52001.4 

7 

1147.6 

50851.9 

51999.5 

8 

1877.1 

50127.0 

52004.0 

9 

381.0 

51622.1 

52003.1 

Boundary 

Conditions 

; 

Node 

X- 

Displ 

y-Displ  Slope 

1 

1 

1 

1 

9 

0 

0 

0 

G)  Solution  Vector: 


Node 

1 

2 

3 

4 

5 

6 

7 

8 
9 


X-Displ 

O.OOOOOOE-^00 

0.166829E-01 

0.650796E-01 

0.140297E-^00 

0.234385E-I-00 

0.336635Ef00 

0.434064E-^00 

0.51S013E-^00 

0.544307E>00 


y-Displ 

O.OOOOO0E4-OO 

-0.190698E-02 

-0.168603E-01 

-0.573560E-01 

-0.134897E+00 

-0.2S9876E-«>00 

-0.442870E>00 

-0.710561E>00 

-0.100920E-^01 


Slope 

O.OOOOOOE4.00 

-0.533238E-02 

-0.107525E-01 

-0.163633E-01 

-0.223074E-01 

-0.288281E-01 

-0.365050E-01 

-0.496258E-01 

-0.415951E-01 
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OPTIMIZATION  15 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters; 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 

FX . 

FY . 

FM . 

FA . 


Youngs  Modulus: 

30000000.0 

Yield  Strength; 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

2 

0.0000 

0.0000 

0.0000 

-100.0000 

D)  Elemental 

Dimensions 

and  Stress 

Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

4.94233 

0.49484 

6.27310 

2.44566 

2 

4.59883 

0.46036 

6.27310 

2.11711 

3 

4.20711 

0.42104 

6.27310 

1.77135 

4 

3.76835 

0.37698 

6.27310 

1.42060 

5 

3.27617 

0.32787 

6.27310 

1.07416 

6 

2.72309 

0.27231 

6.27310 

0.74152 

7 

2.08293 

0.20953 

6.27310 

0.43644 

8 

1.31418 

0.13235 

6.27310 

0.17394 

9 

0.14293 

0.03000 

6.27310 

0.00429 

E)  Objective 

Function: 

Total  ucture  Volume:  55.704273 


Node 

N  >1  Stress 

Bending  Stress 

Total 

1 

.173.9 

50828.3 

52002.2 

2 

1210.8 

50791.8 

52002.6 

3 

1109.9 

50893.1 

52002.9 

4 

996.4 

51006.4 

52002.7 

5 

868.5 

51135.7 

52004.2 

6 

724.0 

51279.7 

52003.7 

7 

555.5 

51446.8 

52002.3 

8 

351.8 

51646.5 

51998.3 

9 

1.8 

4.3 

6.6 

F)  Boundary  Conditions: 


Node 

X-Displ 

Y-Displ 

1 

1 

1 

9 

0 

0 

Slope 

1 

0 


G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOE+00 

2  0.1443S3E-01 

3  0.564283B-01 

4  0.1236S7E-f00 

5  0.210848C4-00 

6  0.309650E-^00 

7  0.408546E-^00 

8  0.492877E-f00 

9  0.557215E>00 


Y-Displ 

O.OOOOOOE^OO 

-0.168651C-02 

-0.146790E-01 

-0.508635E-01 

-0.122671E-»00 

-0.243321E1-00 

-0.428622E-^00 

-0.706927E*0n 

-0.136017Ef0: 


Slope 

O.OOOOOOE-^00 

-0.447172E-02 

-0.932896E-02 

-0.147199E-01 

-0.208707E-01 

-0.2ei898E-01 

-0.375918E-01 

~0.522128E-01 

-0.130850E+00 
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OPTZHZZATZOS  *S 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

NO  of  Design  Var:  id 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  pec  Clement . 

Number  of  Iterations . 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Clements:  6 


9 

27 

6.2731 

1 


C)  Structure  Loading: 


FX .  0.0000 

FY .  -2000.0000 

FK .  0.0000 

FA .  -100.0000 


D)  elemental  Dimensions  and  Stress  Distribution: 


Mode 

Height 

Base 

Length 

Area 

1 

5.82316 

0.58232 

6.27310 

3.39091 

2 

5.56534 

0,55653 

6.27310 

3.09731 

3 

5.25381 

0.52538 

6.27310 

2,76025 

4 

4.88828 

0.48883 

6.27310 

2.38953 

5 

4.46201 

0.44620 

6.27310 

1.99096 

6 

3.96199 

0.39620 

6.27310 

1.56974 

7 

3,35911 

0,33591 

6.27310 

1.12836 

8 

2.27093 

0.32719 

6.27310 

0.74303 

9 

0,40176 

0.51850 

6.27310 

0.20831 

Objective  Function; 

, 

Total  structure 

Volume;  97,474487 

Node 

Normal  Stress  Bending  Stress 

Total 

1 

1433.6 

50561.4 

51995.1 

2 

1457.9 

50537.3 

51995.2 

3 

1378.4 

50617.0 

51995.5 

4 

1284.9 

50710.9 

51995.8 

5 

1175.4 

50821.4 

51996.9 

6 

1046.4 

50951.4 

51997.8 

7 

889.9 

51109.2 

51999.1 

8 

605.0 

51393.6 

51998.6 

9 

940.6 

0.6 

941.2 

1  Boundary  Conditions 

Node 

X-Displ 

Y-Oispl 

Slope 

1 

1 

1 

1 

9 

0 

0 

0 

G)  Solution  Vector: 
Node  X-Displ 

1  o.ooooooe-^00 

2  0.1ie390B-01 

3  0.461727C-01 

4  0.100311C-^00 

5  0.169098e>00 

6  0.244936B4-00 

7  0.317850000 

8  O.37Se26B-v00 

9  0.404874C-»-00 


Y-Dlspl 

O.OOOOOOC-t-OO 

-0.148107C-02 

-0,122074B-01 

-0.414621C-01 

-0.982482C-01 

-0.191026000 

-0.327870C+00 

-O.S19533C-fOO 

-0.815207C-»-00 


Slope 

O.OOOOOOE+00 

-0.371830C-02 

-0.763781C-02 

-0.1ie3l9C-0l 

-0.164007E-01 

-0.215055E-01 

-0.274519E-01 

-0.353806E-01 

-0.533505E-01 
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OPTIMIZATION  «6a 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90 . 000 

Arch  Radius:  32.000 

No  of  Design  var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . . 

Number  of  Iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

1 

C)  Structure  Loading: 


fX .  0.0000 

fV .  -2000.0000 

.  0.0000 

FA .  -100.0000 


Elemental 

Dimensions 

and  Stress 

Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

5.82323 

0.58232 

6.27310 

3.39101 

2 

5.56545 

0.55654 

6.27310 

3.09742 

3 

5.25394 

0.52539 

6.27310 

2.76039 

4 

4.88843 

0.48884 

6.27310 

2.38968 

5 

4.46221 

0.44622 

6.27310 

1.99114 

6 

3.96197 

0.39623 

6.27310 

1.56986 

7 

2.96078 

0.43150 

6.27310 

1.27759 

8 

2.01471 

0.41514 

6.27310 

0.83639 

9 

0.53831 

1.20792 

6.27310 

0.65023 

E)  Objective  Function: 


Total  structure  Volume:  101.764938 


Node 

Normal  Stress 

Bending  Stress 

Total 

1 

1433.6 

50556.2 

51989.7 

2 

1457.8 

50531.4 

51989.2 

3 

1378.3 

50609.9 

51988.2 

4 

1284.9 

50703.3 

51988.1 

5 

1175.3 

50812.4 

51987.8 

6 

1046.3 

50946.8 

51993.1 

7 

786.1 

51211.2 

51997.3 

8 

537.6 

51463.1 

52000.6 

9 

301.4 

0.5 

302.0 

F)  Boundary  Conditions: 
Node  X-Displ 

1  1 

9  0 


y-Dispi 

1 

0 


Slope 

1 

0 


G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOE-^00 

2  0.118375E-01 

3  0.461666E-01 

4  0.100297E4.00 

5  0.169072E-^00 

6  0.244899E-f00 

7  0.318251E+00 

8  0.377833E+00 

9  0.404775Ef00 


Y-Displ 

O.OOOOOOE4.00 

-0.148090E-02 

-0.122058E-01 

-0.414S62E-01 

-0.982332E-01 

-0.190995E-*-00 

-0.328633E-^00 

-0.525529E+00 

-0.799478E+00 


Slope 

O.OOOOOOEt-00 

-0.371783E-02 

-0.763674E-02 

-0.118301E-01 

-0.163980E-01 

-0.215019E-01 

-0.277224E-01 

-0.367195E-01 

-0.474619E-01 
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OPTIMIZATION  17 


OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000  Youngs  Modulus:  30000000.0 

Arch  Radius:  32.000  Yield  Strength:  52000.0 

No  of  Design  Var:  18  No  of  Elements:  8 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 

No  of  Degrees  of  Freedom.... 

Length  per  Clement . 

Number  of  Iterations . 

C)  Structure  Loading: 


FX .  0.0000 

FY .  -8000.0000 

FM .  0.0000 

FA .  0.0000 


Elemental 

Dimensions  and 

Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

0.49766 

0.32385 

6.27310 

0.16117 

2 

2.49460 

0.47039 

6.27310 

1.17344 

3 

3.05960 

0.46175 

6.27310 

1.41278 

4 

2.96505 

0.48983 

6.27310 

1.45238 

5 

2.41709 

0.49892 

6.27310 

1.20593 

6 

0.96211 

0.15996 

6.27310 

0.15390 

7 

3.08578 

0.43109 

6.27310 

1.33023 

8 

4.00244 

0.54499 

6.27310 

2.18129 

9 

4.21890 

0.79543 

6.27310 

3.35583 

Objective 

Function : 

Total 

structure  Volume:  64.558678 

Node 

Normal  Stress 

Bending  Stress 

Total 

1 

51998.6 

0.2 

51998.8 

2 

7361,8 

44637.4 

51999.2 

3 

6359.0 

45640.2 

51999.3 

4 

6185.7 

45813.3 

51999.0 

5 

7163.7 

44836.1 

51999.7 

6 

51734.2 

267.2 

52001.4 

7 

5246.1 

46199.3 

51445.4 

8 

2625.6 

49374.5 

52000.1 

9 

1501.8 

50497.5 

51999.4 

F)  Boundary  Conditions: 

Node  X-Displ  Y-Displ  Slope 

11  1  0 
9  1  0  1 


G)  Solution  Vector: 

Node  X-Displ  Y-Oispl  Slope 

1  O.OOOOOOE+OO  O.OOOOOOC+00  0.307764E-01 

2  -0.149671C+00  0.117776E-01  0.102322E-01 

3  -0.192268C-^00  0.232035E-01  0.336248C-02 

4  -0.193966E+00  0.225997C-01  -0 , 298150E-02 

5  -0. 161968E+00  -0 . 546921E-02  -0 . 100998E-01 

6  -0.785744E-01  -0 . 112048E+00  -0 .273030E-01 

7  -0.161085C-01  -0.234507E+00  -0 . 111216E-01 

a  -0.133252E-02  -0 . 285889E+00  -0 . 514754E-02 
9  O.OOOOOOC^OO  -0.303321C-t-00  0 .  OOOOOOE-^OO 
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OPTIMIZATION  #7a 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters; 

Arch  Angle  :  90.000  Youngs  Modulus:  30000000.0 

Arch  Radius:  32.000  Yield  Strength;  52000.0 

No  o£  Design  Var:  18  No  of  Elements:  8 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 

No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 


FX .  0.0000 

FY .  -8000.0000 

FM .  0.0000 

FA .  0.0000 


Elemental 

Dimensions  and 

Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

0.67370 

0.23920 

6.27310 

0.16115 

2 

2.61798 

0.42995 

6.27310 

1.12561 

3 

3.27094 

0.40726 

6.27310 

1.33212 

4 

3.15268 

0.43628 

6.27310 

1 . 37546 

5 

2.37931 

0.51326 

6.27310 

1.22121 

6 

0.77194 

0.19791 

6.27310 

0.15277 

7 

3.03132 

0.44320 

6.27310 

1.34348 

8 

4.32486 

0.46885 

6.27310 

2.02772 

9 

4.62704 

0.66330 

6.27310 

3.06910 

Objective 

Function : 

Total 

structure  Volume:  61.674786 

Node 

Normal  Stress 

Bending  Stress 

Total 

1 

52005.8 

0 . 1 

52005.8 

2 

7674.6 

44327.7 

52002.3 

3 

6743.9 

45258.6 

52002.5 

4 

6531.3 

45470.1 

52001.4 

5 

7073.5 

44927.7 

52001.1 

6 

52110.1 

1128 . 4 

53238.5 

7 

5193.7 

46613.6 

51807.3 

8 

2824.0 

49177.3 

52001.2 

9 

1641.8 

50360.2 

52002.1 

F)  Boundary  Conditions; 

Node  X-Displ  Y-Displ  Slope 

11  1  0 
9  1  0  1 


G)  Solution  Vector: 

Node  X-Displ  Y-Displ  Slope 

1  O.OOOOOOE+00  O.OOOOOOEt-00  0.285876E-01 

2  -0. 140675E+00  0.1065:5E-01  0.102757E-01 

3  -0.184847Ef00  0.224797E-01  0.386062E-02 

4  -9.190542E+00  0.239361E-01  -0 . 204050E-02 

5  -0.163888000  0.229189E-03  -0 .8e6111E-02 

6  -0.790007E-01  -0 . 1081 36E+00  -0 . 284664E-01 

7  -0.148789E-01  -0 . 233588E+00  -0 . 104461E-01 

8  -0.115786E-02  -0 . 281574Et-00  -0 . 47225  3E-02 


9  O.OOOOOOE+OO  -0 . 297573E«-00  0  .  OOOOOOEt-00 
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OPTIMIZATION  «8 
OPTIMIZATION  SOLUTION 


A)  Problsm  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Nusdier  of  Iterations . 

C)  Structure  Loading: 

rx . 

FY . 

FM . 

FA . 


Youngs  Modulus: 

3000000C .0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

6.2731 

2 

0.0000 

-8000.0000 

0.0000 

0.0000 

Elemental  Dimensions 

and  Stress  : 

Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

3.96647 

0.39665 

6.27310 

1.57329 

2 

2.66054 

0.26605 

6.27310 

0.70784 

3 

2.69799 

0.26980 

6.27310 

0.72791 

4 

3.29169 

0.32917 

6.27310 

1.08352 

5 

3.29281 

0.32928 

6.27310 

1.08426 

S 

2.70298 

0.27030 

6.27310 

0.73061 

7 

2.65240 

0.26524 

6.27310 

0.70352 

8 

3.95112 

0.39511 

6.27310 

1.56114 

9 

4.87425 

0.48694 

6.27310 

2.37347 

Objective  Function: 

Total  structure  Volume:  52 

.992058 

Node 

Normal  Stress 

Bending 

Stress 

Total 

1 

5469.7 

46526.9 

51996.6 

2 

12833.6 

39153.1 

51986.7 

3 

13542.6 

38456.9 

51999.5 

4 

9462.4 

42532.7 

51995.1 

5 

9456.6 

42533.2 

51989.9 

6 

13496.0 

38480.4 

51976.4 

7 

12918.1 

39108.3 

52026.4 

8 

5103.2 

46906.0 

52009.2 

9 

3085.5 

48816.0 

51901.5 

Boundary  Conditions: 

Node 

X-Displ 

Y-Displ 

Slope 

1 

1 

1 

1 

9 

1 

0 

1 

Solution  Vector: 

Node 

X-Displ 

Y-Displ 

Slope 

1 

0.000000E>00  0. 

OOOOOOE-^00 

O.OOOOOOE-^00 

2 

-0.237838E-01  0. 

695637E-03 

0.631554E-02 

3 

-0.683462E-01  0. 

11302SE-01 

0.624288E-02 

4 

-0.894876E-01  0. 

199ie9E-01 

0.336519E-03 

5 

-0.792890E-01  0. 

896548E-02 

-0.506634E-02 

6 

-0.479138E-01  -0. 

329894E-01 

-0.109648E-01 

7 

-0.147824E-01  -0. 

100892E>00 

-0.110874E-01 

8 

-0.821552E-03  -0. 

152604E-»00 

-0.4712e4E-02 

9 

O.OOOOOOEt-OO  -0. 

168974E-^00 

O.OOOOOOE-cOO 
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OPTIMIZATION  #8a 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  90.000 

Arch  Radius:  32.000 

Arch  Height:  2.000 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus:  30000000.0 
Yield  Strengths  52000.0 
NO  of  Elements:  12 


13 

39 

4.1858 

1 


C)  Structure  Loading: 

. . 

. . 

. . 

FA . 


0.0000 

-8000.0000 

0.0000 

0.0000 


Elemental  Dimensions 


Element 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Height 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 


and  Stress  Distribution: 


Base 

Length 

1.72707 

4.18580 

0.94577 

4.18580 

0.27194 

4.18580 

0.65040 

4.18580 

1.09075 

4.18580 

0.78986 

4.18580 

1.06223 

4.18580 

0.75388 

4.18580 

0.25849 

4.18580 

0.93275 

4.18580 

1,75263 

4.18580 

2,65143 

4.18580 

E)  Objective  Functi 

Total 

structure 

Node 

Stress 

1 

51991.36 

2 

49986.21 

3 

46231.32 

4 

51991.33 

5 

49404.85 

6 

51991.04 

7 

51841,79 

6 

42080.21 

9 

51997.34 

10 

51279,07 

11 

52434.15 

12 

51908,94 

13 

51935.26 

F)  Boundary  Condi t; 

Node 

X-Displ 

1 

1 

13 

1 

107.886574 


Y-Oispl 

1 

0 


Slope 


Volume 

14.45838 

7.91762 

2.27661 

5.44488 

9.13130 

6.61236 

8.89255 

6.31116 

2.16400 

7.80861 

14.67236 

22.19675 


169 


G) 


Solution  Vector: 

Node  X-Oispl  Y-Displ  Slope 

1  O.OOOOOOE-^00  O.OOOOOOE-t-OO  0 .  OOOOOOE-«-00 

2  -0.121854E-01  0.457452E-03  0.527930E-02 

3  -0.442765E-01  0.615Q66E-02  0.947556E-02 

4  -0.861533E-01  0.177060E-01  0.907372E-02 

5  -0.112986E-^00  0.297132E-01  0.435047E-02 

6  -0.121058Ef00  0.342979E-01  -0.109393E-03 

7  -0.110454E4-00  0.237526E-01  -0 . 70ia4eB-02 

8  -0.850803E-01  -0 . 622361E-02  -0 . 11S516E-01 

9  -0.537900E-01  -0 .547404E-01  -0 . 155000E -01 

10  -0.252856E-01  -0 . 118485E-«-00  >0 . 155501E-01 

11  -0.744947E-02  -0 . 173190E+00  -0 . 110976B-01 

12  -0.654979E-03  -0 . 209063E+00  -0. 580765E-02 

13  O.OOOOOOEt-OO  -0.222069E-»-00  0 . OOOOOOE-t-00 
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OPTIMIZATION  #9 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  180.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom. . . . 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 

rx . 

FY . 

FM . 

FA . 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  8 


9 

27 

12.4858 

2 


16000.0000 

0.0000 

0.0000 

0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

3.87125 

2.61830 

12.48578 

10.13610 

2 

2.99788 

1.94284 

12.48578 

5.82439 

3 

0.78829 

0.32181 

12.48578 

0.25368 

4 

2.77540 

1.34286 

12.48578 

3.72635 

5 

2.82891 

1.32326 

12.48578 

3.7S337 

6 

2.47737 

1.02727 

12.48578 

2.54492 

7 

1.60172 

1.17975 

12.48578 

1.88963 

8 

0.98586 

0.84212 

12.48578 

0.83021 

9 

1.27732 

1.14827 

12.48578 

1.46670 

E)  Objective  Function: 

Total  structure  Volume:  287.147583 


Node 

Normal  Stress 

Bending  Stress 

Total 

1 

573.5 

51426.4 

51999.9 

2 

1458.0 

50542.0 

52000.0 

3 

51241.7 

770.7 

52012.4 

4 

4166.3 

47833.6 

52000.0 

5 

2239.2 

49761.4 

52000.5 

6 

404.8 

51595.6 

52000.5 

7 

1007.1 

50991.3 

51998.4 

8 

2995.1 

49002.5 

51997.6 

9 

1835.2 

0.1 

1835.2 

F)  Boundary  Conditions: 

Node 

X-Displ  Y- 

Displ 

Slope 

1 

1 

1 

1 

9 

0 

1 

0 

G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOE-^00 

2  0.904933E-01 

3  0.558262E>00 

4  0.883004E-^00 

5  0.909118E-K00 

6  0.923784E4-00 

7  0.109988E-C01 

8  0.169034E-^01 

9  0.266363E-^01 


Y-Displ 

O.OOOOOOE-fOO 

-0.176853E-01 

-0.327630E-*-00 

-0.806151E•^00 

-0.928170E-^00 

>0.854074E*00 

-0.590018E>00 

-0.194612E4.00 

O.OOOOOOE4.0O 


Slope 

O.OOOOOOE400 

-0.130647E-01 

-0.608843E-01 

-0.171659E-01 

-0.266783E-02 

0.134924B-01 

0.348143E-01 

0.717773E-01 

0.833529E-01 
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OPTIMIZATION  t9a 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  180.000 
Arch  Radius:  32.000 
No  of  Design  Var:  18 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  8 


B)  Derived  Constants: 

No  of  System  Modal  Points...  9 

No  of  Degrees  of  Freedom. ...  27 

Length  per  Element .  12.4858 

Niaiber  of  Iterations .  2 


C)  Structure  Loading: 


FX .  16000.0000 

FY .  0.0000 

FM .  0.0000 

FA .  0.0000 


Elemental 

Dimensions 

and  Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

1.07555 

0.19099 

12.48578 

0.20542 

2 

5.49600 

0.69864 

12.48578 

3.83971 

3 

5.93999 

0.96376 

12.48578 

5.72472 

4 

5.98177 

1.04035 

12.48578 

6.22316 

5 

6.00000 

0.54928 

12.48578 

3.29571 

6 

5.57273 

0.58244 

12.48578 

3.24580 

7 

4.24805 

0.49093 

12.48578 

2.08549 

8 

2.17557 

0.52420 

12.48578 

1.14044 

9 

0.24362 

0.66203 

12.48578 

0.16128 

Objective 

Function: 

Total 

structure 

Volume:  344.338989 

Node  Normal  Stress  Bending  Stress 

Total 

1 

53383.0 

1.3 

53384.3 

2 

3451.4 

50161.1 

53612.5 

3 

2907.2 

50644.1 

53551.3 

4 

2812.1 

50764.5 

53576.5 

S 

2854.3 

77667.4 

80521.7 

6 

925.1 

52414.6 

53339.7 

7 

2660.2 

50774.2 

53434.4 

8 

6355.6 

47117.6 

53473.2 

9 

48642.3 

0.0 

48642.3 

Boundary 

Conditions : 

Node  X> 

Displ 

Y-Displ 

Slope 

1 

1 

1 

0 

9 

0 

1 

0 

Q)  Solution  Vector: 

Node 

X-Displ 

Y-Displ 

1 

O.OOOOOOE-i-OO 

O.OOOOOOEt-OO 

2 

0.538764E-^00 

-0.103983E-^00 

3 

0.766928E-«-00 

-0.254801E-^00 

4 

0.867834E4.00 

-0.403592E>00 

5 

0.88S116B-^00 

-0.482743E-«-00 

6 

0.891319E-^00 

-0.450538E4.00 

7 

0.97461lEi-00 

-0.324620E4.00 

8 

0.122814E-f01 

-0.153178E4.00 

9 

0.197489Ei-01 

0.000000E4.00 

Slope 

-0.532524E-01 

•0.253326E-01 

-0.178851E-01 

-0.108023E-01 

-0.246003E-02 

0.696965E-02 

0.161255E-01 

0.301513E-01 

0.765027E-01 
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OPTIMIZATION  110 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  180.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

NusUier  of  Iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

12.4858 

2 

C)  Structure  Loading: 


rx .  16000.0000 

FY .  0.0000 

FM .  0.0000 

FA .  0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

3.87807 

2.57526 

12.48578 

9.98702 

2 

3.14748 

1.72693 

12.48578 

5.43547 

3 

0.88930 

0.29131 

12.48578 

0.25906 

4 

2.92467 

1.14515 

12.48578 

3.34919 

5 

2.96833 

1.05559 

12.48578 

3.13335 

6 

2.03171 

1.03518 

12.48578 

2.10318 

7 

0.77820 

0.60572 

12.48578 

0.47137 

8 

1.64419 

0.87326 

12.48578 

1.43580 

9 

1.68701 

1.15400 

12.48578 

1.94680 

E)  Objective  Function: 

Total  structure  Volume:  256.608276 


Node 

Normal  Stress 

Bending  Stress 

Total 

1 

631.7 

51368.8 

52000.5 

2 

1646.5 

50354.2 

52000.7 

3 

51530.6 

479.7 

52010.3 

4 

4692.9 

47307.6 

52000.5 

5 

2706.4 

49294.1 

52000.5 

6 

579.8 

51421.7 

52001.5 

7 

4780.6 

47222.8 

52003.4 

8 

2050.8 

49951.2 

52002.0 

9 

1637.1 

50364.7 

52001.8 

F)  Boundary  Conditions: 
Node  X-Displ 

1  1 

9  0 


Y-Displ 

1 

1 


Slope 

1 

1 


G)  Solution  Vector: 

Node  X-Displ  Y-Displ 

1  O.OOOOOOE-f-00  O.OOOOOOE*00 

2  0.882631E-01  -0 . 172024E-01 


Slope 

O.OOOOOOE-f-00 

-0.127257E-01 


3 

4 

5 

6 

7 

8 
9 


0.520760E-f00 

0.82313lE-f00 

0.e50524E<»'00 

0.86S174E-»-00 

0.115273E-I-01 

0.165293Ef01 

0.181555E4-01 


-0.303339E-«-00 

-0.747606E+00 

-0.874578E1-00 

-0.800410E-^00 

-0.368882E4-00 

-0.331472E-01 

O.OOOOOOE^OO 


-0.S59400E-01 

-0.172292E-01 

-0.359896E-02 

0.136811E-01 

0.570417E-01 

0.251769E-01 

O.OOOOOOE-^00 
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OPTiniZATION  *11 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  180.000 
Arch  Radius:  32.000 
No  of  Design  Var:  18 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  8 


B)  Derived  Constants: 

No  of  System  Nodal  Points...  9 

No  of  Degrees  of  Freedom....  27 

Length  per  Element .  12.4658 

Number  of  Iterations .  2 


C)  Structure  Loading: 


FX .  0.0000 

FY .  -12000.0000 

FM .  1000.0000 

FA .  0.0000 


D)  Elemental  Dimensions 

and  Stress  Distribution: 

Node 

Height 

Base 

Length  Area 

1 

1.15624 

1.02865 

12.48578  1.19142 

2 

2.00490 

1.12219 

12.48578  2.24987 

3 

2.04000 

1.08818 

12.48578  2.21989 

4 

0.57320 

0.21096 

12.48578  0.12092 

5 

2.42160 

1.36445 

12.48578  3.30415 

6 

1.23817 

0.70722 

12.48578  0.87587 

7 

i.S43:;5 

1.08778 

12.48578  1.67882 

8 

1.43388 

0.95977 

12.48578  1.37619 

9 

1.49848 

1.26079 

12.48578  1.68926 

E)  Objective  Function: 

Total  structure  Volume:  153.076752 

Node 

Normal  Stress 

Bending  Stress  Total 

1 

5944.7 

0.0 

5944.7 

2 

3277.7 

48720.7 

51998.4 

3 

3321.9 

48527.9 

51849.8 

4 

51700.0 

253,5 

51953.4 

5 

1618.5 

50381.3 

51999.8 

6 

6817,0 

45180.4 

51997.4 

7 

4083.4 

47917.5 

52000.9 

8 

4865.8 

47133.8 

51999.6 

9 

3360.4 

48642.0 

52002.3 

F)  Boundary  Conditions: 

Node 

X-Displ 

Y-Displ 

Slope 

1 

1 

1 

0 

9 

1 

1 

1 

G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOEfOO 

2  -0.471806E-fOO 

3  -0.619634E-f00 

4  -0.359468E-:>00 

5  -0.269622E-^00 
$  -0.213488E4-00 

7  0.128494E-01 

8  0.B88700E-01 

9  O.OOOOOOE-^00 


Y-Displ 

O.OOOOOOE4-00 

0.920ei0E-01 

0.189140E-f00 

-0.20647eE-^00 

-0.667966E-^00 

-0.379826E-^00 

-0.370866E-01 

0.160224E-01 

0.000000E4-00 


Slope 

0.458664E-01 

0.241675E-01 

0.415893E-02 

-0.589599E-01 

0.495288E-02 

0.343967E-01 

0.214467E-01 

-0.531602E-02 

O.OOOOOOE-^00 


174 


OPTIMIZATION  #lla 
OPTIMIZATION  SOLUTION 


1 


A)  Problem  Parameters: 

Arch  Angle  :  180.000 

Arch  Radius:  32.000 

No  of  Design  var:  18 

B)  Derived  Constants; 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 

rx . 

FY . 

FM . 

FA . 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  8 


9 

27 

12.4858 

2 


0.0000 

-24000.0000 

1000.0000 

0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

1.25459 

1.11170 

12.48578 

1.39472 

2 

2.56424 

1.39293 

12.48578 

3.57180 

3 

2.61780 

1.33869 

12.48578 

3.50443 

4 

0.76364 

0.31775 

12.48578 

0.24265 

5 

3.03972 

1.74821 

12.48578 

5.31405 

6 

1.59975 

0.89444 

12.48578 

1.43088 

7 

1.98431 

1.34883 

12.48578 

2.67651 

8 

1.85250 

1.19860 

12.48578 

2.22040 

9 

1.96692 

1.46917 

12.48578 

2.88973 

Objective  Function: 

Total  structure 

Volume;  241.778809 

Node 

Normal  Stress  Bending  Stress 

Total 

1 

10133.6 

0.0 

10133.7 

2 

4120.0 

47880.2 

52000.2 

3 

4199.2 

47797.9 

51997.1 

4 

51413.5 

611.5 

52025.0 

5 

2009.1 

49990.9 

52000.0 

6 

8340.7 

43657.1 

51997.9 

7 

5124.9 

46875.8 

52000.8 

8 

6039.8 

45960.0 

51999.8 

9 

4402.2 

47599.3 

52001.5 

Boundary  Conditions 

Node 

X-Displ 

Y-Displ 

Slope 

1 

1 

1 

0 

9 

1 

1 

1 

*• 


G)  Solution  Vector: 

Node  X-Displ  Y-Displ  Slope 

1  O.OOOOOOEt-00  O.OOOOOOEfOO  0.377959E-01 

2  -0.377866E-»'00  0.726540E-01  0.168604E-01 


3  -0.474151E+00  0.134826E+00  0.149644E-02 

4  -0.277207E-*-00  -0 . 167482E-*-00  -0 . 440847E-01 

5  -0.210637E+00  -0 . 513900E+00  0.340668E-02 

6  -0.168729E+00  -0 . 295879E+00  0.2622B7E-01 

7  0.359865E-02  -0 . 330045E-01  0.166200E-01 

8  0.647848E-01  0.107673E-01  -0 . 365121E-02 

9  O.OOOOOOE+00  0  .OOOOOOE-^00  0 . OOOOOOE-^00 
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OPTIMIZATION  «12 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  180.000 

Arch  Radius:  32.000 

No  o£  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

12.4858 

2 

C)  Structure  Loading: 


PX .  -9000.0000 

py .  -17000.0000 

PM .  1000.0000 

PA .  0.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

2.07490 

1.64551 

12.48578 

3.41427 

2 

0.81408 

0.49069 

12.48578 

0.39946 

3 

2.06842 

0.85360 

12.48578 

1.76560 

4 

2.08744 

0.85702 

12.48578 

1.78898 

5 

1.11662 

0.11166 

12.48578 

0.12468 

6 

2.82593 

0.84646 

12.48578 

2.39204 

7 

3.62188 

1.18265 

12.48578 

4.28341 

8 

0.81741 

0.25531 

12.48578 

0.20869 

9 

2.18518 

1.96045 

12.48578 

4.28395 

Objective  Function: 

Total  structure 

Volume:  156.554611 

Node 

Normal  Stress  Bending  Stress 

Total 

1 

1331.2 

50668.2 

51999.4 

2 

13159.1 

38839.6 

51998.8 

3 

3525. 9 

48473.0 

51998.9 

4 

3491.6 

48507.9 

51999.4 

5 

42640.2 

6458.5 

49098.6 

6 

1495.5 

50503.6 

51999.1 

7 

1368.2 

50631.9 

52000.2 

8 

51935.8 

45.6 

51981.4 

9 

2914.9 

49085.0 

52000.0 

Boundary  Conditions 

; 

Node 

X-Displ 

y-Dispi 

Slope 

1 

1 

1 

1 

9 

1 

1 

1 

G)  Solution  Vector: 
Node  X-Displ 

1  O.OOOOOOEt-00 

2  -0.385756E-^00 


Y-Dlspl  Slope 

O.OOOOOOE-^00  O.OOOOOOE-cOO 
0.7S4815E-01  0.480384E-01 


3 

4 

5 

6 

7 

8 
9 


-0.777271E+00 

-0.808788E-f00 

-0.737499E+00 

-0.812760E-^00 

-0.800801E-^00 

-0.415861E-*'00 

O.OOOOOOEfOO 


0.334000E>00 

0.378437E>00 

0.351551E-02 

-0.364503E-f00 

-0.345997e-^00 

-0.859055E-01 

0.000000E-^00 


0.140370E-01 

-0.S38832E-02 

-0.430591E-01 

-0.423870E-02 

0.959411E-02 

0.510140E-01 

O.OOOOOOE-i-00 
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OPTIMIZATION  #13 
OPTIMIZATION  SOLUTION 


A)  Problsm  Pacaneters; 

Arch  Angle  :  180.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 

C)  Structure  Loading: 

FX . 

FY . 

FM . 

FA . 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  8 


9 

27 

12.4658 

2 


16000.0000 

0.0000 

1000.0000 

0.0000 


Elemental 

Dimensions 

and  Stress  Distribution: 

Node 

Height 

Base 

Length  Area 

1 

3.88833 

2.61929 

12.48578  10.18466 

2 

3.06710 

1.88559 

12.40578  5.78328 

3 

0.92232 

0.26590 

12.48570  0.24524 

4 

2.98689 

1.21410 

12.48578  3.62638 

5 

3.00467 

1.32694 

12.48578  3.98700 

6 

1.88382 

0.97124 

12.40578  1.02964 

7 

0.67362 

0.55665 

12.48578  0.37497 

8 

1.49367 

0.78025 

12.48578  1.16544 

9 

1.67967 

1.18504 

12.48578  1.99048 

Objective 

Function: 

Total 

structure  Volume:  265.960205 

Node  Normal  Stress 

Bending  Stress  Total 

1 

536.3 

51465.6 

52001.9 

2 

1412.3 

50589.8 

52002.0 

3 

51994.3 

96.5 

52090.8 

4 

4244.8 

47757.2 

52002.1 

5 

3935.9 

48065.6 

52001.5 

6 

4523.6 

47478.1 

52001.8 

7 

4413.4 

47572.7 

51986.1 

8 

1855.2 

50144.8 

52000.0 

9 

1175.7 

36538.5 

37714.2 

Boundary 

Conditions : 

Node  X- 

Displ 

Y-Displ 

Slope 

1 

1 

1 

1 

9 

0 

1 

1 

G)  Solution  Vector: 

Node  X-Displ  Y-Displ  Slope 

1  O.OOOOOOEfOO  O.OOOOOOE4-00  O.OOOoOOE-^00 

2  0.892425E-01  -0 . 174S56E-01  -0 . 128954E-01 


3  0.527041E+00  -0 . 307447E+00  -0 . 566398E-01 

4  0.833359E+00  -0 . 758307E+00  -0 . 176816E-01 

5  0.862231E+00  -0.894399E+00  -0 . 436349E-02 

6  0.880363E>00  >0 . 814806E4-00  0.140516E-01 

7  0.117817E+01  -0.368085E+00  0.590398E-01 

8  0.168336Ef01  -0 . 291622E-01  0.222087E-01 

9  0 . 182676E4-01  0  .  OOOOOOEf-OO  0  .  OOOOOOEfOO 
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OPTIMIZATION  il3a 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  180.000 

Arch  Radius:  32.000 

No  of  Design  Var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points... 
No  of  Degrees  of  Freedom.... 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus: 

30000000.0 

Yield  Strength: 

52000.0 

No  of  Elements: 

8 

9 

27 

12.4858 

2 

C)  Structure  Loading: 


FX .  16000.0000 

FY .  0.0000 

FM .  1000.0000 

FA .  0.0000 


Elemental 

Dimensions 

and  Stress  Distribution: 

Node 

Height 

Base 

Length 

Area 

1 

5.58951 

1.27357 

12.48578 

7.11864 

2 

5.10116 

0.69418 

12.48578 

3.54110 

3 

1.56517 

0.15652 

12.48578 

0.24497 

4 

4.38397 

0.58589 

12.48578 

2.56854 

5 

4,63376 

0.58228 

12.48578 

2.69813 

6 

3.27838 

0.34568 

12.48578 

1.13328 

7 

1.61894 

0.16189 

12.48578 

0.26210 

8 

2,57960 

0.25796 

12.48578 

0.66544 

9 

2,91537 

0.29154 

12.48578 

0.84994 

Objective 

Function: 

Total 

structure  Volume:  175.647415 

Node  Normal  Stress 

Bending  Stress 

Total 

1 

765.7 

51234.1 

51999.8 

2 

2303.6 

49696.1 

51999.8 

3 

52019.3 

211.4 

52230.6 

4 

5991.5 

46008.2 

51999.7 

5 

5816.0 

46183.7 

51999.7 

6 

7301.1 

44698.3 

51999.4 

7 

6284.3 

35672.9 

41957.2 

8 

3234.0 

48765.9 

51999.9 

9 

2740.5 

47790.4 

50530.9 

Boundary 

Conditions : 

Node  X- 

Displ 

Y-Displ 

Slope 

1 

1 

1 

1 

9 

0 

1 

1 

Solution  Vector: 

Node 

X-Displ 

Y-Displ 

1 

O.OOOOOOE-f-OO 

O.OOOOOOE^OO 

2 

0.S63288E-01 

-0.107647E-01 

3 

0.306688E-^00 

-0.174214E4-00 

4 

0.481161B-:>00 

-0.42S410E4.00 

5 

0.497296E-f00 

-0.493441E4.00 

6 

0.513184E4-00 

-0.431268E4.00 

7 

0.666860E>00 

-0.199686E4.00 

8 

0.932019E4-00 

~0.202692E-01 

9 

0.102734Ef01 

0.000000E4.00 

Slope 

O.OOOOOOEfOO 

-0.813259E-02 

-0.318503E-01 

-0.977707E-02 

-0.126135E-02 

0.955349E-02 

0.294S84E-01 

0.147632E-01 

O.OOOOOOE+OO 
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OPTIMIZATION  #14 
OPTIMIZATION  SOLUTION 


A)  Probltm  Parameters: 

Arch  Angle  :  180.000 

Arch  Radius:  32.000 

No  of  Design  var:  18 

B)  Derived  Constants: 

No  of  System  Nodal  Points. 
No  of  Degrees  of  Freedom. . 

Length  per  Element . 

Number  of  Iterations . 


Youngs  Modulus; 
Yield  Strength: 
No  of  Elements: 


9 

27 

12.4858 

2 


30000000.0 

52000.0 

8 


C)  Structure  Loading: 

rx . 

FY . 

FM . 

FA . 


9000.0000 

-5000.0000 

1000.0000 

0.0000 


'  Elemental 

Dimensions 

and  Stress  Distribution: 

Node 

Height 

Base 

Length  Area 

1 

2.36812 

1.89649 

12.48578  4.49111 

2 

0.84905 

0.08490 

12.48578  0.07209 

3 

3.71044 

0.71004 

12.48578  2.63458 

4 

3.05720 

0.50168 

12.48578  1.53373 

5 

1.56758 

0.15986 

12.48578  0.25059 

6 

2.11368 

0.60867 

12.48578  1.28654 

7 

2.28870 

0.73028 

12.48578  1.67139 

8 

2.09919 

0.65663 

12.48578  1.37838 

9 

0.85662 

0.28601 

12.48578  0.24500 

Objective 

Function: 

Total 

structure  Volume:  121.283012 

Node  Normal  Stress 

Bending  Stress  Total 

1 

251.3 

48722.1 

48973.4 

2 

18883.4 

33152.8 

52036.2 

3 

422.4 

51577.1 

51999.5 

4 

831.1 

51168.1 

51999.1 

5 

9637.3 

42361.6 

51998.9 

6 

2478.0 

49521.0 

51996.9 

7 

2079.4 

49920.1 

51999.5 

8 

2346.2 

49653.4 

51999.6 

9 

12165.1 

0.2 

12165.3 

Boundary 

Conditions: 

Node  X- 

Displ 

Y-Displ  Slope 

1 

1 

1 

1 

9 

1 

1 

0 

G)  Solution  Vector: 
Node  X'Displ 

1  0 .  OOOOOOE't-00 

2  0.427757E*00 

3  0.8193e3E-K00 

4  0.822885E4-00 

5  0.767984E+00 

6  0.028458Ef00 

7  0.S41128E^00 

8  0 .653208E-K00 

9  O.OOOOOOE-^OO 


Y-Displ 

O.OOOOOOE^OO 

-0.853866E-01 

-0.3461B3E-)>00 

-0.351655E>00 

-0.809962E-01 

O.ZBISOBE-^OO 

0.2S2532E400 

0.12811SC^00 

O.OOOOOOE-^00 


Slope 

O.OOOOOOE-vOO 

-0.52S055E-01 

-0.786212E-02 

0.527751E-02 

0.3138598-01 

0.108871E-01 

-0.809329E-02 

-0.271114E-01 

-0.664133B-01 
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OPTIMIZATION  *15a 
OPTIMIZATION  SOLUTION 


A)  Problem  Parameters: 

Arch  Angle  :  180.000 
Arch  Radius:  32.000 
No  of  Design  Var:  18 


Youngs  Modulus:  30000000.0 
Yield  Strength:  52000.0 
No  of  Elements:  8 


B)  Derived  Constants: 

No  of  System  Nodal  Points...  9 

No  of  Degrees  of  Freedom. ...  27 

Length  per  Element .  12.4858 

Number  of  Iterations .  2 


C)  Structure  Loading: 


rx, .  0.0000 

FY .  32000.0000 

FM .  1000.0000 

FA .  100.0000 


D)  Elemental  Dimensions  and  Stress  Distribution: 


Node 

Height 

Base 

Length 

Area 

1 

3.25280 

2.04372 

12.48578 

6.64780 

2 

3.20982 

1.57865 

12.48578 

5.06717 

3 

1.05585 

0.53193 

12.48578 

0.56164 

4 

3.45172 

1.93749 

12.48578 

6.68769 

5 

4.53938 

2.40127 

12.48578 

10.90026 

6 

3.87991 

2.04928 

12.48578 

7.95104 

7 

2.87337 

1.82061 

12.48578 

5.23130 

8 

1.78269 

1.33703 

12.48576 

2.38351 

9 

1.43268 

1.44987 

12.48578 

2.07577 

E)  Objective  Function: 

Total  structure  Volume:  516.579224 


Node 

Normal  Stress 

Bending  Stress 

Total 

1 

3235.7 

48783.9 

52019.6 

2 

3969.0 

48027.9 

51996.9 

3 

28716.2 

2630.5 

31348.8 

4 

1520.8 

50476.2 

51997.0 

5 

574.6 

51424.1 

51998.7 

6 

1021.3 

50976.6 

51997.8 

7 

2359.0 

49638.0 

51997.0 

8 

6361.0 

45636.9 

51998. C 

9 

7781.0 

0.1 

7781.1 

F)  Boundary  Conditions: 

Node 

X-Displ 

Y-Displ 

Slope 

1 

1 

1 

1 

9 

0 

1 

0 

G)  Solution  Vector: 

Node  X-Displ  Y-Displ 

1  O.OOOOOOE+OO  O.OOOOOOE-^OO 

2  -0.800741E-01  0.174874E-01 


Slope 

O.OOOOOOE+00 

0.125050C-01 


3  -0.427959E+00 

4  -0.65460aE-»-00 

5  >0.666896E-^00 

6  -0.6e2870E>00 

7  -0.811784E+00 

8  -0.118582E-)-01 

9  -0.183804E-*-01 


0.254098E-^00 

0.S96950E>00 

0.660397E-I-00 

0.578787E+00 

0.384649E-»-00 

0.132794E^00 

0.000000E4-00 


0.443218E-01 

0.100603E-01 

-0.108817E-02 

-0.114177E-01 

-0.243721E-01 

-0.440123E-01 

-0.579555E-01 
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